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Characteristics and Sound Speed in 
Nonisentropic Gas Flows With 
Nonequilibrium Thermodynamic States' 


E. L. RESLER, JR.* 


Cornell University 


SUMMARY 


A theory is developed extending the method of characteristics 
to flow problems involving lagging heat capacities and/or finite 
chemical reaction rates. Special cases appropriate to present 
high-speed flow problems are considered in more detail—namely, 
heat-capacity lag and dissociation of a gas made up of diatomic 
molecules. The theory is easily adaptable to the usual finite 
wave approximation used for computations 


INTRODUCTION 


Ww MUCH ATTENTION being focused on problems 
relating to both chemical and gas-dynamic 
phenomena simultaneously, it is appropriate to examine 
how chemical processes affect flow parameters or, 
similarly, how flow phenomena can be utilized in chem1- 
cal-kinetic studies. For example, the mechanism of 
excitation of the vibrational degrees of freedom of air 
molecules is of importance in studies of the internal 
aerodynamics of air-breathing jet engines; similar 
problems face the rocket engine designer. The lag 
times of heat capacities, including processes of dissocia- 
tion and ionization, affect the flow about high-speed 
missiles. Also, aerodynamic phenomena have been 
utilized to study basic chemical kinetics of some reac- 
tions—e.g., the impact-tube method of measuring re 
laxation times * or the shock tube used for various 
studies. * This paper develops a theory which treats 
flows that are not isentropic due to phenomena occur- 
ring in the gas itself, such as a net energy flow between 
various degrees of freedom of molecules due to heat- 
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capacity lag or chemical reactions. Emphasis is placed 
on a formulation free of particular assumptions concern- 
ing the chemical reaction rates. Previous theories, of 
course, occur as special cases. The theory is developed 
for one-dimensional unsteady and _ two-dimensional 
steady flows with special attention paid to processes 


of dissociation and ionization. 


ONE-DIMENSIONAL UNSTEADY FLOWS 


Let & be the single space coordinate on which the 
velocity u (in the ~ direction) and the thermodynamic 
variables depend, and let ¢ be the time. Then if p is 
the density, p the pressure, and .1 the area of cross 
sections normal to the coordinate &, the fluid flow equa- 
tions, neglecting friction and other transport phe- 


nomena, are 
Equation of Motion: 
(Ou/Ot) + u(Ou/dOE) = (1/p)(Op/O8) 
Conservation of Mass: 
A(Op/Ot) + (Opu.1/0E) = O 


Although we have only two equations and three un 
knowns, we will proceed as far as possible without 
specifying what our third equation is. If we can com- 
bine Eqs. (1) and (2) so as to define the derivative of 
our unknowns along a particular direction, then that 
direction will be a ‘‘characteristic direction’”® if it is a 
function of &, ¢, and the unknowns. Define the operator 
6/ét as 


6/é6t = (O/Ot) + (u + a)(O/OE 


that is, the derivative in the é, ¢ plane in the direction 
dt/dt = u + a. The velocity a is as yet undefined but 
will turn out to be the speed of sound. Choose wu as the 
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unknown to operate on; then we have 
6u/5t = (Ou/Ot) + (u + a)(Ou/dOE) = 
(Ou/Ot) + u(Ou/OE) + a(Ou/OE) 
Using Eq. (1) and (2) to eliminate the u derivatives, 
one finds 
6u/6t = —(1/p)(Op/0&) + 
a[—(1/p)(Op/Ot) — (ufp)(Op/OE) — (u/.1)(0.1/0€) | 
but since A = A(é) only, 
OA /OE = [1/(u + a) |(64 /ét) 
and defining the substantial derivative D/Dt as 
D/Dt = (0/Ot) + u(O/0E) 
we find 
a(O0/0é) = +(6/6t) * (D/Dt) 
therefore, 
6u/dt = *(1/pa)(6p/dt) F [ua/.l(u + a) |(6.1/6t) + 
[(1/pa)(Dp/Dt) — (a/p)(Dp/Dt) | 
Now if we set the last bracket of this expression 
equal to zero, we have defined a and have a differential 
expression for our unknowns in a particular direction; 
so for 
a? = (Dp/Dt)/(Dp/Dt) (3) 
bu = F(1/pa)ép F [ua/A(u + a) |b (4) 


The quantity a so defined can now definitely be identi- 
fied as the speed of small disturbances. For a small dis- 


turbance theory with 64 = 0, and replacing the 6 by 
differences, we have from Eq. (4), Ap = *¥ podyAu 
along directions d§/dt = +d), where 


ay = [(Op/dt)/(Op/dF) |y' ” 


and the subscript 0 denotes that the quantities are to 
take on the values appropriate to the undisturbed 


state. Then in the direction di/dt = dy) the quantity 
P = p + podou is propagated unchanged and, in the 
direction d&/dt = —a), Q = p — podou is constant. 
Notice that 


p= (P+ Q)/2 and u = (P — Q)/2pod 


Consider a disturbance moving into a region where Q 
= constant; then as P is propagated unchanged with 
velocity ad, so are u and p and, therefore, the disturb- 
ance. Eq. (3) which has now been identified as the 
speed of sound was given by M. Munk.’ As Munk 
points out, the speed of sound depends upon the rate of 
change of p with p along a particle path. This defini- 
tion is more general than the usual more specialized 
relation a? = dp/dp |isentropie- 

Also note that for 64 = 0 we have 6p = ¥ paéu. 
That we actually have defined the velocity of wave 
propagation can also be recognized from this expression. 
Across the wave we have a change in pressure 6p, which 
is the force changing the velocity by amount 6u of an 


amount of fluid per unit time given by pa, the amoy 
of fluid affected by the wave per unit time. 

Suppose we now introduce the gas law into our syst, 
of equations in the form 


p= pRT 


Here 7’ is the temperature appropriate to the trang, 
tional degrees of freedom of the gas atoms or molecule 
This temperature may or may not be the same as th 
temperature appropriate to describe any of the othe 
degrees of freedom the molecule may possess. We wi 
no longer concern ourselves with the area term in E 
(4), as we may take it into account at any time; so fy 
6A = 0, 


bu = *(6pRT/pa) = *(6RT/a) * (RT/a)é6 In p 


When chemical reactions occur in the flow field the 
affect the velocity u only insofar as they change th 
pressure p. This is done by changing either th 
translational temperature T or the gas constant R by 
means of the reaction. To keep account of energy e 
change between the translational degrees of freedon 
and other degrees, it is convenient to define the entrop 
n of that part of the system that is always in equ: 





librium with the translational temperature 7. The de-| 


grees of freedom not in equilibrium with 7 we will call| 


the lagging degrees of freedom; these are to be con-| 


sidered as a heat reservoir carried along with each unit 
mass of the flow, which can change the entropy of th 
unit mass by energy exchange with that part of the in 
ternal energy in equilibrium with the temperature / 
How much heat is exchanged between the lagging an 
following degrees of freedom is governed by the laws o/ 


chemical kinetics and depends on the thermodynamic | 


state of the unit mass as it proceeds through the flow 
field. Of course, if the chemical composition of the unit 
mass changes, this must also be taken into account 
In defining 7 in this way we have assumed that there 
always exists a translational temperature 7° which is de 
fined by means of the gas law p = pRT. 

If E is the internal energy that is in equilibrium at all 
times with the translational temperature 7°, we cat 
write, with 


E = K(T) 
din p = (dE/RT) — (1/R)dn ‘ 


To compute 7 we must follow the particles in their m 
tion, in which case if Ul’ is the lagging internal energ) 


we can write 
Dyn/Dt = —(1/T)/(DU/Dt) a 


The specification of DU/Dt amounts to specifying 
the needed chemical relations. So, given DU/Dt, the 
specification of our problem is complete. Unfortunately 
in most cases of aerodynamic interest the form of the 
function specifying DU’/Dt has not been determined. 

Summarizing then, we have, in the most general case 
the following set of equations to solve: 
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+ [6( + RT)/a| + 
(T'/a)én * [ua/A(u + a) 6A 
(Dp/Dt)/(Dp/Dt) 


a 
p pRT 
E E(T) 
din p (1/RT)dE — (1/R)dn 
Dn/Dt —(1/T)(DU/Dt) (10) 
DU/Dt sum total of appropriate chemical 
rates of components with lag- 
ging internal energy 
R = function of the concentrations 


of the components making up 
the fluid 
Perfect Gas, Constant Composition, Lagging Heat Capacity 
Consider the general case where l’ = >~U;, where U; 
1 
is the energy in lagging part 7 of the internal energy. 
Then 


Dn/Dt = —(1/T)>(DU,/Dt) 


and if EK = K(T), 


RT)(dE/dT)(DT/Dt) + 
(1/RT)>((DU;,/Dt) 


Din p/Dt = (1 


To find a”, we use the gas law; viz., 


2 = (Dp/Dt)/(Dp/Dt) 


(p/p)|(D In p/Dt)/(D In p/Dt) | 
p (D1n p/Dt) + (D In T/Dt) 
- p (D In p/Dt) 
= RT X 


1/R)(dE/dT)(DT/Dt) + (DT/Dt) + (1/R)>(DU,/Dt) 


(1/R)(dE/dT)(DT/Dt) + (1/R) (DU ;/Dt) 


Now let C:, denote dE/dT; then if U; = U;,(0;) where 
0, is a temperature typifying L’;,, and if C»; denotes 
dU,/d0,, the formula for a? becomes 


(Cop + R)(DT/Dt) + YOCo(DO;/Dt) 


a? = RT (11) 


Co>(DT/Dt) + >3(De,/Dt) 


Also, 


bu = (12) 


(Cop + R)/alé7T + (T/a)éin if 6.1 = 0 


0; here defined need not be interpreted as a temperature 
in the usual thermodynamic sense, but merely as a 
means of accounting for the energy in the “reservoir” 
It is 
convenient, however, in many cases to identify it in such 


associated with the lagging internal energy lL’; 


a way that when 0, becomes equal to 7, LU’; takes on its 


equilibrium value. Such a procedure is followed below. 


Lagging Rotational Heat Capacity of a Diatomic Molecule 

In this case Co, = (3/2)R and Co; = R, at least at 
temperatures such that the rotational heat capacity is 
fully excited. Let Org be the rotational temperature. 





NONISENTROPIC 











GAS 





FLOWS 


Then expression (11) for the speed of sound becomes 


_ A DT/Dt) + 2(De0r/Dt) = 
an = BrRT 
3(DT/Dt) + 2(D0r/Dt) 


a* = (13) 

Note that if D6r/Dt = Owe get the monatomic gas 
result—i.e., a? = (5/3)R7T—and if Org = T so that 
DT/Dt = D®pr/Dt, we get the diatomic gas result—i.e., 
5)RT. 


we get appropriate intermediate sound speeds. 


a= (7 For 9p lagging 7 by various amounts 


Working with expression (12), we find 


6u = *(5/2)(R/a)éT + (T/a)éin 
+ 6u = —(5/2)(R/a)é(a*/BrR) + (T/a)in 
5[(5a/Br) + u| = (5/Br)ba + (5a/2BR*)bBr — 


(5/Br)6a — (5a/Br*) Br + (a/RBR) pn 


6[(5a Br) + u | = 


(a/Br)|—(5/2)6 In Be + (On/R)] (14) 


Using Eqs. (13) and (14), and given that D6,/Dt = 
f(p, 7, Or), one can proceed as in the usual case to con- 
struct solutions. Notice that if D@p/Dt = O, then 
Br = 5/3 = constant and 6n = 0, so that 6(3a + u) = 0, 
in agreement with the Riemann invariants for the case 
of a perfect monatomic gas with constant heat capaci- 


ties. If Og = T everywhere, then Dn = —(R/T)DT or 
n = —RIinT + m along particle paths, where mp is the 
entropy, say at ¢ = O, of the particles. If all the par- 


ticles were initially in the same thermodynamic state, 
then m is the same for all particles, and we can use the 
relation 7 = —RI1n 7 + m anywhere in the field. Thus, 
for Or = i Br = 7 S. 6Br = 0, and 
(a/RBr)én = [a/R(7/5)|[—R(6T/T)] = 
—(5a/7a*)2a6a = —(10/7)éa 

so that Eq. (14) becomes 

6} [5a/(7/5)] + u} = —(10/7)é6a 
or 6(5a + u) = 0 
in accordance with well known theory for a perfect 
diatomic gas with constant heat capacities. 


Lagging Vibrational Heat Capacity of a Diatomic Molecule 


This is similar to the previous case. Suppose the ro- 
to 


equilibrium with the translational degrees; then if Us; = 


tational degrees of freedom are assumed be in 
(vibrational energy in degrees of freedom 7), one finds, 


since & = (5/2)RT, 


7(DT/Dt) + (2/R)>-(DU>,/Dt) 


a? = R7 =BRIT (15) 


5(DT/Dt) + (2/R)>3(DU>;/Dt) 


6[(7a/B,) + u| = 


(a/B,)|—(7/2)6 In B, + (6n/R)] (16) 


It is appropriate for some applications to neglect the 
lag of rotation and not that of vibration, even though 
one knows that in reality all the degrees of freedom lag 


behind the translational degrees of freedom. One must 
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decide first what mass of fluid he is willing to regard as 
the smallest mass of interest to him in any particular 
circumstance, then find the typical volume and thus 
length, from the chosen mass and a typical density. 
Finding the time, given by this length divided by the 
speed of sound, one can neglect the adjustment times of 
degrees of freedom which adjust faster than this time. 
Of course, if the particular circumstances in question 
are such that the phenomena cover a large range of 
density, the criteria as to when times are negligible or 
not will change with the density and temperature. For 
example, in ordinary acoustics we are quite content to 
assume that air is a diatomic gas with no lagging rota- 
tional heat capacity, which we know is not true. This 
is quite a good approximation, since the human ear 
does not respond to frequencies above 20,000 cycles 
per sec., or to a time scale of about 50 usec. The rota- 
tional degrees of freedom come to equilibrium with the 
translational degrees of freedom in at most 100 colli- 
sions. Under atmospheric conditions the time between 
collisions is about 10~° sec., so that the time scale for 
the adjustment for rotational degrees of freedom is 10 


DU/Dt = q(Da/Dt) + (D/Dt)(E,) = gqo(Da/Dt) + (dE,/d8,)(D0,/Dt) = q(Da/Dt) + C 


where C,,;, denotes the vibrational heat capacity, which may be a function of 6,. 


( 
a? = (1 + a)ReT ~ 
= BRT 


5) [(7 + 3a)/Bla + “ui = 


AERONAUTICAL 


(5 + a)(DT/Dt) + T(Da/Dt)} (2q0/RoT) — 1} 


SCIENCES NOVEMBER, 1957 

sec. Thus the human ear has a frequency response 5) 
times too low to be able to detect any rotational hea 
capacity lag effects in acoustic waves. 


Dissociating Diatomic Gas With Vibrational Heat-Capacit, 
¢ 
In this case we must allow for a chemical change j; 
species, due to the dissociation as the particles mo, 
about. If a@ is the fraction of the original diatom; 
molecules that have become dissociated, we have fo 
the gas constant R = (1 + a)Ro, where Ry is the ga 
constant appropriate to the undissociated particles 
We shall assume that the rotational degrees of freedoy 
of the diatomic molecules are always in equilibrium wit! 
the translational degrees of freedom. Then, per uni 
mass, we can write the energy 


— a)RoT + 
(3/2)(2a)RoT = [(5 + a@)/2)Ri 


B= (6/201 


and U = awt E,(O,) 
where q is the dissociation energy, /¢, the energy in th 
vibrational heat capacity (only one mode for a diatomi 
gas), and 0, is the vibrational temperature. Thus, 
,(DO,/Dt 


Proceeding as before, 


7 + 3a)(DT/Dt) + T(Da/Dt)}1 + (2q0/RoT)} + (2Csin/Ro)(DO,/Dt) 


+ (2C,i»/Ro)(DO,/Dt 17 


(3a/2B)6a — |[(7 + 3a)/B7]a68 + (a/BR)én 


= (a/8)\(3/2)6a — [(7 + 3a)/B]6B + [5n/(1 + a)R,y|! 


Equations similar to Eqs. (17) and (1S) can be derived 
in like fashion for the case of ionization of a monatomic 


gas. 
STEADY Two-DIMENSIONAL FLOW 


To discuss steady two-dimensional flow, we will choose 
as coordinates the distance along a streamline, s, and 
the normal to the streamline, ». Then our hydrody- 


namic equations are 


Continuity : 
Od /On = —(1/p)(Op/Os) — (1/w)(Ow/Os) (19) 
Motion: 
0d0/Os = —(1/pw?)(Op/dn), 
w(Ow/Os) = —(1/p)(Op/ds) (20) 


where w is the velocity of the fluid and ¢ the angle that 
the velocity vector makes with some reference direction. 
Define the directional derivative 6/6/ as 6/6/ = (O0/Os) + 
a/V Me — 1)(0/On), where M = w/a. 
is not yet determined, as a is not yet defined. Operating 
on the variable ¢, we have 


The direction 


53/81 = (08/ds) = (1/V M2 — 1)(d8/dn) 


Using Eqs. (19) and (20), 
63/5/ = —(1/pw?)(Op/On) + (1/V M? — 1)(1/p) X 
(Op/Os) * (1/V M? — 1pw?)(dp/d 


and using the identity 





d0/dn = VM? — 1[¥(6/dl) + (0/ds)| 
we have 
63/6] = +(V M2? — 1/pw?)(5p/dl) + 
(VM? — 1 pw)(Op/Os) + (1 VM? — 1) X 
(1/p)(Op/Os) ¥* (1 VM? — 1 pw*)(Op/ds) (2 


Setting the terms with the streamline derivatives equ 


to zero, we find 


a” = (Op/0s)/(Op/Os) 2 





) 


in agreement with Eq. (3), since derivatives along th 


particle path. 
to show that a so defined is the speed of propagation 
Also, 


streamline direction are now derivatives along tht] 
Again, as in the unsteady case, it is eas! | 


small disturbances —1.e., the speed of sound. 


that 


But 


so th 


sponse 5) ; 
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NONISENTROPIC 


8/6l = +(V AM? — 1/pw*)(6p/5l (23) 
l'sing the gas law, Eq. (23) becomes 
1/6 V M? — 1/pw*)(6pRT/5/) = 
£(V M? — 1/w?) X 
6RT/6l) + (VM 1 RT/w*)[{(6 In p)/él 
id using Eq. (S 
1/5 V Me? — 1/w*)[6(E + RT) /6/] 
VM 1 /w?)T(6n/6/ 24) 
In this case, 
On/Os -{1/T)(OU/oOs (25 


Os amounts to a description 


of the chemical kinetics. 


Summarizing the two-dimensional relations, we have 


the following set of equations to work with: 


F) V MVP? — 1/w°)6(E + RT) F 
(V M? — 1/w?)T bn 
Op, Os) /(Op/Os 
pRT 
E(T 
In p 1/RT)dE — (1/R dn 
dn/0s 1/T)(Ol/ds (26 
OU’/ds sum total of appropriate chemical 
rates of components with lagging 
internal energy 
| wea 
R function of the concentrations of the 


components making up the fluid 
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at any point along a stream tube the stagnation e1 
thalpy will be constant. Then to the set (26) we can 
add the equation 

h=E R1 1/2) 27 


where //,, the stagnation enthalpy, is constant along a 


streamline. 


Nonisentropic Steady Two-Dimensional Flow Starting 


From Uniform Conditions 

For a flow starting from uniform conditions, such as 
a body flying through the atmosphere, as all stream 
tubes originate in a region of constant stagnation en 
thalpy JJ), 7, will be a constant everywhere in the flow 


; Basa field. In this case we can differentiate Eq. (27) in any 
In this case we require another equation, which involves ere aS ee ee 
: 4 ; arbitrary direction; so 
an assumption about the flow field. We assume here , 
that there is no heat transfer across streamlines, so that 0 = d(E + RT) + wdw + dl 
and 
My t V i = l qw?) | — WOW = él - (VV AJ? _ l w)T dn 
= +(V AF — 1/2)6 In [w?/(2H, — 2U)| + V AP? — 1f6U/(2 — 2U))| = (VAP w*)dl 
VM w?) 7 én 
But 
w?/(2H, 2U) = M?/) M? + [2(E + RT)/a* |i 
so that 
l l l 
3 : | —e Us 
97? —-1 UF (I? - 1d/f1+ [2(E + RT)/a*}}) 1+ (M@? — 1) 
V MP 1 (E+ RT) VM? —1 { (2A + RT)/a*| |. VM _ 
; ana 0 * ae 6l° F 1 én 
M? + [2(E + RT)/a’ a Ww" (i? + [2(£ + RT) /a*\! t 
ind finally, if we let E + RT = II, we find 
£6) = 0/O.\?);V 1 + (2H/a*) tan-! V [1 + (2H/a?) |? — 1) — tan-! VAP? — 161 - 
iV M 1/[.M? + (211 /a?)|{6(1/a2) — (VM? — 1/w?)} (2H /a?) (MP? + (2H /a?)\j6U — (VM w*) 1 6n 
Now suppose we define a property angle v as 
y= V1 + (2H/a?) tan-' V[1 + 2H /a “AM? — 1) — tan '' VM (28) 





790 JOURNAL OF THE AERONAUTICAL SCIENCES NOVEMBER, 1957 


so that v = v( MJ", 277/a"); then, 


6(v + #) = (2 [dv 0(2/T/a*) \yy2 — VANE — 1/[l/*? + (2H a*) |} )6(H7 a*) = 
(V M2? — 1/w?)) (2H/a*)/[M? + (2H /a?)|{6U — (VAP w-) Ldn 


Dissociating Diatomic Gas With Vibrational Heat-Capacity Lag 


Here, as before, R = Ro(1 + a), EF = [((6 + a@)/2)RoT, and U = aq + E,(0,), so that 
' RT (7 + 3a)(OT/Os) + T(Oa/ds)} 1 + (2ge/RoT)} + 2(C,i,/Ro)(00,/0s 
a> = (1+ a)Ry — = oe : 
(5 + a)(O7/Os) + Tlda Os)} (2qo RT) — ly + 2(C,3,/Ro)(00,/Os 
= BRT 
Also, | Ce 
he | (5 + a)/2]/+ (+ a){ RoT 7+ 3a rist 
a® BR(1 + a)T ~ 2801 + a) supe 
and 6U = gqda + C,60 vert 
ni 
Although we could eliminate ///a* from Eq. (29), it is more convenient te leave it and prepare a graph of » foe Strat 
lit ¢ 


> 


various values of 2/7/a*, which varies from 3 to possibly 10 depending on the circumstances. Fig. 1 gives a chart 


v versus .\/ for various values of 2/7/a*. Eq. (29) for this case is 


Ov . (2 VM? — 1 j FT 217 (qoda + Cyn 60,)) Tin p omidy 
d(v + dv) = 6 —} — : — a5 _ — Vi? - 1 3 ving 
O(2H/a?) Jy: \ a? M2 + (22H/a) U \e a? w? f w givin 


Note that for a = 0, 9, = 0, 8 = 7/i and H/a? = 
5/2; then for 7 = coistant we have the well known re- 





sult for a diatomic gas—namely, 6(v + ¢) = O along servations of flow phenomena can lead, via the theon 
. ‘ : | one = to conclusions concerning the chemical kinetics 

the directions dn/ds = *#1/V JA/? — 1, where vy = V6 - 

tan-! V (1/6)(M7? — 1) — tan-! V MW? — 1. The case REFERENCES 
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ANTONIO FERRI,* JOSEPH H. CLARKE, ** 


Favorable Interference in Lifting 
Supersonic Flow’! 





Svstems 1n 


AND LU TING*** 


Polytechnic Institute ef Brooklyn 


SUMMARY 


Certau imple integral relations are applied to the qualitative 
effects 


various spatial arrangements of wings and bodies in a 


ind ecaleulation of the favorable interference 


ssessiiit 


rising in 


supersonic stream. The arrangements include slender bodies in 
proximity to wings, wings with underslung pylons or other 
vertical surfaces, and wing-fuselage combinations. The ele 


msiderations used in these examples serve to demon 
strate to the designer how the pressure drag due to volume and 


lift can be reduced by arranging the several components com 


configuration so that they exhibit mutually advan 


prising the 


tageous shapes and orientations. Vertically antisymmetric 
nodifications of the fuselage shape in the neighborhood of the 
midwing are suggested which reduce the drag due to lift of mid 
wing-body systems The considerations are different from those 
giving rise to the supersonic area rule, according to which verti 
illvy symmetric fuselage indentations are made in order to reduce 
the drag due to thickness or volume, and might have comparable 


signincance 


SYMBOLS 


| = area 


| = reduced aspect ratio 


V MW l 
span 
( drag coefficient, based on exposed wing area 
( lift coefficient, based on exposed wing area 
es cl ord 
ti - 2u/V, pressure coefficient 
D = drag 
F = lineal source intensity 
surface source intensity 
f image source intensity 
h semiheight of fuselage 
L lift 
== Mach Number of undisturbed stream 


dynamic pressure of undisturbed stream 
dimensionless height of wedge; see Fig 
5 cross-sectional area 
u component of perturbation velocity 
velocity of undisturbed stream 
y-component of perturbation velocity 
component of perturbation velocity 
streamwise Cartesian coordinate 
lateral Cartesian coordinate 
rmal Cartesian coordinate 
ngle of attack of wing 
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Annual Meeting, IAS, New York, January 28-31, 1957 
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6,, 6; = inclination of bottom and top of body, respectively 
positive when opening body lL, 2, 

¢ = normal Cartesian coordinate 

n = lateral Cartesian cocrdinate 

r = 02/Ox; streamwise inclination 

\,, A; = inclination ef lower and of upper sides of body, respec 
tively; positive when opening body 1,2,3 

im = Mach angle of undisturbed stream 

é = streamwise Cartesian coordinate 

o = thickness ratio of slender body 

T = thickness ratio of wing 

o = velocity potential 

bu = elemental width 

refers to isolated body 
(); refers to inverted system 


( refers to minimum value 
refers to isolated system 
denotes two-dimensional value 
( refers to wing 
denotes incident disturbance 
denotes additional disturbance 


denotes additional disturbance 


INTRODUCTION 


oe LITERATURE exists concerning the 
favorable interference in arrangements of wings 
and bodies in a supersonic stream. In reference | it 1s 
shown that it is possible to establish a wave system 
between two infinite rectangular wings in_ biplane 
arrangement such that the internal wave drag is zero. 
A number of authors have subsequently considered the 
three-dimensional interfering biplane. An axially sym 
metric counterpart of the biplane, consisting of a body 
of revolution enveloped by a cowling, is proposed in 
reference 2 and analyzed by other writers. In reference 
3 it is shown that the wave drag of an array of slender 
bodies can be reduced through use of orientations 
wherein each body is in the buoyant field of another. 
Additional examples of such favorable interference have 
been given by others. 

Recent developments make it possible to study 
spatial distributions of sources, horseshoe vortices, and 
multipoles which simulate, in principle, any array or 
combination of wings and bodies to the linear approxi 
mation. Hayes’ theorem, permitting the construction 
of equivalent lineal distributions of singularities by 
moving disturbances along certain oblique planes, is the 
fundamental device employed. It has become possible, 
for example, to construct a complementary line of 
to a given body which 


multipoles corresponding 


eliminates the wave drag of the combination, and to 
consider distributions of singularities which provide 


minimum wave drag for a given total lift and volume in 
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Fic. 1 Flow field of semi-infinite source 


a prescribed region. References 4, 5, and 6 (and the 
further research cited therein) offer accounts of these 
methods. It is now clear that unconventional arrange- 
ments of suitably shaped wings and bodies make possi- 
ble significant reductions in pressure drag. In the case 
of symmetric combinations of wings and slender bodies 
at zero lift, it is possible to relate the singularities 
studied to the geometry of the system and there results 
the supersonic area rule. When the bodies are non- 
slender oi the system develops lift, the geometry simu- 
lated by the singularities is not usually known or con- 
sidered in detail. Whereas the geometry can, in 
principle, be determined in the last step, no assurance 
usually exists that it conforms to engineering exigencies 
or even is real; there is encountered the fiction of 
negative volume, for example. This factor constitutes 
a difficulty in a significant body of research. 

In the present paper, certain simple integral relations 
for the flow field of an elemental supersonic source are 
stated in the first section and then applied in subsequent 
sections in the design and analysis of a number of lifting 
systems exhibiting favorable interference. These are 
(1) a slender body under a lifting wing, (2) a wedge or 
wedges in the vertical plane of symmetry of a lifting 
wing, and (3) a pylon under a lifting wing. From the 
integral relations, certain properties of the chordwise 
distribution of pressure in wing-fuselage combinations 
are then deduced. These are used in the design and 
analysis of (4) a lifting wing-fuselage combination, the 
fuselage of which is antisymmetrically contoured about 
the plane of the midwing. The proposed advantageous 
contouring for lifting combinations is related to the 
familiar symmetric contouring for nonlifting combina- 
tions. All configurations analyzed exhibit a single 
method for the reduction of the pressure drag of lifting 
systems, the concept being developed synthetically in 
the several examples. Whereas the drag of each of the 
arrangements considered is certainly not a minimum in 
the sense discussed above, significant reductions in drag 


‘ 


relative to the chosen standards of comparison ay 
found in all cases, the physical processes are reasonab) 


clear, the geometric boundaries are explicit, and som. § 
BS inde 


new avenues of study are perhaps suggested. 


INTEGRAL RELATIONS FOR SOURCE FLOW FIELp 


Consider a semi-infinite horizontal elemental surfa 
QR, with width 6w, which is aligned in the streamyjs 
direction at a distance ¢ from the x-y plane as shoy 
in Fig. 1. Let there be distributed over this surface 
source sheet with surface intensity /(é). The trace 
the Mach cone from Q in the x-y plane is denoted | 
ACB. At an arbitrary distance downstream of 
construct in the x-y plane a lateral line segment 4 
which runs between the two branches of the hyperbo! 
It may be verified by direct calculation that the |i 
integral along AB of either the x or the s component 
the perturbation velocity, induced by the line of dis 
turbances, is given by the product of the element; 
width 6w and the two-dimensional value of the corr 
sponding component which is associated with the sour 
intensity at the point P. That is 


2B 
| 6u dy = —([nf(P)/Blbw = u7(P) bw 
JA 


2B 
Ow dy =—_— 1f(P)éw = w(P) bw, ¢ +0 Z 
1 


The point P is defined by the intersection of the ek 
mental surface QR with the forward-leaning Mac 
plane containing the line AB (see Fig. 1). The reaso: 
for the appearance here of the two-dimensional value: 
suggested by the following observation: the integr 
of the potential due to a point source at Q(é, 7, 6), 


V (x — &)? — B?[(y — n)? + (2 — 67] 


with respect to y over the limits defined by the vanisl 
ing of the radicand is the same as the integral of th 
potential with respect to 7 over similarly defined limits 
since the roles of y and 7 are interchangeable, the tw 
dimensional value appears in each case. Eqs. (1) at 
2) give some insight into how disturbances spre 


laterally in three-dimensional flow. 

Provided that the distance ¢ is nonzero, results 
and (2) are unchanged if the surface element QR 
rotated about its streamwise axis to an arbitrary pos 
tion, or if the distribution of disturbances is regarded @ 
a line distribution instead of an elemental suriat 
distribution. In the latter case, one lets 


S(E)bw — F(E) 


where F(é) is the lineal intensity of the source distribu 
tion. 

If the source intensity is independent of £, the int 
grals of induction, along the line segments defined ! 
successive Mach planes and the hyperbola in the 
plane associated with Q, are the same. A rectangul! 
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FAVORABLE INTERFE 


element with area 6f6w (See Fig. 2) may be formed from 
two such uniform distributions by subtraction. From 
Eqs. (1) and (2), there follows immediately that the 
induction due to the disturbances distributed over this 
rectangular element is such that 


6udA = —(nf B)btdw (5) 
JA 
6w dA = Ff dté6w (6) 
JA 
»\ 
budy = 0, x > Bi (7) 
JI 
*\ 
ébwdy = 0,x > Be (S) 
Ji 


In the above, A is the area defined by the two hyper- 
bolas and the lateral tangent shown in Fig. 2. The 
lateral line segment L.V is any line downstream of A 
which extends from one branch of the upstream hyper- 
bola to the other branch. The plus sign in Eq. (6) is 
taken when ¢ is negative. This subtraction process is 
used to avoid singularities connected with the induction 
between the two hyperbolas when, in the limit, the two 
merge into one.’ As mentioned above, the element 
may be rotated about its streamwise axis to an arbitrary 
position or reduced to a line supporting an elemental 
lineal source distribution; as will be shown later, the 
integral relations can be used to study diffraction around 
a corner in connection with wing-body interference 
problems. 

The relations (5)-(S) provide a simple, physical 
means to assess and calculate representative examples 


When 


the geometry is appropriate, they make possible the 


ina number of classes of interference problems. 


computation of the lift force, and sometimes the drag 


force, from a knowledge of the source distribution 
satisfying the boundary conditions without computing 
the induced flow field. A number of such applications 
is presented in the following sections of this paper. In 
the procedure used, the induction of each elemental 


source is integrated before consideration is extended to 
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Fic. 3. Slender body under wing 


all elements of the distribution. This reversal of the 
order of integration may be interpreted in terms of a 
finite-step approximation to the actual source distribu- 
tion, the limiting process being reserved for the final 
operation.” It is possible to generalize the integral 
relations to the case where the line segments, such as 
AB in Fig. 1, are swept. Various features and applica 
tions of these relations are discussed in references 7, 8, 


9, and 10. 


FAVORABLE INTERFERENCE BETWEEN SLENDER Bopy 
AND LIFTING WING 
Consider a rectangular wing having a chord c, a span 


b, a symmetric diamond profile with thickness ratio 7, 
and an angle of attack a which is sufficient to generate a 
prescribed lift. Assume that the reduced aspect ratio 
A is greater than two, so that the tip effects do not over 
lap. Suppose that a given symmetric body of revolu 
tion, with thickness ratio o, is to be positioned under 
the wing such that the atterdant rise in pressure drag 
is as small as possible. For simplicity of discussion, 
assume that the length of the body is equal to the 
chord of the wing (see Fig. 3). 

Consider the chordwise load distribution induced on 
the wing by the body, as the body is moved beneath the 
This loading in any 
and (7) [for 


wing in a streamwise direction. 
position is given immediately by Eqs. (5) 
Eq. (1)] in terms of the known source distribution 


representing the body (in this application, the rec 
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Drag rise for body under wing. 





Fic. 4 


tangular element in Fig. 2 is reduced to a line element) ; 
it is seen that each elemental source of the forebody 
makes a positive contribution to the load, whereas each 
elemental sink representing the afterbody makes a 
negative contribution. (Away from the surface of the 
body it is sufficient to use the linear form of the pressure 
coefficient.) Evidently, the value of the normal force 
is a maximum when the body is positioned so that the 
nose Mach cone strikes the shoulder of the wing profile 
and the midbody Mach cone strikes the trailing edge 
of the wing. Furthermore, the interference force for 
this position corresponds to a lift and, when 7 > a, toa 
thrust on the afterwing. Since the boundary shape of 
the wing is to be preserved in the presence of the body 
upwash, the integral of which is given by Eqs. (6) and 
(8), the load distribution must be doubled in order to 


account for the reflections. Using Eq. (4) in Eq. (5), 


noting Eq. (7), and integrating from zero to c/2, one 
obtains 
| udA = —(r B) | F(&)dé (9) 
eo & /7 0 
After use is made of the relations 
F(x) = (V/2n)S'(x) (10) 
Cp = —(2u/V) (11) 
and account is taken of reflections, there results 
| Cp dA = (2/B)Snaz (12) 


for the integral of the pressure induced by the slender 
body over the lower afterwing. The force components 


associated with Eq. (12) are 
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Liyo/C = (2/B) Swaz 


—(2/B)(r — @)Spnaz 


The orientation of the body with respect to the wing 
is not completely fixed by the above considerations 
It is seen from Fig. 3 that the body may still be trans 
lated along the Mach waves until the Mach wave fron 
the wing leading edge strikes the midpoint of the body 
and the wave from the wing shoulder strikes the tail, 
the body. The two-dimensional positive pressure fiel 
induced by the forewing thus imparts a thrust to th, 
body which is given by the relation 


p... g = —(2/B)(r a) Sma 





For simplicity, the uniform downwash from the wing 
is permitted to camber or kink the body at its mid 
This does not alter the transverse cross section 
The self-induced drag of the body in th 


point. 
of the body. 
presence of the wing is the same as the drag of the un 
cambered isolated body (the source distribution retains 
the requisite ‘‘smoothness”’ in the presence of the dis 
continuous downwash). In the case of a Sears-Haack 
body, this drag is 


D,/q = (92/2) (S*maz/C?) I( 


The lift of the wing is given by the sum of its self 
induced lift, corresponding to the angle of attack a 
and the lift imparted to the wing by the body, given by 
Eq. (13). Because the lift of the system is prescribed 
the introduction of the body permits a diminution oi 
the wing angle of attack, with a corresponding decrease 
in its drag due to lift. The drag of the system is giver 
by the sum of the isolated drag of the wing, correspond 
ing to an angle of attack a and thickness ratio 7, and the 
components given by Eqs. (14), (15), and (16). If the 
angle of attack is eliminated between these two expres 


sions, the drag coefficient in polar form results. Sub- 





traction from this expression of the drag of the isolated 
wing, for the same value of the lift coefficient, gives th 


following expression for the drag rise, relative to th 
isolated drag of the body: 


(Cp — Cpr)/Cp = 1 — (32/9277)(7/Bo?) (S/Qn*) X 


slr Bo*)(BC,/r) — (w/4A)]) i 
\ [1 — (1/2A)] 


In the above, the lift coefficient is based upon the wing 
area; A denotes the reduced aspect ratio of the wing. 





For the case 
BC;/r = 28 Is 


the variation of the relative drag rise with the paramete! 
V B/r o is given in Fig. 4 for values of the reduce 
aspect ratio of both two and infinity. In order t 
facilitate interpretation, a second scale for the absciss 


which reads directly in terms of o for the case \/ = 2 


and + = 0.07, is also provided. In the latter case, it® 
noted that the pressure drag rise is zero for a bod! 
thickness ratio of about 15 per cent. For thicknes 


ratios less than this value, a great deal of thrust appeals 
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FAVORABLE 


Even for a thickness ratio of 20 per cent, well outside 
the range of the theory, the drag rise appears to be only 
half of the isolated value. It is noted from Fig. 4 that 
the influence of aspect ratio on the relative drag rise is 
slight. The numerical calculations indicate that the 
contribution by the body to the wing lift is responsible 
for about 35 per cent of the reduction in drag, for the 
value of the lift coefficient considered. 

A body under an elliptic wing with zero thickness and 
constant loading is studied in reference 4, using Hayes’ 
theorem, in connection with the use of multipole dis- 
tributions to reduce the wave drag. Vortex drag is not 
considered. Because the loading of the wing in such 
analyses is specified, rather than the geometry, the body 
flow field deforms the wing instead of loading it. 


THREE-DIMENSIONAL BIPLANE 


A number of three-dimensional biplanes of the Buse- 
mann type are analyzed in reference 8. The con- 
figurations considered involve wings with swept leading 
edges and straight, unswept trailing edges. The high 
and intermediate ranges of aspect ratios are treated. 
Positioned over each wing is a small, flat reflecting plate 
which intercepts forewing disturbances and _ reflects 
them to the afterwing, providing a thrust contribution. 
For the thickness distributions considered, it is found 
possible to compute the interference effects, without 
determining the flow field, by means of Eqs. (5)-(8). 
In these applications, the rectangular element of Fig. 2 
is used in the horizontal position—1.e., parallel to the 
plane of the wing. Results indicate that the thickness 
drag of the wing can be greatly reduced, or that the 
thickness of the root section of the wing can be 
considerably increased without penalty in drag. The 
details may be found in reference 8. 


VERTICAL SURFACES USED TO INDUCE LIFT ON A WING 


In the case of the slender body under the wing, some 
of the reduction in drag resulted from the lifting pres- 
sures induced over the wing surface by the forebody. 
It appears to be of interest to study further the possi- 
bility of offsetting the drag rise due to adding a given 
body to a lifting system by using an arrangement 
wherein the body induces lifting pressures on the wing. 
This approach would be of particular interest in cases 
where the components of the system do not exhibit 
inclinations which offer the possibility of obtaining 
thrust contributions from remotely induced 
disturbances. More generally, it might be conjectured 
that the drag of a lifting system can be reduced if 


direct 


primarily nonlifting components, such as the fuselage, 
are designed to produce lift indirectly in this way. 

A simple system which is suitable for such a study is 
Shown in Fig. 5. It consists of a wedge, with a semi- 
included angle \;, mounted beneath a sweptback flat 
wing, with angle of attack a, in its vertical plane of 
The wedge is truncated at the tip as 
shown, in order to study the effect of its height, and 


symmetry. 
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along a Mach wave forming its trailing edge, so that 
every part of its surface will influence the wing. The 
leading edge of the wing is swept at the Mach angle so 
that every portion of the wing will be influenced by the 
wedge. The wing trailing edge is formed by the hyper- 
bolic trace of the Mach cone from the point of inter 
section of the lines defining the leading and trailing 
edges of the wedge, and by an arbitrary lateral line. 
The flow field induced by this configuration, and by 
similar ones, can be determined by superposing the flow 
fields due to (1) the wing alone at angle of attack a, and 
(2) the wedge in the presence of the wing when at zero 
angle of attack. Because no communication between 
the top and bottom of the wing exists, both subprob- 
lems are readily handled with wing theory. For the 
case ¢ cs = 1, the integral over the lower surface of the 
wing of the pressure coefficient induced in the plane of 
the wing by the wedge (second subproblem) can be 
written down immediately by use of Eqs. (5) and (7); 
the integral is equal to twice the two-dimensional pres- 
sure coefficient, corresponding to the inclination Aj, 
times the elevation area of the wedge. The factor two 
accounts for the image of the wedge in the plane of the 
The rectangular element of Fig. 2 is used verti- 
Thus the tft, as 


wing. 
cally in this application (see Fig. 5). 
well as the drag, induced on the wing by the wedge is 
readily determined and is evidently appreciable. As 
pointed out in reference 7, the self-induced lift and drag 
of the flat triangular wing can also be written at once in 
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Fic. 5. Wedge under wing. 








796 JOURNAL OF THE AERONAUTICAL SCIENCES NOVEMBER, 1957 





} Coy CoCo, 





GIVEN C, >0 





wm 











Fic. 6. Variation of drag with inclination 


the same way. When & = 1, the flow field induced by 


the wedge is evidently two-dimensional; from a con 


sideration of this field, it may be seen that a greater 
amount of lift can be obtained from a given wedge if the 
trailing edge of the wing is extended rearward to the 


hyperbolic boundary of the two-dimensional region 
(see Fig. 5). The trailing edge of the wing was designed 
as specified above in order to provide this feature. 

The computed drag Cp of the configuration may be 
compared to the drag C,, of the isolated system .e., 
the sum of the isolated drag of the same wedge and the 
isolated drag of the wing when generating the same lift. 
The drag of the configuration can also be compared to 
the drag Cp, of the wing alone when producing the same 
lift. The latter provides a standard for comparison 
when the wedge is regarded to exist primarily to pro- 
duce the lift. 

The force coefficients for interfering configurations of 
the type discussed evidently assume the forms 


Cr = Aga®” + AgyQAy + Aygar + Ay)Ai* (19) 


aa 
Cr = Aaa® + AyoA1 (20) 


where the a’s are all positive (here functions of ¢ ¢ 
and or R) and where 


darn * Qa 


in general. Evidently the two ad; drag components, 
as well as the \;" component, tend to offset the favorable 
lift contribution of the wedge to a certain extent in the 
present simple design. From Eqs. (19) and (20), the 
expression for the drag in polar form is 


Cy = (Cx?/daa) + [dex — Ara) CaalaCr + 


[(daarn — Garda) /Aaaldr.? (2 


As long as Bye > Aan 2 


so that the coefficient of A4,C;, is negative, some adya; 
tage due to the interference will always appear. (Th 
means that the life due to the wedge can outweigh ¢ 
aX, interference drag terms. If the reverse of Eq. (2 
were true, then a corresponding advantage wou! 
appear when the wedge is placed above the wing.) 7) 
design of the interfering configuration of Fig. 5 
adequate to fulfill Eq. (22). Eq. (21) has the form 
the expression for the polar drag of a warped plar 
wing without thickness. The drag coefficient C), 


the isolated system assumes the polar form 
Cos = (Czr?/dew) + @'yA1 9 
where, in the present example, 
a’y, < ayy 


The drag coefficient Cp, of the wing alone is given | 
Eq. (21) with A; 0. 

The properties of the three configurations are cot 
veniently exhibited in terms of two diagrams. Fig. | 
shows, for a given value of C,, the variation of the drag 
coefficient with the wedge inclination \, for the thre 
cases. In the case of the specific configuration dis 
cussed, Rand c c) are regarded as fixed. Evidently th 
drag of the interfering configuration is less than that oi 
the isolated wing over a range of positive inclinations 
Its drag is less than that of the isolated system over the 


entire range of positive inclinations of interest, and 


greater than that of the isolated system for all negative 


inclinations. This diagram provides the basis for 
selecting a wedge angle at some design value of (; 
Fig. 7 shows, for a given positive value of \j, the var 
ation of the drag with C, for the three cases. Also 
shown is the drag coefficient Cp, of the configuration 
formed from the one in Fig. 5 by removing the positive 
wedge from the bottom and placing it on top of the 
wing. The curve for Cp, is evidently the reflection 0! 
the C, curve about the abscissa. All of these parabolas 
are congruent, being merely displaced relative to one 
another. The interfering configuration exhibits less 
drag than its inverted counterpart for all positive 
values of C,. As C, increases from zero the con 
figuration successively overtakes and excels first thé 
isolated system, and then the isolated wing. TW 
quantities are suggested by these diagrams which are 0! 
special interest. The first is the minimum value of 
Cp/Cps corresponding to an optimum value of the ratio 
\1/C,; these two quantities are denoted by the sub 
scripts m and sm, respectively. The second is the 
minimum value of Cp/Cp, corresponding to an opt 
mum value of the ratio \; C,; these two quantities art 
denoted by subscripts m and wm, respectively. 


sketches of Figs. 6 and 7. Fig. 8 shows, for the spect 


configuration discussed, the variation of (Cp Cp 
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Fic. 7. Variation of drag with lift 
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Drag reductions for wedge under wing. 
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Fic. 10. Pylon under wing 
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and (Cy Cpw)m with R for values of ¢/c) of 1.0 and 1.5. 
Both ratios decrease with increasing R. Evidently the 
effect of sweeping the trailing edge is beneficial. It is 
noted that the value of (Cp Cp,)» falls as low as 0.68 
and (Cp Cpw)m drops to 0.88. 

It appears that the drag rise due to adding a vertical 
surface such as a pylon can be considerably reduced by 
establishing interference of the type discussed. Fur- 
thermore, one might interpret the wedge as representing 
the sides and bottom of the neighboring fuselage surface 
of a high wing-body configuration. It would appear 
then that drag reductions can be obtained by bulging 
the fuselage section beneath the wing of a high wing- 
body configuration or by converting a low wing con- 
figuration into a high wing configuration. 

An obvious modification of the interfering configu- 
ration shown in Fig. 5 is to add a negative wedge to the 
top of the wing in order to double the lifting force 
induced on the wing by the vertical surfaces. The con- 
figuration then becomes vertically antisymmetric, and 
the coefficients d,,, d,, and d,, in Eq. (19) and (20) 
become twice their former values. Its characteristics 
are qualitatively the same as previously discussed. 
The only standard of reference which appears to be of 
interest is the isolated wing producing the same lift, 
since the total change of volume due to the addition of 
the wedges is zero. In Fig. 9 the variation of (Cp 
Cpw)m With R for values of c/cy of 1.0 and 1.5 is plotted. 
In the latter case, (Cp /Cpx)m falls to 0.73 when R = 1.0. 
The wedges in this configuration might be interpreted 
as the antisymmetric contouring of the sides, top, and 
bottom of the fuselage near the wing in a midwing-body 
combination. These results suggest the possibility of 
important reductions in the drag of such lifting systems 
and will receive further attention in a later section. 

Configurations of the type discussed herein were 
first proposed in references 11 and 12 and are also 
studied in reference 13. A configuration composed of a 
circular half cone mounted on the lower side of an 
arrowhead wing with sonic edges is considered in 
reference 14. The cone induces lifting pressures on the 
lower wing and supports some lift over its own surface 
as well. The drag of the configuration Cp, is compared 
to that of its inverted counterpart Cp,, to the drag Cp, 
of a configuration employing a symmetrically disposed 
cone of equal volume, and to the drag Cp, of the wing 
The polar diagrams of these four configurations 
For given values of the Mach 
form parameter, the 


alone. 
are as shown in Fig. 7. 
number, cone angle, and wing 
values of C, for which the ratios Cp/ Cpr, Cp’ Cps, and 
Cp Cp, are a minimum may be computed. The best 
values found in reference 14 for these ratios are 0.64, 
0.89, and 0.95, respectively. It is again found that the 
effect of sweeping the wing trailing edge is favorable. 


FAVORABLE INTERFERENCE BETWEEN PYLON AND WING 
By use of the ideas developed in the previous Section, 


it is possible to design a closed pylon positioned below a 
flat plate wing at incidence in such a way that consider- 
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able favorable interference is established at the desig 
Mach Number. Such a pylon might be used for 4 
support of an engine or other installation. 
be interpreted as any other vertical underslung surfa 
In Fig. 10 is shown a rectangular wing in combinat 
with a simple pylon, each side of which is formed | 
A C0 


stant positive streamwise inclination A, appears 


five planes. The four swept edges are sonic. 
either side over the region OAD (in this paper a bo 
inclination is counted positive when it tends to open t 
body, and is counted negative when it tends to close t 
body). This region induces lifting pressures on { 
wing which are readily computed by the integ 
relations. This lift is produced in a relatively efficie 
manner because the region does not support win; 
induced interference drag. The pylon is gradual 
closed over the region CDBE by means of two plat 
on either side, the inclinations of which are —) 
The region CDBE induces no negative lift on the win; 
and receives a thrust contribution from the ty 
dimensional pressure field induced by the lower win 
surface. The streamwise inclination of regions O( 
and ABD is zero. 

The configuration may be analyzed by the meth 
discussed in the last section. The calculations are ver 
simple; they are given in reference 12. The expres 
sions for the lift and drag assume the forms given j 
Eqs. (19), (20), and (21). The configuration consider 
is of special interest for the following reason: — the drag 
of the body in the presence of the uninclined wing is 
virtually the same as that of the isolated body, asi 
indicated by the asymptote in Fig. 11; also, the tw 
interference drags cancel. When the isolated py! 
and wing are brought together, the increment in lit 
without sensible concomitant e 
traneous effects. When the aspect ratio of the wing: 
regarded to be large, the performance of the configur 


therefore appears 


tion may be evaluated according to the relative dra 
rise criterion used previously in the case of the interfer 
ing wing and slender body. (This parameter plays t! 
role of the parameter Cp/Cp, but is more suited t 
interfering systems with disparate dimensions at 
loads; the latter parameter must tend to one as th 
aspect ratio tends to infinity.) 

Fig. 11 presents the variation of (Cp Cox) /C 
with the characteristic parameter \;/ BC, for the cas 
where the reduced aspect ratio is infinite. As _ pre\ 
ously remarked, the influence of aspect ratio on this 
variation is small. Shown on the curve is the value 
(Ai BCz)wm, Corresponding to which the ratio Cp C 
is a minimum. It appears from the curve that, whe 
the thickness ratio of the pylon is on the order of typi 
wing thickness ratios, then the relative drag rise can! 
significantly less than one. When the thickness of th 
wing is taken into consideration, it may be expecte 
that pylons which exhibit even more favorable interfer 
ence can be designed. This point is illustrated by th 
wing-slender body configuration previously discussed 

It is of interest to present here results for the cas 
where the pylon is mounted beneath a triangular wit! 
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Fic. 12. Diffraction of ¢p around right corner 


with sonic edges. It is found that 


(C. (Cod = 690 


0.713 (25) 


(Cp Cy lm = 


When a negative pylon is placed above this wing 1n 
addition, we pass to the case of antisymmetric contour- 


ing in the vertical plane. In this case, 


O.S49 (26) 


(Cp Cou \m = 


These results may be compared to those given for the 


configuration illustrated in Fig. 5 for R = ! 


and ¢/cy = 
|. The comparison gives an indication of the effect of 
closing the wedge(s) in the presence of favorable inter- 
ference. 

ANTISYMMETRIC CONTOURING OF LIFTING WING- 
FUSELAGE COMBINATIONS 


It has been established that appreciable reductions 
in the pressure drag of lifting cruciform systems are 
possible through use of appropriate antisymmetric con- 
touring in the vertical plane. The contouring serves 
either to induce lifting pressures on the wing or to 
receive thrust contributions from the wing-induced 
flow field. 
posed of wings midmounted on semi-infinite quast- 


We now turn to lifting configurations com- 


prismatic bodies of essentially rectangular cross sec- 
tion—i.e., the body surfaces differ so slightly from the 
surface of the basic rectangular prism that the boundary 
Only 


the case where the prism is at zero angle of attack will be 


conditions may be imposed on the latter surface. 
considered. By extending the method discussed, the 
bodies will be antisymmetrically contoured about the 
plane of the wing in order to reduce the pressure drag at 
given lift. 


Determination of Flow Field 


The flow field of this system may be determined by 
superposition, once the solution to the problem of the 
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diffraction of disturbances around a semi-infinite right 
convex corner is available. The solution of this prob- 
lem is obtained in reference 15 by an extension of the 
familiar diaphragm method of reference 16. The pro 
cedure may be briefly indicated with the aid of Fig. 12 
The disturbance field incident to the right corner 1s 
given by the velocity potential ¢), which crosses the 
corner atx = 0. The Figure suggests the wave pattern 


in a plane x = constant > 0. The potential ¢, is 


defined for z > O and satisfies the boundary conditions 


(Ody OZ) (x, v < O, O) VA(x, y) 27) 
(Ogo/Oz)(x, vy > 0, O) = O IS) 


The resulting discontinuity in pressure in the x-y plane 
for vy > 0 is relieved by the back-to-back source distribu 
tions —/(x, y) and f(x, y), and by the distribution 


f(x, vy); these are indicated in the Figure. The satisfac 


tion of the boundary condition 


(0g /Oy)(x, 0, s < 0) = O 29) 


is assured by the distribution f(x, y), which is the image 
of f(x, y) and therefore satisfies the relation 


F(x, y < 0) = f(x, —y) (30) 


Appearing inside the Mach cone springing from the 
corner are the additional potentials ¢), defined for z > 0 
Q and 
The additional poten- 


and induced by —/(x, y), and @», defined for s - 
induced by f(x, vy) and f(x, y). 
tials are determined from the integral condition that the 
aforesaid pressure discontinuity vanish. 

It is convenient here to regard the solution to the 
wing-body problem to be given by the superposition of 
the solutions to the following three subproblems: (1) 
the wing panels supporting deflections in combination 
with the body when free of deflections, (2) the wing 
when free of deflections in combination with the body 
when supporting only side deflections, and (3) the wing 
when free of deflections in combination with the body 
when supporting deflections only on the top and bot- 
tom. Each subproblem is subject to further decom- 
position. The solution of each subproblem is begun 
with the solution to the planar problem associated with 


the deflections considered, these incident disturbances 
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giving rise to systems of horizontal and vertical dia- 
phragms at the four corners of the body. 
problems considered, no interference exists between the 
flow field above the wing plane and the flow field below 
the wing plane, the wing may be regarded as a plane of 
symmetry in each half space and the image artifice 


employed. 


Integral Relations 


The integral pressure relations for an elemental 
source [Eqs. (5) and (7)] yield considerable information 
about the chordwise distribution of pressure in the sys- 
tem without determining the interference flow field. 
for example, the diffraction phenomenon 
shown in Fig. 13. Assume that all faces of the body 
are free of inclination and that the flow field of the 


Consider, 


planar wing panel has been determined for the case 
where the side of the body shown extends vertically to 
infinity, behaving like the vertical plane of symmetry of 
the wing (first subproblem). This field corresponds to 
the incident disturbance at the corner D of the body 
and gives rise to the vertical diaphragm D/A inside the 
Mach cone from the point D. An oblique elemental 
source on either side of the diaphragm is shown at point 
P. It gives rise to the two illustrated hyperbolic traces 
on the top of the body and on the plane of the wing. 
Let 77 be a lateral line segment on the top of the body 
which extends from the Mach wave from PD to the 
corner. Let WG be a lateral line segment extending 
from the wing-body junction to the hyperbolic trace 
from D in the plane of the wing. Let the line VG be 
displaced downstream from the line 77/ a distance Bh. 
It follows immediately from the integral relations that 


“ °G ~J 
{ udy = — | 2u,dy = —(r B) | f dz (31) 
J1 Ju I 


The subscripts correspond to the additional potentials 


shown in Fig. 12. If the plane of the wing may be 


regarded as a plane of symmetry, the additional induc- 


tion thereon is 2u,. Other incident disturbances 


which generate vertical diaphragms arise from subse- 





When, in the 
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quent reflections and refractions; the associated adj 
tional disturbances bear the identical relationship, 
An incident disturbance propagating across the t 
of the fuselage and giving rise to a horizontal diaphrag, 
indicated by D/A, is shown in Fig. 14. This 
field might be due to the far wing panel or to ineliy 
tions appearing on the top of the fuselage. Assup 
that the plane of the wing is free of other disturbang 
serving only as a plane of symmetry. Let 77A be 
horizontal, lateral, elemental strip extending from o 
Mach wave from D in this plane to the other. Let Jj 
be a similar strip, displaced downstream from the fir 


j 
' 


by a distance Bh, which extends from the wing-bo 
junction to the hyperbolic trace from D in the plane 
the wing. Then the integral pressure relations, applic} 
to the three diffraction source distributions, indicat 
that the load induced over strip /A by the incide 
disturbance reappears over strips 77 and IVG when th 
additional disturbances are introduced to relieve this 
This fact may be demonstrated with the ai 
By definition, the additional dis 


loading. 
of Figs. 12 and 14. 
turbances must be such as to fulfill the condition 


nA +A 
ax | ugly + dx { u(—f)dy = 


I I 
7A 


oA 
dx | u(f)dy + ax | 
1 Ji 


e 


u(f)dy (32 


where the notation u(—/f) denotes the value of the 
component of the perturbation velocity induced by th 
distribution —/(x, vy) over the strip, etc. Use is noy 
made of Eqs. (5) and (7) for the case where ¢ in Fig.? 
is zero. The induction in the plane of the diaphrag: 
due to the distributions —/ and / is such that 


*/ +A 
av | u(—f)dy + dx } u(—f)dy = 


1 J / 


74 
—(mr B)dx | (—f)dy (3: 
I 


e/ | 
dx { u(f)dy + dx | u(f)dy = 


7 I 
+A 
—(m B)dx | I dy (34 
Bi 
Also, by symmetry, 
a4 2] 
dx / u(f)dy = dx | u(f)dy 3 
J] BS 


Substitution of Eqs. (33)—(35) into Eq. (32) yields 


24 e/ 
dx } udy = dx | u(—f)dy — 
JI J1 
el 
(7 B)dx | fdy — (r Bde | 
J7 


Jl 


vA 


/ dy 0 


Use is now made of Eqs. (5) and (7) for the case ¢ = 
With application of vertical symmetry and a changet 
the subscript notation of Fig. 12, there results 


A 2] 2G 
dx [ Uy dy = dv | mdy + dx | 2u.dy (di 
Bl 1 on 
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Zecause the plane of the wing is a plane of symmetry, 
the reflection of the additional disturbances therefrom 
has the effect of doubling the pressure induced over the 
strip IG by the diffraction distributions f and f. The 
pressure 1S then proportional to 2m. Therefore, Eq. 
37) is the desired result. 

When the geometry is appropriate, the chordwise 
distribution of pressure in each of the three subprob- 
lems enumerated above follows immediately from Eqs. 
21) and (37), Consider the 


laterally symmetric configuration shown in Fig. 15. 


applied. 


repeatedly 


The basic rectangular prismatic body is at zero angle of 


attack. The wing leading edge LE is everywhere 
forward of the hyperbola /7Y from the point D. The 
shape of the trailing edge is at first arbitrary. Con 


sider a forward-leaning Mach plane which intersects the 
wing-body juncture at right angies and which is every- 
where upstream of the trailing edge. 

|) Let the wing panel planes support an arbitrary 
distribution of thickness and angle of attack and let the 


body surfaces be undeflected. Then 


aS e/ I 
u dy + | udy = { u, dy 
JM J ( ed ( 


the integral of induced pressure along the 


(3Sa) 


That 1s, 
horizontal lines of intersection of the Mach plane with 
the combination is equal to the integral of induced 
pressure (~ u,) along the lines of intersection with the 
wing panels in the case where the body sidewalls are 
The latter integral is 
the 


leading edges of the wing panels are supersonic, this line 


extended vertically to infinity. 
determined from the planar wing problem; if 


integral can be written at once with use of Eqs. (5) and 


(7), applied in a plane. Eq. (88a) may then be written 


2S »/ 
udy + udy = 
VW ( 


« 


| 
—(mr B) { f dy (88b) 
ed ( 


where / is the source intensity appropriate to the planar 
wing problem. If the trailing edges of the wing are 
straight and unswept, it follows from the above (con- 
sidering both upper and lower half spaces) that the lift 
supported by the system, forward of the two Mach 
planes which contain the trailing edge, is the same as 
that supported by the wing panels in the case where the 
body sidewalls are extended vertically to infinity. 
When the leading edges are supersonic, this is the two- 
value. to carry 
forward some of the ideas presented in reference 7. 

(2) Let the sides of the body support deflections. 


dimensional These results appear 


Let the wing panels and the top and bottom of the body 
be free of deflection. If the additional disturbances do 
not extend forward to the wing leading edges, Eq. (38a) 
forward to 


again applies. If no disturbances extend 


the wing leading edges, then 


2] »] 
uay } Uu dy = | 
e ( a ( 


u,dy = 


I, 


—(7/B) | dz (39) 























Flow field of wing-body combination 


where uw, and / are, respectively, the value of u and the 
source intensity which are associated with the side 
deflections for the problem where the sidewalls extend 
The 


statement with respect to the lift is also similar to case 


vertically without further deflections to infinity. 


(1) above. 

(3) Let the top of the body support deflections which 
induce additional disturbances not extending to the 
leading edges. Let the wing panels and the body sides 


be free of deflection. Then 


as al 7\ 
| udy + | udy = } u,dy = 
JM a ( JM 


~(x B) | fdy (40) 
JM ; 


where “, and / are, respectively, the value of u and the 
source intensity associated with the top deflections for 
the planar problem where the top of the body extends 
laterally without further deflections to infinity. 

When no disturbances extend forward of the wing 
leading edges, these results for the chordwise distribu- 
tion of pressure can be summarized for the complete 
wing-body system by the statement that the integral of 
the pressure along the lines of intersection of the three 
horizontal surfaces with any Mach plane forward of the 
trailing edges assumes the two-dimensional value based 
on the integral of the source intensity along the lines of 
intersection of all of the surfaces with this Mach plane. 
The source intensity is to be related to the local surface 


deflections by the planar wing relation. 
Role of Antisymmetric Fuselage Contours 


The results for the chordwise distribution of pressure 
supply information on the lifting action of deflected 
portions of the wing-body surfaces and provide some 
basis for the antisymmetric contouring of the fuselage 
in order to reduce the pressure drag of the combination 
for a given lift. The information is, however, incom- 
p’ste because it does not extend sufficiently far down- 
In order to determine the drag, the surface 


stream. 
pressure distribution must still be computed, in general, 
but the labor involved when the method of reference 15 
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Fic. 16. Wing in cotnbination with antisymmetrically contoured 


fuselage 


is applied to this end is substantially reduced by use of 
these results. 

Consider the lifting problem where the wing has zero 
thickness and, for simplicity, supersonic leading edges. 
Then the sides of the body may be antisymmetrically 
deflected about the plane of the wing in order to induce 
lift indirectly on the horizontal surfaces or to receive 
thrust therefrom. The lifting action is evidently the 
same as for the cruciform systems previously consid- 
ered. Similarly, the top and bottom of the body near 
the wing would be deflected downward in the direction 
of the wing chord, so that the body is kinked near the 
leading edge and again near the trailing edge. In the 
sign convention employed, the lower deflection is posi- 
tive and the upper negative. These deflected surfaces 
then induce lifting pressures over the wing in accord- 
ance with the mechanism just discussed in terms of 
oblique Mach planes. But these surfaces are dis- 
tinguished from the sides by their ability to support 
direct lift. In fact, forward of Mach planes defining 
straight trailing edges, they behave like high aspect 
ratio wing surfaces. It is then evident that the drag 
of the system considered, at given lift, may be reduced 
by providing antisymmetric deflections of the top and 
bottom of the body near the wing which are equal in 
magnitude to the wing angle of attack. This simply 
increases the effective wing area and therefore amounts 
to passing the wing through the body in some way. A 
second method for the further reduction of the drag of 
lifting systems, which is possibly coordinate with the 
first method utilizing indirect lift, now appears. This 
involves providing a deflection distribution to the top 
and bottom of the body and to the wing which achieves 
between the three horizontal planes the same favorable 
interference effects which are derived from the camber- 
ing and twisting of a planar wing with unchanged lift. 
Because each deflected element supports only part of 
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the lift it induces, the second method also exhibits th, 





indirect lift feature. 
A fuselage with an essentially rectangular cross se 





me 


tion is used in the present study in order to divide th 
body into regions which support direct lift and regions 
which do not support direct lift. Fortunately, th 


rectangular cross section is also analytically simpler 


7 SP pT RM 


A fuselage which is deflected on the top and _ bottop 
near the wing perhaps suggests an uncontoured fuselag: | 
The upwash field of th 


which is at an angle of attack. 
fuselage at angle of attack is known to impart a 
apparent angle of attack to the wing; the resulting lif 


RT aEgS 


can be estimated by passing the wing through the bog 
in some way and considering the so-called blankete 
area. Whereas the wing is effectively passed throug! 
the body in both cases, the two mechanisms are con 
siderably different. 


Example 


Fig. 16 illustrates a particular wing-body combina 
tion which has been subjected to analysis. Each wing 
panel has zero thickness, an angle of attack a, a chord, 
a sonic leading edge, and a straight, unswept trailing 
edge. The body has an essentially square cross sectior 
with sides c B. 
c/2 from the trailing edge of the panels so that the 


It has been terminated at a distance 


domain of influence of the trailing edge does not extend 
to the top and bottom of the body. Base effects are not 
considered. Body contouring is provided which is 
considered to strike a balance between effects desired 
and analytic complexities. Beginning directly below 


the apexes of the triangular panels, six plane rec 





tangular surface elements, with deflections 6; to 4, 
The top 
of the body (not shown) is antisymmetric with deflec- | 


' 
appear on the bottom of the body, as shown. 


Below the wing on either side of the | 
i 
I 


tions —6, to —ég. 
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Wing in combination with asymmetrically contoured 
fuselage 


Fic. 17. 
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two rectangular elements and three triangular 
of 
The inclination of the sides 


body art 
elements the dimensions and orientation these 
appearing 1 the Figure. 
of the body in the domain of influence of the wing wake 
is zero. The upper sides of the body are contoured 
antisvmumetrically. The plan-form area of the wing 
panels is c? B. The plan-form area of the deflected 
horizontal portions of the fuselage is 3c?/2B. 

The flow field is given by the superposition of the 
flow fields of the three subproblems. Each subproblem 
is divided into equivalent symmetric problems, these 
then being further decomposed. Each problem is 
solved by the method of reference 15 with the aid of the 
integral relations. lift the 


minimum drag of the combination, corresponding to 


For a given coefficient, 


optimum deflections of the planar contour elements, 


mav be computed. There results for the minimum drag 


Cp/B C,? = 0.1083 (41) 
corresponding to the optimum deflections, 


B Cr = 
6. B Cr = 


a/B C, = 0.0573 & 0.0579 


1253 0.0596 


Ay B Cr ™ Q 


\/BC, = 0.1179 63/B C, = 0.0724 


(42) 


3/BC, = 0.1034 6,/B C, = 0.0895 


0.0159 6, B C, = 0.1068 


06 B CL = U0 


\4/BC, = 
—(.0047 1248 
The deflections above the wing are the negative of these 


values. For the case 


the drag is found to be 


Cp B C,? 0.1291 (44) 


for which a/ BC, = 0.1276 (45) 


When the body is entirely free of deflection, the drag 


of the combination is 


C,»/B C,* = 0.1830 (46) 


for which a/B C, = 0.2139 (47) 


The drag of the optimum design is 16.1 per cent less 
than the drag of the configuration in which all hori- 
zontal elements are equally deflected and all vertical 
elements are undeflected. Further, the drag of the 
optimum design is 40.8 per cent less than the drag of the 
If the 
horizontal surfaces assume optimum inclinations, the 


wing in combination with the undistorted body. 


vertical inclinations being suppressed, there results a 
value for the drag which is only slightly less than the 
value given by Eq. (44). The inclinations are, how- 
ever, quite different from the value given by Eq. (45). 
It is believed that Eq. (44), and not Eq. (46), provides a 
reasonable standard to which the drag of the optimum 
design can be compared. The modification of the body 
according to Eq. (43) is an obvious improvement of the 
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configuration considered which merely increases the 
The 16.1 
improvement in the drag, through use of antisym 


effective wing area. per cent additional 


metric distortions, is evidently significant. 


Relationship With Symmetric Contouring 


It is of interest to indicate for the midwing configu 
ration the relationship between antisymmetric contour 
ing to reduce the drag due to lift and symmetric con 
touring to reduce the drag due to volume or thickness. 
The now familiar body distortions of the latter type are 
given by the supersonic area rule (see reference 17) in 
the case of wings in combination with slender fuselages, 
and by more elaborate considerations when the fuse 
lages are not slender (see, for example, references 1S and 
19). 


considered. 


Fig. 17 shows the same basic configuration just 
However, the wing has been provided 
with a thickness distribution corresponding to sym- 
metric diamond sections, the shoulder falling at the 
2/5 chord point. The thickness ratio of the section ts r. 
The deflection of each planar contour element above 
the wing is now distinguished from its image below the 
wing by a bar (see Fig. 17). The flow field and drag 
for the symmetric, zero-lift problem where a 0 may 
be calculated within the framework already established. 
For a given value of the thickness ratio 7, the minimum 
drag for optimum deflections of the planar contour 


elements may be computed. There results for the 
minimum drag 
CpB 7? = 3.891 (48) 
corresponding to the optimum deflections 
Ai/r = —1.2028 67 = 0.1720 
As tT = —0.098S 6./Tr = 0.0828 
A3/T = 0.3391 62/7 0.0323 
(49) 
\y/ Tt = 0.1452 6, 7 = —0.0199 
As 7 = 0.2888 6, 7 = —0.0253 
d6 7 = —0.0052 
and Ai, =A, (¢ = 1-5) 
. (90) 
6; = 6; (@ = 1-6) 
by symmetry. When all body deflections are sup 
pressed, the drag is given by 
CpB/r? = 5.583 (51) 


The results for the deflections given by Eq. (49) as well 
as Hayes’ theorem, show the importance of the four 
leading planar elements for the plan form considered 
and suggest that a subdivision into smaller elements is 
inorder. In spite of the fact that the coarse contouring 
is not especially appropriate upstream for the thickness 
problem, Eq. (48) indicates a 30.3 per cent reduction 
with respect to Eq. (51). 

Consider, finally, the combined problem where the 
asymmetrically contoured configuration in Fig. 17 is to 
generate a prescribed lift and the wing has a prescribed 
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thickness. When the plane of the wing is the hori- 
zontal plane of symmetry of the surfaces which support 
the boundary conditions, as is the case here, an im- 
portant simplification appears. It is clear that the 
optimum deflections for minimum pressure drag in the 
combined problem are linear functions of the prescribed 
C,, and 7, and that a is a linear function of C; alone; 
these deflections are then given by the sum of the 
corresponding optima already calculated for the sepa- 
rate lift and thickness problems. The optimum flow 
field is, therefore, the sum of the two optima. 
already been noted, the optimum body deflections for 


As has 


the lifting problem, corresponding to an arbitrary dis- 
tribution of wing deflections, are antisymmetric (odd 
potential), and the optimum deflections for the thick- 
ness problem are symmetric (even potential). When 
an odd and even potential are superposed, the pressure 
corresponding to the first produces no drag on the 
boundaries corresponding to the second, and vice versa; 
the drag corresponding to the combined potential is 
then the sum of the drags corresponding to each taken 
separately, as in the planar wing problem. In _ par- 
ticular, we find that the minimum drag for the com- 
bined problem is the sum of the two minima given by 
Eqs. (41) and (48). More generally, if either Eq. (48) 
or (51) is added to Eq. (41), (44), or (46), there results 
the drag of the system corresponding to the two super- 
posed boundaries. We then see at a glance the effect 
on the drag, for given lift and thickness, of the several 


contour combinations. 


CONCLUSIONS 


The integral relations for the flow field of an elemental 
supersonic source provide a useful means for the assess- 
ment and calculation of the interference effects arising 
in various spatial arrangements of wings and bodies, 
and in wing-body combinations. The relations imple- 
mented the analysis of a number of special problems. 
It appears that slender bodies and pylons can be 
located under lifting wings with small, or even negative 
rise in pressure drag. Antisymmetric modifications of 
the fuselage near the midwing can provide significant 
drag reductions in lifting wing-body combinations. In 
the example given, the pressure drag for given lift was 
reduced 16.1 per cent below the drag of the combination 
with the wing, in effect, passed through the body, and 
40.8 per cent below the drag of the combination with an 
undistorted body. In these antisymmetric modifica- 
tions, vertical surfaces of the body are distorted either 
to induce lift on the wing or to receive thrust from the 
flow field induced by the wing. Since horizontal sur- 
faces of the body near the wing can be deflected to 
support self-induced lift, they can be provided with a 
distribution of angle of attack related to the wing dis- 
tribution of angle of attack in such a way as to provide 
mutual favorable interference analogous to that ob- 
tained by warping planar lifting systems. The anti- 
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to the 
drag due to lift that symmetric contouring bears to th, 


symmetric contouring bears the same relation 


zero-lift drag, and appears to have comparable signif 
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Rotating Stall in Axial Flow Compressors 


JEAN FABRI* ann RAYMOND SIESTRUNCK** 
Office National dad ktudes el de Recherches eronautiques 


SUMMARY 


[he classical theories of linearized aerodynamics yields a com- 
lete description of the flow field due to the rotating stall in a 
two-dimensional cascade of airfoils, when the extension of the 
stalled region is small. The corresponding physical scheme can 
be used for larger stalled regions, without the limitations of the 
inearizing assumptions, if no through flow is assumed through 

stall zones. Experimental researches on a single rotor con 
firm the results of the theoretical calculations, where the mean 


radius section of the rotor is described by means of a two-dimen 





SYMBOLS 


Db = § 


tatic pressure 


P total pressure in the absolute motion 
motion relative to the rotor 
S motion relative to the stall 
tangential velocity of the rotor 
J = tangential velocity of the stall 
if = xial velocity of the flow 
local upstream S-perturbation velocity for linear 
ized theory and actual upstream S-velocity for 
nonlinearized theory 
local harmonic parts of the downstream S-pertur 
bation velocity of linearized theory 
= asymptotic downstream rotational S-perturba 
tion velocity for linearized theory 
direction normal to the cascade (itself defined by 
= () 
y = direction parallel to the cascade 
, = R-angle of the uniform motion 
re = critical R-angle 
= R-exit angle from the cascade 
{ = stagger angle of blades 
6 = downstream S-streamline angle at infinity 
0 = angular distance between two hot-wire anemom 
eters 
= stalled fraction of a period 
p = density 
¢ = flow coefficient 
y = total pressure coefficient 
¥ = static pressure coefficient 
the sign ~ designates all quantities corresponding to the ideal 
point of retating stall onset 
(1) INTRODUCTION 
1.1) Occurrence of Rotating Stall 
WwW N THE FLOW RATE in an axial-flow turbo- 
machine decreases, due, for instance, to the 


throttling of the engine, the angles of attack of the 
blade rows increase and eventually the machine stalls. 
Very often this stall appears as a periodic disturbance 
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of the general flow and is characterized by the uniform 
propagation, from pressure faces toward suction faces 
of the blades, of symmetrically disposed stalled zones 
made up of one or several groups of stalled blades. 

The flow and pressure fields observed during such a 
rotating stall undergo periodic variations, the steep 
fronts and high intensities of which show the non 
linear character of the process. However, as long as 
the extent of the stalled zones is not too large—i.e 
at an operation close to the point where rotating stali 
appears —linearized theories may apply and yield a 
mathematical scheme for the physical aspects of the 
phenomenon. 

Such a linearized theory was given by W. R. Sears, 
who thoroughly examined all the basic assumptions 
and the results they can yield. The very important 
part played by the assumed behavior of the stalled air 
foils was clearly emphasized, as well as the eigenvalue 
character of the speed of stall propagation, which ex 
perimentally seems to be a perfectly determined param 
eter of the phenomenon. This character appears in 
the theory only if the pressure versus the angle of at 
Neither the in 
troduction of a time lag between the attack 
and the lift of the airfoil, in order to take into account 
the unsteadiness of the phenomenon, nor, as suggested 
by H. W. Emmons et al.” and later by A. H. Stenning, 


the airfoil theory by a _ blade 


tack performance curve is not linear. 


angle of 


the replacement of 
channel theory, would avoid these difficulties. 

On the other hand, F. E. Marble* succeeded in de 
scribing by means of an aerodynamic theory the main 
properties of the rotating stall, by assuming that the 
pressure versus angle of attack performance curve is 
discontinuous for the critical angle of attack. 

A result which seems new to us and which consti 
tutes the main point of interest of the linearized ap 
proach® given briefly hereafter, is that all the main 
kinematic properties of rotating stall may be obtained, 
qualitatively as well as quantitatively, without using 
any definite performance curve for the operation of the 
stalled airfoil cascade, as long as it can be assumed that 
both the stalled and unstalled performance curves can 
be linearized and that there is an actual discontinuity 
between these two operating conditions. 

Then the conclusions of our linearized theory are the 
following: when rotating stall appears in an axial 
flow compressor, the angle of attack of the unstalled 
blades remains constant and equal to its maximum 
value, compatible with a nonseparated flow, whereas 
practically no flow goes through the stalled zones. 
to a 
where 


one- 
the 


further generalized 
the flow 


These results are 


dimensional representation of 
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linearizing assumptions are no longer necessary—which 
gives the corresponding results for stalled operations 
with large stall cells. 

Both representations are described in the case of 
two-dimensional cascade flow only. 


1.2) Mathematical Representation of the Flow 


Let us consider a two-dimensional cascade of in- 
finitely close airfoils placed along the y-axis of the 
xy-plane. This cascade moves along the y-axis with 
the velocity —U(U > 0). A uniform incompressible 
flow (density p, velocity IV” > O directed along the 
x-axis) attacks the profiles. In the R-motion, relative 
to the blades, the ratio U/W defines the angle of attack 
a (tan a = U/JWW) and as long as a does not exceed its 
critical value @, the whole motion is uniform and 
permanent. 

However, when a > 4, equidistant stalled regions 
propagate along the cascade with the absolute ve- 
locity —V(V > 0). (See Fig. 1.) Then both the ab- 
solute motion and the R-motion are unstationary. On 
the contrary, the S-motion, the motion relative to the 
stall, is stationary if an infinitely close blade spacing is 
assumed. 

For the sake of simplicity, the distance separating 
two stalled zones is measured by 27, each stalled zone 
extending over a fraction v of the whole period. It is 
further assumed that the R-exit angle ax may be con- 
sidered in first approximation as constant and inde- 
pendent of the angle of attack.*® 


(2) LINEARIZED FLOW FIELD CORRESPONDING TO THE 
ONSET OF ROTATING STALL 


2.1) Description of the Flow Field 


The upstream S-motion is obtained by superimposing 
on the uniform (IV, V’) flow an irrotational perturbation 
field (w, v) periodic with respect to y and vanishing at 
the upstream infinity. 

A similar description of the downstream field can 
be obtained only in the case when the extent 27 of 
each stalled zone is small in comparison with the 
period of the phenomenon (vy < 1). Then the rota- 
tional downstream S-flow is made up of the main flow 
(W, W tan 6) which assures the mass flow rate and gives 
asymptotic streamlines with a slope tan 6, of a har- 
monic perturbation field (w’, v’) which vanishes at 
the downstream infinity, and of a rotational perturba- 
tion flow the local value of which can be derived from 
the asymptotic perturbation velocity (@, @), where 
d = @ tan 8, by calculating (@, 3) for the local value of 
y — x tan 6, since in the case of the linearized theory the 
downstream .S-streamlines coincide nearly with the 
straight lines 


y — x tan @ = const 


2.2) Correspondence Between the Perturbation 
Velocities on the Two Faces of the Cascade 


As all perturbations are periodic in y it follows that 
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W(tan ax — tan@é) = U-— V (1 


With the subscript 0 designating the values of the per- 
turbations at x = O, the following relations are ob- 


tained: 
9 nm”? 4. ££ a1, a wo «€ * 9 
Wo = Wo TT Wo, ra) + Yo = WW, tan a é 


which express the conservation of the mass flow rate 
across the cascade and the constancy of the R-exit 
angle. 
3y eliminating the rotational perturbation (@p, i = 
@) tan @) there follows the relation of correspondence 
w(tan ax — tan 0) = vy’ — wy’ tan 6 (3) 
between the two harmonic fields (w, v) and (w’, 0’ 
along their common boundary (x = 0), from which 
there results the conjugate relation (see Appendix A 


vy (tan ax — tan 0) = wy’ + w’ tan 0 f 


2.3) Energy Equations 


Now, in the unstalled regions, the energy equations 
can be written in their usual form, providing that only 
steady motions be considered and also that the same 
streamline be followed. 

In the upstream field, the S-motion is steady and, if 
P designates the absolute upstream stagnation pres- 
sure and pp the pressure at a point of the upstream 
front of the unstalled region, the Bernoulli equation 
follows in the form 


(P/p) + (V?/2) = (po/p) + [(W + wo)?/2] + 
[((V + w)?/2] © 
As the flow passes through the rotor, a steady F- 


motion can be assumed; thus, p’ being the pressure at 
the corresponding point of the downstream face of the 
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cascade, one ¢ ybtains 
p) Tt (W+ Wo)" 2] a (V+ a)- 2] = 


po’ p) + [(W + Wy)" 2) (1 + tan- ax) (3 ) 


Do 


Finally, in the downstream S-field, the asymptotic 
pressure P is related to the conditions at the cascade 
exit by the equation 
p/p) + [((W + wo)?/2] + 

[((W tan 6 + wy tan ax)*/2]} = 
(p/p) + ((W + ®)? 2] (1 + tan? 6) (5”) 

By linearizing Eqs. (5), (5’), and (5”)—that is, by 

cancelling all squares of the perturbation terms—-there 


follows 

P-pbd V? Ww? . 

=> (1 + tan- ax) = 

Ww, tan ax (tan ax — tan 6) — m(U — V) + 
Wa (1 + tan? @) (6) 


») 


Then by Eqs. (1), (2), (3), and (4), the pressure in- 


crease through the cascade is: 


P — p)/p] + (UL? 2) — (W?/2) (1 + tan? ax) = 
W{w)(1 + tan? ax) — 2u(tan ay — tan @)] (7) 


2.4) Determination of the Flow Field 


For the sake of simplification, Eq. (7) will be written 


as 


— vwcot B = b (S) 


«uy 


valid in the unstalled fraction of the cascade front; 
there is of course a very similar relation 


Wo — Ucoty =c (9) 


valid in the stalled fraction, the exact values of the 
constants c and y depending upon the assumption made 
concerning the energy transfer through the stalled 
blade channels. 

The analytical function w — w defined in the x < 0 
half of the x + zy plane by means of Eqs. (8) and (9) 
is given by (see Appendix B) 

btan 8 — ctan y + 2(c — 5b) tan B tan y 


tan 6 — tan y 


e- tur Zz _— 1 7 pb T 
Ae a} =| (10) 
re ‘3 


il the stalled zones are assumed to occupy the fractions 
(2k — v)xr<y < (2k + »)z;k = +1,2+2,...] 


of the cascade front and also if minimum singularity 
conditions are assumed. The parameter \ is deter- 
mined by the condition of regularity at the upstream 


infinity (w — iv = Oforx = —©@). 


2.5) Stall Propagation Velocity 


At both ends [y = (2k + v)z] of the stalled zones 
the R-angle of attack of the blades is necessarily equal 


to its critical value @ and there follows 
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tan@ = tana, b= Ucota — W (11) 


both relations being independent of the assumption 
concerning c and y. Thus, the rate of stall propaga- 


tion, 
V = U — W(1 + tan? ag)/(2 tan a) (12) 
and the pressure increase through the stalled cascade, 


(b — P)/p = (U?/2) — (W?/2) (1 + tan? as) — 
W2(1 + tan? ax) }[U/(Wtana)] — 1; (12’) 


(which derive from Eqs. (7), (8), and (11), without any 
assumption on the actual behavior of the stalled blades) 
give the following very simple expression for the rate 
of stall propagation at the onset point of the rotating 
stall (U’ = JW tan @): 











Fic. 2. S-streamlines of rotating stall 
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V = U}(1/2) + [(p — P)/pU?}} (13) 


Furthermore the R-angle of attack of the blades re- 
mains constant all along the unstalled blade front, the 
incidence of the flow being ‘‘saturated’’ to the critical 


value a. 


(2.6) Example of a Flow Field 


As an example of the flow field induced by the rotat- 
ing stall, Fig. 2 shows the S-streamlines corresponding 
to the operation of the mean section of the test com- 
pressor A, the performances of which are described in 
Section (4). 

In Eq. (9) c and y were taken to be zero—.e., no 
through flow was assumed for the stalled blade chan- 
nels. The following aerodynamic parameters define 
the flow: 


W/U = 0.47, tana = 1.78, 


tanas = 1 (»v = 0.14) 


The straightness of the downstream S-streamlines 
justifies a posteriori the use of linearization and in 
particular the assumption that the rotational S-per- 
turbation velocity remains constant along the lines 


y — x tan 6 = const 


The whole S-flow field can be described as made up 
of no-flow zones attached to the moving blades, some- 
what analogous to the potential flow in radial compres- 
sors where small flow rates induce high angles of attack 
on the blades and lead to dead fluid zones attached to 
the blades.’ 

(3) INDUCED FLOW FIELD IN THE CASE OF FINITE 
PERTURBATIONS 


(3.1) A Nonlinearized Representation of the Rotating 

Stall 

As the throttling of the flow increases, the mean R- 
angle of attack also increases and so does the extent of 
the stalled zones. Then the intensity of the pertur- 
bations becomes important and the above assumption 
of linearized theory no longer applies. 

The following assumptions allow then an extension 
of the linearized theory to the case of disturbances 
with finite amplitude, (w, v) designating at this time 
the components of the S-velocity. (a) No flow is as- 
sumed through the stalled blade channels,* the stall 
cell being a kind of constant pressure wake: 


w =O for -—yt<y<ver 


(b) This constant pressure wake is bounded at least 
at one side by a regular jet line along which the velocity 


remains constant: 


w =O at y= —vm and/or y= vr 


(c) The angle of attack of the blades in the R-motion is 
equal to its critical value @ at the regular edge of the 
front of the stall cell. 
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Thus, from the general definition of the local R-ang| 
of attack, 


tan ag = (% + U — V)/w% l4 


at the regular edge of the stall zone there follows 


Then from the Bernoulli equation, written along th, 
upstream S-streamline which eventually passes throug! 
this regular edge of the stall zone one obtains 


2) = (p/p) + [((V — U)?/2] (I 


(P/p) + (V? 


Since no pressure rise is given at this point by the rotor 
p designates the downstream uniform pressure —that 
is, the pressure in the wake. 

There follows the value of the speed of stall propa 


gation 


V = U{( ~P)/pl}} (13 


2) + [(p 


identical to the one obtained in the case of small per 
turbations, but valid here for finite disturbances also. 


3.2) Extent of the Stall Zones 


A one-dimensional approach gives then the extent 
of the stall zones since it can be assumed, by generaliz- 
ing the results of the linearized theory, that the R- 
angles of attack along the unstalled blade front are 
close to the critical angle @. Due to the conservation 
of the tangential momentum, there is no mean induced 
tangential velocity —that is, 

Y,m = V (17 
where the subscript m designates a mean value along 
the unstalled front of the blade row, where 


(2k + vy)r<y< (2k+2 —)2 
The mean axial velocity is very simply related to the 
rotor speed as 


Wo.m = U cot & (18 


Finally, the mass flow conservation equation yields 


the extent of the stall zones as 


vy = 1 — tan @/tana (1 


(4) COMPARISON BETWEEN THEORETICAL PREDICTIONS 
AND EXPERIMENTAL RESULTS 


4.1) Simplified Theory in View of Experimental 

Verification 

The mean conditions (17) and (18) at the unstalled 
front of the cascade imply also that the conditions at 
the front of the unstalled zones are independent of the 
extent of the stalled regions. It follows that the pres 
sure increase through the unstalled blade channels 8 
independent of the extent of the stall, and the stall 
speed of rotation remains constant along the whole 
organized stall performance curve. Thus for all e 
perimental verifications, the static pressure p neces 


sary for the calculation of VU’ should be measured 
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Fic. 38. Experimental compressor A 


during the unstalled operation, at the extrapolated inter- 
section of stalled and unstalled performance curves. 


4.2) Experimental Compressor A 


Compressor A is a one-stage single rotor with a 


casing diameter of 2R = 600 mm. and a hub tip ratio 
Ry R 


untwisted blades with chords proportional to the dis- 


0.625. The solid body blading consists of 
tance y to the compressor axis, the solidity being thus 
constant along the blade. The stagger can be char- 
acterized by a single angle ¢ and the stall always takes 
place along the whole blade. 

On such a compressor the measurements made at the 
V (R? 


closely the mean performances of the whole machine, 


mean radius R,, = + R,*)/2 represent very 
and in particular the local axial velocity is very little 
different from the mean axial velocity. Thus a two- 
dimensional cascade theory may apply at this radius. 

The experimental data are taken at three stations 
see Fig. 3). Station 1 represents the upstream in- 
finity where all perturbations are already damped, and 
permits correct measurement of the mean mass flow 
rate. Station 2 is situated immediately behind the 
rotor and gives the aerodynamic data (velocity, pres- 
sure and angle) of the flow at unstalled operation. A 
hot-wire anemometer placed in this section records 
during the rotating stall the periodic flow fluctuations. 
Station 3 represents the downstream infinity, where all 
downstream perturbations are damped out, but due 
to the absence of a stator, the tangential velocity com- 
ponent is not recuperated, and the pressure increase 
between stations 2 and 3 is due to the recuperation of 
the axial velocity component only. 

The complete operation of a single rotot compressor® 
is represented in Fig. 4, for a stagger angle of ¢ = 41.5 


to which corresponds the mean R-exit angle ax = 
30.9"). The performances of the rotor are character- 
ized by means of the following classical flow and pres- 
sure parameters: 


g¢=W/U, yp = 2(P3 — Pi)/pU? (20) 


where the reference velocity LU is the tangential velocity 
of the rotor at the mean radius. 

The point ¢ = 0.67 corresponding to the onset of the 
rotating stall is obtained in Fig. 4(1) by the intersec- 
tion of the extrapolated stalled branch of the perform- 
The 
Static pressure, measured only for the unstalled oper- 


ance curve (yg, W) with the unstalled branch. 
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ation, is represented by the static pressure coefficient 


¥. = 2(po2 — P;) pV? shown in Fig. 4(II). Its value, 


y, = 0.22 corresponding to ¢, determines the theoretical 


speed of the stall, namely 


V = (U/2) (1 + 9) 0.611 (21) 


The experimental values of | l’, obtained from the 
hot-wire oscillograms taken at station 2, are represented 
i*/ea~ 


is the period of the rotor given by a rotating 


in Fig. 4 (III); they correspond to |" 
where 7” 
contact mounted on the shaft, 7” the period of the stall 
and m the number of stall zones (7 = 1 or 2); the time 
basis tr = 0.01 second is also recorded on the oscillo- 
grams. 
The extent of the stall is represented on Fig. 4(IIII) 
where the theoretical value 
y= 1-9/3 22 
derived from Eq. 19 is compared to the experimental 


measurements vy = 7°’'/7” where 7”’’’ is the duration 


of the stall. 


4.3) Downstream Pressure Evolutions 


Behind the rotor the stalled and unstalled zones 
undergo a complex mixing process, which in a compres 
sor with large hub tip ratio is neither a constant sec 
tion nor a constant pressure mixing. In compressor 


A, where the tangential velocity component is not 


recuperated, the pressure increase between stations 2 
and 3 can be expressed by 

2 99 

Y 5 Y a <)) 


during the unstalled operation. During the stall, with 
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Extent of the stalled zones. 


Fic. 5. 


all conditions at the exit of the unstalled blade chan- 
nels remaining constant, the pressure variations at 
downstream infinity depend on the extent of the stall 
zones only; a constant section mixing would give 


¥-¥=—-(e— 9) (24) 
whereas a constant pressure mixing gives 
y-y=-2+°¢° (25) 


A more general mixing process, which would ensue 
from a combination of both elementary processes, could 


be described by 


¥—-y= —-ey — 2). - 
(l—e) (¢?— ¢?), O<e<1 (26) 

In the case of compressor A the straight performance 
curve of stalled operation [Fig. 4(1)] shows that both 
processes seem to have equal importance (€ = 0.5). 
It can be verified from Fig. 4(1) that the slope of the 
stalled y vs. curve is really close to ¢. 

In some compressors where the hub/ tip ratio is close 
to 1 (see reference 10) the mixing is nearly a constant 
section process and the value « = 1 applies to the 
downstream pressure recuperation. 


(4.4) Extent and Speed of Rotation of the Stall 


The experimental variations of » for various stagger 
angles (¢ = 41.5°, 45.5°, 49.5°, 53°) are represented in 
Fig. 5. The values of ¢ are obtained by correlating the 
stalled curves (¢, 
slope ¢. 

The correct determination of v is difficult due to the 


y) by means of a straight line of 


interference between stall zones and blade wake, exact 
values of vy being measurable only downstream of the 
rotor. In view of this fact, it may be considered that 
the divergence of the experimental points from the 
theoretical (vy, v) lines is not great. 

The corresponding values of y, and the theoretical 
and experimental V’/ U ratios are given in Table 1. 
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TABLE | 


4.5) Upstream Decay of the Disturbances 

The decay of the disturbances can be measured }y 
means of hot-wire recordings at various stations ly 
tween | and 2. 
compared, for one zone pattern, to the results of th 


In Fig. 6, the hot-wire recordings a 


calculations under the conditions of paragraph (2/ 
(transformed here into the equivalent electronic signa 
by means of the calibration curve of the hot-wir 
anemometer). This comparison is made by expressin; 
all upstream distances, measured from the axis of t 
blade, in terms of its ratio to the mean radius R,, of th 
rotor. 

This mean radius R,, is the characteristic length 
the decay of the upstream disturbances due to ther 
tating stall with one stall pattern. It is interesting t 
note that this length is somewhat greater than th 
characteristic length R/e, of the decay of the thre 
dimensional disturbances of a potential flow in t 
same annular space between the hub and the casing; 
the compressor, the eigenvalue of which is €, a functio 
of the hub/tip ratio only and generally close to 2x (s 


reference 11). 


(4.6) Experimental Compressor B 

Experimental compressor B, which is quite similar t 
compressor A, has also a solid body single rotor, but wit 
a hub/tip ratio of Ry/R = 0.79; the casing diameter i 
2R = 703 mm.; and there are 60 blades with a mea 
R-exit angle of 60 Due to the great number of blades 
































Fic. 6. Upstream decay of the disturbances. 
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sitates however the consideration of the three-dimen 
sional equilibrium of a flow containing uniformly rotat 






ing stalled cells, made up of nearly dead fluid. 













APPENDIX A-—-HARMONIC FLOW FIELDS RELATED BY 
MEANS OF A LINEAR COMBINATION ALONG THEIR 


COMMON BOUNDARY 









Let us consider two periodic analytical functions 





-asured by 






w — wand w’ — ww’, valid for x < 0 and x > 0, respec 
0 by 


ations bk 
rdings ar 
alts of th 


raph (2/ 






tively, vanishing at infinity and related along x 










Wy = awe’ + Db’ ’ Al) 









nic signal Let us assume, without loss of generality, that w’ 






hot-wir iv’ can be represented by Fourier series such as 































EX Pressing 


ixis of ¢] Fic. 7. Determination of the number of stall zones w' — 1’ = . (A,’ + 1B,’) A2) 
J] 1=0 
R,, of th 
the rotating stall appears always with a great number then the following relations result: 


length gf of zones, thus giving a less well defined periodic phe- 





to the 1. nomenon. wy’ = >> (A,’ cos ny + B,’ sin ny); 
resting t The determination of the exact number n of these 3 o’ = > (A.’sin nv — B,’ cosny) (A3) 
than the zones is obtained by means of two hot-wire anemom- ha : : 
he thre. eters placed at angular distances of 9 = 27/n’. Fig. 
Ww in i shows the hot-wire recordings obtained for flow uw = > [(a’A,’ — b’B,,’) cos ny + 
casing off coefficients ¢ = 0.46 and » = 0.48 with n’ = 3,4, and ‘ (a’B.' + b’A.") sin nv] (A4) 
i functio—f— 5. Itisseen that for g = 0.48 there are four stall zones ‘i 7 
‘oO 2a (se and for ¢ = 0.46 there are five stall zones. The speed os we SS etd’ = Re 
of rotation remains constant in both cases, VU = ” " a 
0.59, to which corresponds the value V/U = 0.61 de- i(a’B,’ + b'A,’) |e AS) 
rived from the performance curves. 
similart It can also be shown that for incompressible flows 3 s 
but with the usual similarity laws of compressors apply.'* This ; 
ameter is 'S obvious from the equations of motion and can be % 
1 a meng Verified by means of the experimental measurements >, 
of blades @ Shown in Fig. 8 (performance curves of stalled and un- a 
stalled operation of compressor B at different speeds 
of rotation; the stalled operations are represented by 1} 


means of solid points) 


(5) CONCLUSION ’ . 
N= 2200 rpm¢y 
above necessitates very few assumptions on the be- 1800 O 
havior of the stalled blade channels. The physical 1400 raN 
scheme on which it is based, saturation of the R-angle 1000 


The one-dimensional theory of rotating stall given 





a 


O 


Nol 


of attack (angles of attack relative to the rotor) and 
no through flow in the stalled blade channels, seems to 
have a good physical background. 

The predictions of some properties of the rotating . . 
stall—such as relation between the rate of stall propa- 





gation and the static pressure behind the rotor, extent 
of the stall zones as a function of the flow coefficient, 


we 








—a upstream decay of the disturbances, and downstream 
mixing of the flow—are well verified by experimental 
measurements, at least as long as the number of blades 1 ; 
is high enough to justify the assumption of an infinite 
number of blades. 
ae The aerodynamic theory of rotating stall verifies 0 , 
4 c: 3 4 2 ~ 6 


correctly the experimental results on stall propagation 


along a single rotor. The prediction of the stall neces- Fic. 8. Similarity rules. 
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APPENDIX B 

The transformation Z 

| (Fig. 9). 
linear discontinuous syst 


AB: 


or Eq. (S), and 


BA: Wo 


or Eq. (9) relates the 
complex velocity w — iv. 


the above linear system. 
AB 
BA 


1 Sears, W 


pp. 429-455, 1955. 


w= 


REF 


2 Emmons, H. W., Pearson, C. E., 
pressor Surge and Stall Propagation, Transactions of the Ameri- 


= ¢ 


em 


— % cot po = 


The constant 


b tan B — c tan y sin [(vr — y)/2]) ‘ 
Wy = - + AX cos 7 ) - | : fl | 
tan @ — tan 7 lsin [(vr + y) 2]! BT 
(Di 
b tan 8 — ctan 7 jsin [(v — vr) 2], 
w= + A cos B ; 
tan 6 — tan y isin [(y + vr) /2]I 
ERENCES fee 


Y¥ cot y = 


THE HILBERT PROBLEM 


C 


b 


and Grant, H. 


half-plane R(s) < 0 on the interior of the circle 


i 


-s-/ i 
8 
/@\ 2 
—_ 4 AY ) 
+ 4—--— +—~—_+—_- + 
| [- \ x 
P » 
# SL 
Pd 
L- 
Fic. 9 
vw = > [(a’B,’ + b’A,’) cos ny — 
0 
(a’A,,’ — b’B,’) sin nv] 
(A6) 
> [a’'(B,’ cos ny — A,’ sin ny) 4 
F 0 
b'(A,’ cos ny + B,’ sin ny) | 
Thus, Y% = —a’'ty’ + b’w’ (A7) 


P represents the upstream 


PA = 


The images of the stall zones are the arc AB 
= 2vr. At the boundary x = 0 of the half plane, the 


(B2) 


real and imaginary parts of the 


An evident particular solution is obtained by solving 


R., Rotating Stall in an Axial Compressor, Zeit 
schrift fiir angewandte Mathematik und Physik, Vol. 6, No. 6, 


Com- 
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Fy = Wo Vo = 
b tan 8B — ctan y + i(c — 5) tan B tan 7 


tan 8B — tan y 


which is obtained, satisfies on the whole circle t 


conditions (B1) and (B2). 
conditions results from the superposition on this py 


The flow defined by th 


ticular solution of the more general solution of t 
homogeneous system (c = 6 = Q). 


In this case, the analytic function F(Z) = w — 
defined by the condition that its argument tar 


(—v/w) be equal to —y on AB or —8 on BA. 7 
basic function of which the imaginary part is dise 
tinuous in A and B is log [(B — Z)/(A — Z)| andi 
value at the origin is 27v7 whereas its imaginary par 


takes the value (v + 1)z0n AB and v7 on BA. 


In the combination G = FE log [(B — Z)/(A — Z)}- 


1H where FE and H are constant coefficients, the correc fi 


values —y and —8 along the arcs AB and BA are 


tained if 


E=(6@—-—y)/x, HA 3 Eker I 

Thus, 
‘ B-—Z Bg 
= ( 
A Z. 

is a solution of the homogeneous problem. Of cours 
such a solution is not unique and any real functior 
regular inside the circle, could multiply e°. How 
ever, such a function is necessarily singular in A a1 
B, and is a real function \ (see reference 15) of tl 


"™(B — Z)/(A Z), the argument 


to the above calculations, is altern 


variable ¢ = e 
which, according 
Thus the general solution 


tively equal to 0 and z. 


the problem is 
+ rA(E)eO% Bi 


But generally, 0 < (y — 8)/m < 1/2 (it can be seen] 
in Section (2) that 8 = @ and y = 7m 2 in the particul 
case when no through flow is assumed in stalled bla‘ 
channels). If infinite velocities are assumed at th 


front of the stall (point B), which seems to be a corret! 
assumption, the solution with minimum singular 
corresponds to A = const, which gives Eq. (10). 

In particular, along the blade front (x = 0), 


of Mechanical Engineers, Vol. 77, No. 4, 


can Society pp 
169, 1955. 

3’ Stenning, A. H., 
Journal of the Aeronautical Sciences, Vcl. 21, No. 10, pp. 71 
712, October, 1954. 
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A Simple Method of Matric Structural 
Analvsis: Part II—Effects of Taper 
and Consideration of Curvature 


BERTRAM KLEIN* 
National Rocket Corporation 


SUMMARY 


The matric method of structural analysis presented previously 
in Part I is extended to include the effects of panel geometry and 
gage taper, stringer taper, and a consideration of curved elements 
The resulting matrices are of the same form as the previous matrix 
ind of comparable simplicity. Several problems are presented to 


illustrate the mechanics and simplicity of the extended method 


SYMBOLS 


| effective cross-sectional area, or area 
a,b shear panel dimensions 
D point in shear panel 
E,G = Young's and shear moduli 
F shear force 
5, f = functions 
h = rib spacing 
c = length of centroidal line of curved element 
VU bending moment 
P = axial load, or load 
panel shear flow 
R = radial force 
= radius 
S = stiffness 
coordinate 
distance 
1 = tangential force 
u, t = displacements 
x = coordinate 
8 = angle subtended by curved element 
= parameter 
A = determinant 
6 = angular rotation 
k, A = parameters 
é = effective area parameter 
o = stringer stress 
T = shear stress 
d = (8/2 
v = function 
Subscripts 
( = centroidal 
_ a = numerical indices 
max = maximum value of 
min = minimum value of 
0 = outer 
= panel 
= due to shear flow 
= rib 
tt, ete. = denoting various stiffnesses 


Superscripts 
= average value of 
= effective value of 


Received March 12, 1957 
“Formerly, Chief of Stress, National Rocket Corporation 
Now, Head, Structures Methods Group, Convair (Astronautics) 
Division of General Dynamics Corporation 





INTRODUCTION 


A PREVIOUS PAPER,' the fundamentals of a method 
of matric structural analysis have been presented. 
In the present paper, the work is extended by the deriva- 
tion of the equations necessary for treating tapered 
structures and structures containing curved elements, 
both of which are found extensively in practice. 

The aim of this series of articles is to present a method 
of matric structural analysis in which the setup routine 
is reduced to such simplicity that the structural engi- 
neer can proceed directly and rapidly in the initial 
The 


task of obtaining final numerical values is assigned to 


stages of the solution of a complicated problem. 
the computer engineer. The size and form of the 
matrix therefore is not the direct concern of the struc- 
tural engineer, although he must be cognizant of the 
Jimitations of the science of automation. 

In the proposed method, once a structure has been 
idealized, the necessary final equations of equilibrium 
and compatibility are written down directly and 
rapidly in final matrix form without recourse to any 
artificial devices (see numerical examples). In many 
applications, it is of importance to obtain a complete 
both deflections and 
for this requirement 
Also, with the in- 


solution to a 
Allowance 


problem—4.e., 
stresses. is made 
directly in the proposed method. 
creasing importance of thermal loadings, and combined 
thermal and regular loadings, an efficient method should 
be able to treat such problems with no additional com- 
plexity. Allowance is made for these cases directly in 
the proposed method. 

It is admitted that the proposed method is not in- 
tended to solve all structural problems, or lead to 
mathematically exact solutions, or is the best method to 
use in all cases. However, it is believed that, in a 
variety of problems, the method should prove to be 
efficient. 


(1) TAPERED PANELS AND STRINGERS 


It has been learned through experience that, in many 
practical cases, it is permissible to idealize a tapered 
reinforced sheet structure by assuming the structure to 
be composed of two types of elements—panels taking 
only shear, and reinforcing elements, assumed to bound 
the panels, taking direct stresses. Therefore this as- 


sumption is maintained. However, when a shear panel 
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A Min’ 


Fic. 1. Shear panel tapered in plan form and thickness 


is tapered in plan form, such as the one shown in Fig. 1, 
the shear flow in the panel is found to vary in order to 
maintain a state of equilibrium. The magnitude of the 
shear flow is largest at the smaller end of the panel, and 
denoted there aS mar. The panel is assumed to have 
two parallel edges. The shear flow is assumed to be 
constant at any section parallel to these parallel edges 
i.e., fora constant value of x. The thickness of the panel 
is assumed to vary approximately linearly, increasing 
in the x direction, and to be approximately constant at 
any section parallel to the parallel edges. 


(a) Force-Displacement Equations 

For the assumptions made, the necessary expression 
relating the deflections of the four corners of the panel 
of Fig. 1 to the shear in the panel is [see Appendix (1) | 


(M1 +. U2)Ay—2 + (V1 + V3) Daa 
(v2 + v4)Omin — (Uz + Us)ds—4 = (2A,/Gtp)g (1) 


where g is the average shear flow in the panel and is 
J = Ymar(Omin/ Omar) (la) 


t, is the thickness of the panel at the point D, located at 
the intersection of the diagonals of the panel as shown in 
Fig. 1, and A, is the area enclosed by the panel. Note 
the inversion of the multipliers of the terms containing 
the v displacements. 

A stringer is assumed to be tapered so that the effec- 
tive cross-sectional area increases in the x direction, as 
shown in Fig. 2. In accordance with the assumptions 
made regarding the state of shear stress in the adjoining 
panels, the desired relation between end displacement 
in the direction of the stringer, the force at the smaller 
end of the stringer, and the average shear flows in the 


adjoining panels is [see Appendix (1) | 


(EA/a)(u; — Uj41) = 
Py + @ Dd def ((Omin/Omaz)» (Amaz/Amin)] (2) 


n 








i 5 | 
R ce SS I Pi, irs 


where A is the harmonic average value of .1— ice, 


and f is a function dependent on the type of stringy 
taper. 5 is the value of s at which A = A. For 
ample, the values of f and § if the variation of 
linear with respect to s are shown plotted in Fig 
The value of 5/a is read from the vertical where (4, 
Dace) = 1 on the right. Thus if A,,-/Acu 2 at 


mar 


) 


the variation is linear, 5 a = 0.44, from Fig. 3. 7 
positive sign for g acting on the stringer is that show 
in Fig. 2. The value of ” usually does not exceed 
Curves similar to the ones presented in Fig. 3 may} 
drawn for other types of stringer taper variation, su 
as quadratic, etc. 

Since it is assumed that the shear flow is constap 
along a rib, the expression relating end displacement 
and end force is of the same form as the equation give; 
previously in reference 1 (Eq. 1). For example for rij 
1-3 of Fig. 1, 


(2EA,/bmin) (01 — U3) = Pi + Ps 


(b) Equilibrium Equations 

For a stringer such as shown in Fig. 2, the equilibriun 
equation is simply [see Appendix (1) | 

| re nas lar mei Dan 

For a rib the desired relation is of a similar form. For 

the ribs shown in Fig. 1, for example, one may write 
Ps— Pr= Qo bngn, P2 — Ps = DL dnd, 
nN n 

The correct sign of each average shear flow acting on the 
stringer or rib must be used in Eqs. (4) and (5). Not 
here the inversion of b,,;,, and 0,,., such as occurs in E 
tt). 
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Fae. 3: Tikbams ices), (Awesr/ A )} and $/a for linear stringet 
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4 SIMPLE METHOD OF MATRIC STRUCTURAL ANALYSIS S15 
——— (c) Numerical Example 
F SO: EL: AEE AO. AS sacl ; 
7 ; 8 Consider the plane sheet-stringer-rib tapered panel 
4 depicted in Fig. 4. For purposes of simplicity, the struc- 
884 4, 8 m ture is assumed to have an axis of symmetry, the loading 
is taken to be symmetrical about that axis, and the 
4\_\ ESS structure is assumed to be completely fixed at the sup- 
i \. 76 —_—_-9y--— 95 ia 73 ae . : 
d - 7. / ports. The indicated panel dimensions are numbers 
884 _ 8 which could be read from a drawing. The necessary 
q ee values of the parameters appearing in the force-dis 
2 \ oy / placement equations are taken as 
a SS 
--375-+\— 5 —4— 5—- 4 a 9 9_* 9 . 
Lies stringer 2 20 5 a 16 
' qu t - 
c = ¢ ; O.S ) 
RB B EA/a l 1.5 1.2 
Fic. 4. Geometry and loading of tapered panel rib 1-3 2-4 1-1 2-2 
2EA,/b » t ) 7.5 
All the necessary expressions needed for solving panel l 2 
tapered quasi-flat reinforced panel structures are now 2A ,/Gtp 100 100 


known. It is of interest to note that the coefficients of 
the unknown quantities appearing in the above equa- 
tions are panel dimensions or quantities that may be 
read directly from drawings of a structure or auxiliary 
curves (such as Fig. 3) or readily calculated from such 
data. 
to illustrate the workings of the method. 


A simple numerical example is presented below 


l 0 0 000 8.84 0 

| l 0 000 0 8.84 

0 0 l O00 S.1 0 

0 0 l 100 0 8.1 

0 0 0 010 7.5 0 

0 0 0 00 1 5 10 

0 0 0 000 5.04 0 0 0 0 

] 0 0 000 0 1 86 0 0 0 

0) 0 0 000 +. 62 0 I l 0 

0 oO 1 000 0 —~4 46 0 1.5 0 

0 ft) 0 O11 0 0 0 0 0 

0 0 0 O00 ] 0 0 0 0 0 

0) 0 0 010 0 a) 0 0 0 

0 0 0 00 1 0 0 0 0 0 
100 0 8.1 8.1 —8 

0 — 100 0 §. 1 O 


The solution is 
P, = 0.60892090 P; = —0.88963159 
P, = 0.46945728 fe = 0. 12282344 


—0.056117545 


P; = 0.45455211 a = 
—0.015776428 


Py = 0.58234118 i 


As regards the check of equilibrium, it is interesting 
to note that equilibrium at any section in the direction 
of Po checks almost exactly provided it is assumed that 
the stringers are inclined such that the angles of in- 
clnation with the direction of Py are exactly arcsin 
‘ 8.84) and aresin (8/8.1).. The check equations for 
this assumption are 


(S'8.84)P; + (8/8.1)P3 = 0.999999992 


(8/8.84)P. + (8/8.1)P, = 


S84 


FG sn! Oieed, (Aen | Ais | 0.87 0.55 


The external load, Po, is taken as one, and broken up 
into two components, one acting in the direction of 
stringers 3-4-6 and of magnitude, 8.84/8 = 1.105, and 
one acting in the direction of rib 3-1-1’ and of magni- 
tude 3.75/8 = 0.4685. The matrix defining the prob 
lem in the 16 unknown quantities may be written 


P 1.105 
P, 0 
P 0 
P, {) 
P 0.4685 
P,; 0 
0.8 q 1.105 
= q 0 
0 My 7 0) 
0 Us 0 
0 5 0 5 0 u 0.4685 
0 0 —7.50 4.5 u 0 
0 5 0 0 0 0 
0 0 7.50 0 0 
—S8.S84 7 5 5 r 5 5 () 
—S8§_ 84 0 10 0 10 7 0 t 
u,; = 0.60920638 “1 = 0. 17792632 
u. = 0.34994332 v = —0.016376459 
us = 1.4712490 V3 0. 26210263 
Us = 0.44353955 vy = —0.032752919 


On occasion, it may be of benefit to introduce extra, 
actually nonexistent, stringers and ribs, thereby sub- 
dividing actual panels into smaller panels with conse 
quent increase in the size of the matrix. However, this 
artifice does not change in any way the basic equations 
or mechanics of the method. (This procedure is the 
opposite of the more customary method of “lumping” 


in which the size of the matrix is reduced when existing 
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are considered. Refer to Fig. 5 for the basic geometry 
g etr 


stringers and ribs are combined together thereby in- 
creasing the sizes of actual panels. and loading of such an element. Notice that each end 
the element has associated with it three generalized fore, 

(2) CURVED ELEMENTS and three generalized displacements (henceforth calk 

‘ ; ss , simply forces or displacements). Consequently, ther 

In this section, uniform (untapered) elements of con- See is 3 é me: 
exists for the element three equations of equilibriy; 

stant curvature and acted upon by constant shear flows : : ; ai 
: and three force-displacement equations. These equ 


tions are 


Equilibrium: 
sum of vertical forces: (T; + Ti41) tan @ + (R; + Riv1) = O j 
sum of horizontal forces: (7; — 741) cot @ — (R; — Ri+i) = 297, 
sum of moments about 0: (7; — Ti+1) + (1/r-)(My — Mi41) = (7,/r-)29r, oh 


where 7, is the radius of the outermost fiber of the element, and r, is the radius of the centroidal line of the elemen; 
(see Fig. 5); and force-displacement [see Appendix (2) |: 
(M, + Misi) = (65 — 0:41) Sim + (0: + 9441) Simr + (Ui — i441) Sint 
(R; — Ri+1) = (0, + 0411) Sim (vo; — 0441)S,, + (Us + Ue41) Sr — 2K, 
ee + i +1) T= (6; — 644 he —- {Py + v i der + (i; —= Us ites 


-_ + 


where R, is the radial force acting on the left end of the element in Fig. 5 when shear only is considered—i« 
u=v=60= 0. 
The values of all the stiffnesses, S, may be expressed matrically in terms of but two stiffnesses—i.e., S,, and S,, 


as follows [see Appendix (2) ]: 


Sa = (2) Ly) Sar Olas r.{(tan @/¢) — 1]S,, 0 0) r,[(tand/¢) — 1] — tan®@ l | 
be i Jt 1 0 a 0 —r,tan @ l tan @ § 
oie a 0 Oo Sa r-{(tang @¢) — 1] — tan¢d 


The values of S,, and S,;, may be calculated from the equations 


S_(r3/ EDT) = 2(1 + cos B)/(xB — X sin B) 
' rape . ; i) 
Sil(re? ED) = 2(1 + cos B)/[k8 + A sin B — 4(1/B)(1 — cos B)] 4 
where 
K«=14+(1/y)(18 +1), ’=1- (/y)[0/9 -— 1], y = Ar2/], & = A*/A (9a 
A = axial cross-sectional area of element, J = centroidal moment of inertia of cross section of element, and A* = 
effective shear area of cross section of element (see reference 2). The value of X, is given by the expression 


: [x(8 — sin 8B) + B(1 + cos B) — 2 sin B] 
«G6 — AX sin B 


R, = 4 (10 
provided the effect of eccentricity may be neglected—i.e., e 7 is very small. The more correct expression including 
eccentricity is given in Appendix (2). 

Frequently it is permissible to take y infinitely large so that x = A = 1. 

It is convenient to rewrite Eqs. (9) in the following forms: 


2(ET /L*)(1 + cos B) 


S; (9a 
[(1/6) — (1/5) 6? + (1/7)B4 — ...] + U/AL?) (4 /A*) [1 + 


18) sin B] + [1 — (1/8) sin B]} 


2(EI L*)(1 + cos B) 
(1 B7)W(B) + TAL?) [((4/A*) p(B) + xr(B) | 


= (9b 


. 2(FA/L)(1 + cos B) 
S., = bs (9c 


(AL? 1I)5(B) + (4/A*)W(B) + Y-(B) 


where (8) = >) (—1)"*+182"[2n/(2n + 4)!] (9d 
n 1 
= 67[(1/360) — (1/10080)6? + (1/604800)64 — .. .} (9e 
Eq. (10) may be written as 
R, = qr\(1 2)¥i(8) [1 + (S,,/St)] + £(8)(S..L3/ED } (10a) 


where /(8) is given in Appendix (2). 


i 
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Infinite 


If 7 is taken to be infinitely large, the above equations become 


5 2( EI L*)(1 + cos B) 2(EIT L*)\(1 + cos B P 
O- = P (y") 
(1 /6* W(8) F(B 


A very close approximation to X, is given as 


R, = qr,[0.18" — (e/r.)(1 + 0.48? 10’ 
for infinite y and values of 8 less than about one radian. This expression includes the effect of eccentricity 
Straight Element 

When 8 = O, the element becomes a straight element, and 
S,. = (2EI/L), S,, = — (2/L)S,,, = (24EI/L*)/({1 + 12(7/A*L?)|, S, = 2EA/L\ r 
R, = —eq/[1 + 12U A*L?*) { ; 
All other S values become zero. Eqs. (6) and (7) become 
R,+ Rai = 0, T, — Tiwi = gh, M, — Misi t+ LR; — Ri+i) = 2eqgl 6’ 
M, + Misi = (0: — 0:41)(2EI/L) 

(v, — Vit1) — (L/2)(0; + 0441) (2447 /L*) + 20eq 

R, — Rui = * | = 7’ 
1+ 12([/A*L?) 


:. | (ere = (u — uj41)(2EA 5) 


All the necessary equations for treating untapered curved and straight elements taking combined loading in their own 


planes and under constant shear flow are now known. 


Numerical Example 
For purposes of analysis, only one 


Let us consider the circular ring loaded antisymmetrically, as shown in Fig. 6. 
2, and 2’-3, referred 


half of the ring need be considered. This half is broken up into two curved elements, segments | 

to as | and 2 hereafter. The ring cross section is assumed to have the following properties: 
A = 0.250, J = 0.1666, ¢ = 0.16, EI = 10°, e/r =0 

The values of « and \ are computed from Eq. (9a) tobe « = 1.0120785, A = 0.991253. 

The following values of the stiffnesses times 10~* have been computed using Eqs. (8) and (9 














oe Ms eT t 214.4147079 — 19. 78466685 8. 195077336 
ae \, So He = — 209982022 0 .0434887007 0. 104991011 
Sim Sn, Sm | = 8. 195077336 =—0.942802111 0.390521421 | 
eT [ 251.8798715 — 47 .0855601 113.674598 
ae a Sits = — 37 .6270642 $.5419SS44 1.SS8135321 
5, te et | Beene —42.9721225 103.743881 
The following values of 2(X,) have been computed from Eq. (10 2(R,)1 = 25.26276318, 2(R,)2 = 2.770052856. 


The ring is assumed to be fixed at point 1. Twelve unknown quantities are found to exist. These unknowns consist 
of seven forces—R}, M2, R2, T:, Ry», Tx», and R;—and five displacements—2, v2, U2, 93, and U3. 


i.e., Eqs. (6) and (7) written for each element. The matric equation de- 


The necessary twelve 


equations consist of six for each element 


fining the problem is 


tan ¢ 0 1 1 r’, U 
—cot ¢; 0 —l1 |] VW 2qr 

—] —(0.05 0 O Ry 29rd, 
0) 0 0 0 tang 1 | R, 0 

0 0 0 O cot ¢d —1 l T 2qr 

0 0.05 0 0 l 0 O R 29qr@ 
0 1 0 O 0 0 8s — S,., . R 0 

0 0 —1 1 0 0 .66hCUS --S., Sy 0. 2(R,)1 
0 0 0 0 . 6 &. -S,, 3 v 0 
—] 0 0 0 0 8 Ss. Sn Sun ~Sum Sus Us 0 
0 —] see ie aa a Sits As 2( R, 
—! 0 0 S, % Ss ~S, ~Se U3 0 

J L L J 
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Fic. 5. Curved element geometry and loading. 


Let P = 100 and q = 1.1253953925. The solution is 


T. = —28.89457266 R; = —36.62899535 
M, = —482.768716 = —2.91185135 
kK = 18.3551S868 ie = 94 5672443 
ke = 51.4024805 UW. = —93.2107074 
ie = $1.81610529 & = 1. 23966596 
| 19.30819743 us = —137.333729 


The check of vertical equilibrium of half the ring is 


R, — R3 + 2gr = 18.35518868 + 36.62899535 + 
15.01581570 = 99.99999973 


The calculations are repeated assuming that x = A = 
I—i.e., assuming y to be infinitely large. The final re- 
sults are found to be 


T, = —28.93162 Rs; = —36.57660 
M, = 482.0277 6. = —2.947629 
Ri = 18.40758 ™ = 91.460459 
R, = 51.43953 uw = —93.66575 
Ty = 41.77905 i, = 1.1977498 
Ry = 19.27114 us = —136.1715 


The resulting forces are seen to differ but little from the 
forces calculated previously, and the deflections are of 








Fic. 6. Circular ring under antisymmetric loading, 
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approximately the same magnitude as those found be. 
fore. 

Certain problems involving three-dimensional reip- 
forced cylindrical structures containing curved ele. 
ments and curved panels may be solved using the data 
derived herein for curved elements and allied equations, 
However, for certain problems, such as flexible frames 
loaded so that the induced panel shear flow varies 
radically over relatively small arcs, the method pre. 
sented here may have to be modified. 


APPENDIX (1) 

First consider the derivation of the force-displace- 
ment for a shear panel simultaneously tapered in plan 
form and thickness in one direction (x) as shown in Fig, 
1. The required relation may be derived by using an 
elastic energy principle as follows. Equate the work 
done in shear by the four external loads, F\—», Fi—s, F; 
and F;—4, to the energy stored in the panel due to shear 


(F, o/2)(1/2) (41 + ute) + (F; 3/2)(1/2)(%1 + v3) 


s)(1/2) (tg + m4) = 


(F, 4 2)(1 2) (ve + tay = (Fs; 
(1 2G) | TOUS = (GnaxOmin/ OG) He 
| (1/t)(Omin/b)8dx  (A-1 


0 


since 7 = (g/t) = Gmaz(Omin/0)*(1/t) 


Substitution of b = by in + (Omar — Omin)(x/h) and in- 
tegration yields, for constant f, 
[Gmaz?/(2Gt) \(/2) X (bmar + Smin)(Omin/Omaz)? = 
[g?/(2GH)]A, (A-2 
When ¢ varies approximately linearly in the x direction, 
the following close approximation applies: 
(g2/(2Gt,) ]A, (A-3 
where ¢, is the panel gage at point D, as shown in Fig. | 
at the intersection of the diagonals of the panel. 
The following examples are given to illustrate the 
accuracy of this approximation. 
First let (0,../Oniu) = 2, and Upec/tein) = 2. Te 
correct value of the integral is 


h h 
[ (icin /t)(Onsn/O) dx = [av [1 + (x/h)}! = 0.292h 
70 0 


In the approximate method, the thickness ¢, occurs at 
(x/h) = 0.31 in this case, and the corresponding ap- 
proximate value of the integral is 
(h/2)3(1/2)?/(1 + 0.31) = 0.286h 
in error by about 2 per cent. 
If (Omar/Omin) = 3, and (tmar/tmin) = 2, the results are 


by direct integration, 


| dx —— 
o [1 + («/h)][1 + 2(x/h)}? 


and by the approximate method, using the thickness at 
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0.25 in this instance, Refer to Fig. 2. 





panels. The axial strain at any section 


; ao ‘ a located by the coordinate s is 

h/2)4(1/3)2/(1 + 0.25) = 0.178h y the coordinate s is 

in error by about 3 per cent. du/ds = o E = (p a > tals) LE (A-5 
Since the tapers in the latter example are extreme, the olin lias 

approximate method is considered accurate (as long as Integration yields 

the thickness variation is approximately linear and only 


. : eal” iw *s 7 
in the x-direction). Uu,; — U1 = | | (. + pe | ads AEs (A-6 
J0 J0 


Return to the evaluation of Eq. (A-1). The forces 
appearing on the left-hand side of the equation may be Each q may be expressed in the form 
expressed as follows by taking moments about appro- 
- GQ = ¢ Is. /(s + s)/° ,.= [(h h — | 
priate points: f {maz | ti i ai C/N 
) « 
A-i 
F, 2, => F, 3(1— 2 b, wr) = GQ-»2 2 : : : : ; 7 
rp F,_.(b h ab (A-4 s; 18 the distance from the intersection of the sloping 
‘o—4 = £1—-3\9min/ mar) = YG Omir saa) A 7 . ‘ 
ing « B-chte-a/lnes) © Ses edges of the element to point 7. Integration yields 
where | gds = Qmar|SiS/(S¢ + 5 (A-8 
/ V0 
F, 2 * e2 n = Qn vel ie . ae z= g b, - = : P ‘ ; ; ; 
Final evaluation of the integral in Eq. (A-6) results in 
Substitution of these results into Eq. (A-1) and division a _ 
(EA/a)(uy — ini) = Py + ad, Graf, (A-9 


by g 4 yields Eq. (1) of the text. 
Consider next the derivation of the force-displace- — <— A isk si 
ment relation for a tapered stringer bounded by m shear i the variation of 4 1s linear with s—1.c., 4 


= wn t+ (Amer — Amt (S/@Z), 
a, ee j | a= Amin) — ‘| log (Omaz/Omin) ) 
p= - 41 —- ‘ ; (A-10 
[ (Omar Buen) — (Awa And ' (Ones bmin) — 1 Jjlog Mec! asta { 


There remains to discuss the derivation of Eq. (4) of the text. This result follows directly from the relation 


ea 
| ds = ag 
/7 0 


obtained by letting s = a in Eq. (A-8) above. 


APPENDIX (2) 


Eqs. (7) of the text and the stiffnesses, S and R,, appearing therein may be derived with the aid of the concepts and 
equations presented in reference 2. Using the notation of reference 2, for example, one may write 


T, = Udty + vdty + Odny — Ussilty + Viitrr — Ou itnp + 1, ( (A-1] 
Tit, = —Uipitty + ve-atry — Oi4atny + Udte + vdrr + Otne — 1,5 ; 
Adding, one obtains 
Ti + Titi = (41 — Ui+1) (thy + tte) + (Ue + Viti) (7 + tre) + (6; — 0:41) (tmy + inp) (A-12) 


Therefore, the following relations hold: 
S,, = tty + ttr, Sy = tty tire, Sm = tty + tn (A-13 
In a similar manner, the remaining two equations of Eqs. (7) may be derived and the following relations established 
mm = MNy + ntp Sm = Ny + ntr Sm = hy — Hert ‘14 


Sm = My + Nr; S,, = %y — rr St = thy — rtp f 


The task of the determination of the actual value of S, therefore, resolves itself into the evaluation of the sums (or 
differences) of the quantities, /f,,, ftp, etc., presented in reference 2. The complete development of these evaluations 


islengthy and will not be presented here. First one eliminates the quantities with subscript, F, by writing the matric 


relation [from Eq. (34) of reference 2]: 
ln ae. See 2 r. sin B r-(1 — cos 8) | NNy Ny nly 
Ss. S.. S&.1=1)0 (1+ cos #) sin 8 fie We Te A-15 
Soin ey es 0) —sin 8 (1 + cos B Iny thy thy 


Next, the quantities with subscript .7 are expressed in the form of the ratio of the determinants given in Eqs. 22 of 


reference 2 
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(A/r.?) = (8/4)(Bx — Asin B)[B« + Asin B — 4(1/B)(1 — cos B j : eB 
(r./EI)(—A,, sin B + A,,(1 + cos B)|] = (8/2)(1 + cos B) (Bx — Asin B A-l 2 
(r,/EI)[A,,(1 + cos B) + A,, sin B] = (8/2)(1 + cos B)[B«k — A sin B — 4(1 B)(1 — cos B) \ 
By taking the necessary quotients of the above re- 
sults, Eqs. (9) are derived. In a similar manner, the 4 
relations among the S’s given in Eq. (8) eventually may 
be found. The equilibrium Eqs. (6) also prove of use , was 
- os ‘ where e = (ry. — r,.) aSS "i tig. 6. ; 
in establishing some of the relations. For example, — 0 -) as shown in Fig. 6, and 
from the first equation of Eqs. (6), it follows that 83f(8) = @ + }—1 + (1/2)~1(8)(tan o/¢)!tan ¢ = 
: . , D0? | ( + ( 3780)B? + (323 /72576 )B 
‘oa <S as °[(1 1440) + (1/3780)6? + (3237257600) 8 
= - > ‘ 7 F — ” o(| = (| So”) | (te ‘ating | — ’ 
Che value of RX, is derived with the aid of Eq. (19b) of g() | Sg’) [(tan ¢ ¢) L | 
7 Sar ‘ : (1/2 + (1/60)82 + (17/252)@4 + 
reference 2 (with the right-hand side set equal to zero) 1/24) + (1/60)6 ly 202)6 
and the equilibrium equations (6) of this paper. The When ¢ is set equal to zero, the expressions given in Eqs rl 
derivation is lengthy and not given here. The result, (10) are obtained. rbit 
including the effect of shear flow eccentricity, is unst 
; REFERENCES — 
R, = gr, K(B — sin B) + (7/re) X ee 
[81 oe cos 3) — - sin B]t (xB one r sin 3) (A-17) ' Klein, B = A Simple Method of Matric Structural Analys ilong 
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Diffraction of Disturbances Around a Convex 


Right Corner With Applications in Acoustics 
and Wing-Body Interference’ 


LU TING* 
Polytechnic Institute of Brooklyn 


SUMMARY 


The problem of the diffraction of a small disturbance or an 
irbitrary pulse by a convex right-angled corner is reduced to two 
unsteady problems in half space. The matching of the boundary 
conditions of the solutions to these two problems leads to a non- 
homogeneous integral equation for the pressure distribution 
long the boundary The solution is obtained by the method of 
successive substitution. It is shown that the series given by this 
procedure converges faster than the power series of one half. The 
convergence of the series is accelerated by regrouping the terms 
so that the solutions corresponding to any one of the ordinal 
numbers of substitution possesses the correct behavior at the 
boundary of the domain. 
shock, 


the solution corresponding to the second substitution is already 
I 


When the given disturbance is a plane weak (acoustic) 
in good agreement with the linearized conical solution. For the 
other two nonconical fundamental problems considered in this 
paper, the series converges so rapidly that there is no necessity to 
go beyond the second substitution 

By superposition of these fundamental solutions, the pressure 
distribution on a rectangular barrier after it is contacted by a 
weak shock is obtained. The problem of interference of a pris 
matic body of rectangular cross section with a planar wing in super 
sonic flight can be reduced to the problem presented in this paper 
ind the one amenable to solution by Evvard’s method. The 
pressure distribution is obtained for four different wing-body 
configurations. The value of the integral of the pressure dis 
tribution agrees with the value obtained directly by means of a 


certain integral relationship 


SYMBOLS 


1, Ao, Ay, A = constants 

B =VIi?- 1 

c = speed of sound 

G defined in Eq. (11 

HH,“ defined in Eq. (23) 

2L = width of the barrier or the body 

V = Mach Number of undisturbed stream 
N = an integer, 1, 2, 3, 4, 

R defined in Eq. (28 

S = incoming disturbance 

di, 52 = additional disturbances 

g defined in Eq. (13 

ho, hh, hh, h defined in Eqs. (12), (13), (18), (19), (20 
n = an integer 0, 1, 2, 3, 


= pressure 


= time 
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> = Cartesian coordinates with ¢ in the direction 
of the free stream 
= 3/B 
a, B, a 'B = characteristic coordinates 
A = Laplace operator 
t,7,¢ = variables in place of x, ¢, and z, respectively 
p = density of the free stream 
€ = pressure jump across a weak shock 
6 = inclination of the upper surface of the wedge 


with reference to the free-stream direction 


INTRODUCTION 


motion and 


y I ‘iz FLOW FIELDS produced by the 
diffraction of infinitesimally weak (scalar 


dimensional waves are governed by the wave equation 


two- 


Age — (] C?) Cy = 0 


whereas linearized three-dimensional steady supersonic 
flows obey the differential equation 


Ace — B’o.. = 0 


which is equivalent to the wave equation. Therefore, 
problems in these two fields can often be solved by the 
same technique or may even be represented by the 
same type of solution. In this paper an analysis is 
given of the flow field associated with the diffraction 
which comes in contact 


of a two-dimensional wave 


with a convex right-angled corner. The solutions ob- 
tained are applied to the analysis of the equivalent 
steady flow problem of the flow field produced by the 
interference of wings and prismatic bodies of rectangular 
cross section. 

The diffraction of a two-dimensional wave by a wedge 
was analyzed by Sommerfeld.' Since then numerous 
investigations of diffraction problems have been re- 
ported. An extensive account of them can be found 
in reference 2. Much of the work cited deals with time- 
harmonic problems in which the incident wave is time- 
harmonic or “‘monochromatic.’’ For these problems 
the simple wave equation reduces to the time-free 
elliptic equation Ag + k*g = OQ. In principle, the 
solutions of these problems can be extended with the 
aid of the Fourier integral to produce solutions to dif- 
fraction problems in which the incidence wave is not 
time-harmonic. However, the inverses of these Four- 
ier-transform solutions have not yet been determined. 

For linearized three-dimensional supersonic flow, 
the above technique is not practical because the flow 


field to be analyzed is usually nonperiodic with respect 
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Fic. 1 Pressure distribution at an instant ¢ > ( 


to the coordinate (z) in the direction of the undisturbed 
flow. A different approach was introduced by Buse- 
mann*® who pointed out the frequent existence of conical 
fields in supersonic flow. In conical fields the govern- 
ing equation for the velocity components or pressure 


two independent variables. Busemann’s original treat- 
ment has been applied to delta wings by Stewart* 
and to quadrilateral wings with dihedral by Snow,° and 
has been developed in detail by many others (see 
reference 6). 

The general problem of lifting planar wings in super- 
sonic stream cannot be resolved into conical fields. 
Two different approaches to this problem were pre- 
sented in 1947. Gunn’ applied the method of Laplace 
transform. Evvard* ° expressed the solution in terms 
of the equivalent source distribution. The simplicity of 
Evvard’s method has resulted in its being preferred 
usually for practical applications (see p. 184, reference 
6). Nevertheless, Gunn’s method has certain ad- 
vantages for wings of low aspect ratio. 

Busemann’s method of conical flow has been applied 
by Keller and Blank!’ to the case of a uniform plane 
pulse incident on a corner or wedge. Diffraction prob- 
lems of this type correspond to the problem of wings 
with dihedral in supersonic flow.° 

The problem of a uniform plane pulse incident on a 
thin barrier (or an infinite strip) was analyzed by Ting?! 
and Miles!* on the same physical ground as the super- 
sonic flow over a rectangular plate of infinite chord 
length. In the treatment of reference 11, the time- 
history of the diffraction problem is subdivided into 
equal time intervals. At the end of each interval the 
disturbance originating at one edge of the strip reaches 
the other edge, and the response in each case is equiva- 
lent to the basic problem of diffraction by a half plane. 
This basic problem is solved by Evvard’s method and 
detailed pressure variations on the thin barrier (or the 
strip) are obtained. On the other hand, in reference 
12 the method of Laplace transform by Gunn’ and Fox" 
is used and the complete time-history of the total 
force acting on the strip is obtained. 


distribution can be reduced to a Laplace equation of 











In reference 11, it has been pointed out that the time 

history of the diffraction of a pulse by a rectangula; 
barrier (or block) can also be subdivided into tin, 
intervals. At the end of each time interval the dis 
turbance originating at one corner of the block reache 
the other, and the response in each case is equivalent 
to that of the basic problem of the diffraction of ay 
arbitrary pulse by a convex right corner. The approach 
to this basic problem suggested in reference 11 may be 
briefly indicated with the aid of Fig. 1. The incident 
pulse and its reflection by one side of the corner (the 
negative x-axis) is given by the pressure disturbance 
So(x, y, t), which is defined in the region y > 0 for ¢ < | 
and crosses the corner at the instant ¢ = 0. If the same 
expression of Sy) is extended to the same region y > | 
for ¢ > 0, a discontinuity of pressure Ap = So(x, 0,1 
will be created along the positive x-axis. This is q 
violation of the physical conditions of the problem 
The discontinuity is removed by the superposition oj 
two addition disturbances S; and S» created at the 
corner at the instant ¢ = O for the region above and 
below the x-axis—t.e., 

Si + So = Se (1 
on the positive x-axis. In addition to the requirement 
that the resultant pressure be continuous, it is neces- 
sary'! that its normal derivative vanish on the wall and 
be continuous along the positive x-axis. The solutions 
S; and S: to the simple wave equation in half space can 
be expressed in terms of their normal derivatives on 
the x-axis—i.e., S),(x > 0, 0, t)—and consequently 
Eq. (1) yields an integral equation for Sj, or S», along 
the positive x-axis. The solution to the integral 
equation is not presented in reference 11. 

A more general problem, the case of an arbitrary 
pulse incident on an infinite wedge of any angle, has 
been considered by Kay.'! Kay’s approach differs 
from that of reference 11 in the following respect 
The definition of the incident wave together with the 
reflections from the side of the wedge for t < 0-—i,, 
before crossing the corner located at the origin—is 
extended in the same expression only to the region 
outside the characteristic cone ct = Vx? + y® in the 
space-time continuum. The problem is then to find 
a solution of the wave equation inside the cone. For 
this solution the condition of continuity of the solution 
across the cone determines a certain nonhomogeneous 
boundary condition on the cone in addition to the homo- 
geneous boundary condition already imposed on the 
wedge. Hence for the special case of a wedge angle equal 
to 90°, Kay’s solution inside the cone should be identical 
with the solutions in reference 11—1.e., with Sy) + 9 
and S2 in the regions y > 0 and y < 0, respectively 
whereas both S; and S» fulfill a homogeneous boundary 
condition on the characteristic cone. Kay has ob- 
tained, in principle, the solution to his problem. How- 
ever, the intricacy of his method precludes the poss: 
bility of obtaining the final results of any problem othet 
than the conical problem treated by Keller and Blank.’ 
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DIFFRACTION 


Turner” proposed to determine the Green’s function 
for the wave equation for the case of the wedge by 
“two transformations and two inversions.” The pro- 
cedure, in principle, is straightforward. Neverthe- 
less, it is only illustrated by the problem of diffraction 
of a cylindrical pulse by a half plane. 

Graphs of the theoretical pressure distribution'’ re- 
sulting from the diffraction of a weak shock by a rigid 
right-angled corner are compared with graphs of meas- 
ured pressure distribution by Keller.!6 The 


ment is found to be satisfactory except at the corner 


agree- 


of the wedge and at the ‘‘triple-point’’ where the inci- 
dent (or reflected) pulse front is tangent to the char- 
acteristic cone. 

At the corner of the wedge, the velocity components 
based upon the linearized theory or acoustic theory 
become infinite although the pressure remains finite. 
The experimental result shows the existence of a small 
region of finite vorticity near the corner. Theoretical 
investigations of the relevant vortex generation have 
been made by Howard and Matthews’? and Rott.’ 
In the vicinity of the “‘triple-point’’ Lighthill's tech- 
A semiempirical approach has been 


nique! °° fails. 
advanced recently by Sun.*! 

In the present paper, the problem of diffraction of dis- 
turbance or an arbitrary pulse by a convex right corner 
is again considered. The classical acoustic approxi- 
mation for weak pulses is used, while the generation 
of vortices is ignored and improvement of the acoustic 
solution near the characteristic cone is not included. 
The main purpose is to present an elementary and in- 
tuitive method which permits final numerical results 
to be obtained conveniently, rather than to derive a 
solution purely in principle. The approach is based 
on the same physical ground outlined in reference 11. 
Instead of an integral equation for the normal deriv- 
ative of pressure along the positive x-axis, an integral 
equation for the pressure distribution along the x-axis 
issetup. The necessity for doing this will be explained 
A solution to this new integral equation ‘is pre- 

The final results can be obtained by simple 


later. 
sented. 
quadratures. 

Three fundamental types of incident pulses are con- 
sidered and the response obtained in each case. Two of 
them do not possess the properties of a conical field. 

By superposition of these fundamental solutions, 
many diffraction problems and wing-body interference 
problems can be solved. illustrated in the 
present paper by the treatment of the diffraction of a 


This is 
weak shock by a rectangular structure. Likewise, the 
corresponding problem in supersonic flow (the inter- 
ference problem of a prismatic body of rectangular 
cross section with a planar wing having supersonic 
In addition a method of solution 


edges) is treated. 


for wings with subsonic edges is presented. This type 
of wing-body interference problem was first considered 
by Lagerstrom and Van Dyke.”? Discussions of addi- 
tional applications of the results to be presented here 


can be found in references 23 and 24. 








DISTURBANCES 





(1) FORMULATION OF THE DIFFRACTION PROBLEM 


In Fig. 1, the 90 
the third quadrant of the x-y plane. 


convex corner is represented by 
The given two- 
dimensional pressure disturbance, So(x, y, t), which is 
a solution of the simple wave equation, is defined in 
where C is 
fulfills 


the region y > 0 and vanishes when x > Ct, 
It is understood that S 
the boundary condition (0/Oy)So(x, 0*, 1) = 0. 


the speed of sound. 


When ¢ > 0, the pressure disturbance S)(x, y, £) 


creates a discontinuity in pressure, AD = Sp(x, 07, ¢), 
along the portion of the x-t plane where 0 < 4 ce.” 


The additional disturbance, which will be confined in 
side the cone, 
F(x, y,t) = Ct-—- Vx? + y? =0 

should be such that the total pressure disturbance, p(x, 
y, t), and also its first derivatives should be continuous"! 
across the portion of the x-/ plane where 0 < x < Cr 
Let Si(x, y, t) and S2(x, y, t) be the additional pressure 
disturbances in the region y > 0 and y < 0, respectively. 
The problem now is to find S; and S». Both of them 
are solutions to the simple wave equation 

(0*p/Ox?) + (07p/dOy?) — (1/C*) (0?p/dt?) = 0 (2) 


and fulfill the following initial conditions: 


Si(x, y,£ << 0) = 0 
(3) 
So(x, y, ?< 0) = 0 
and boundary conditions: 
Sox > 0,0-,¢ >0) — Si(x > 0, 0+,¢>0) = 
ols > 0,07, > 9 (4) 


(0/Oy)Sox(x > 0,0-,¢ >0) = 
(0/doy) Si(x > 0,0*,t>0) (5) 


(0/Oy)Si(x < 0,0*,¢>0) = 0 (6) 
(0/Ox)S(0*, y << 0, >0) = 0 (7) 
S; = So = 0 on and ovtside the cone F = 0 (8) 


(2) THE INTEGRAL EQUATION AND ITS SOLUTION 


It is clear that Eq. (7) can be replaced by the condi- 
tion of symmetry, 


So(x, y < 0,¢ > 0) = S(—x, y < 0, t > O) (9 


Due to the boundary conditions (5), (6), and (9), the 
solutions of the wave equation in half space, S; and S, 
can be expressed in terms of the following integrals: 


: l G(&, r)didr 
5 == 7 
7 Vit— 7)? — [we — €)? + ¥)//C 
(10 
ey G(\é|, 7) dédr 
So = 11) 
7 V(t — 7)? — [(x — £2? + v*|//C? 


* The case in which the discontinuity is created in the super 
sonic region is treated at the end of Section (5). 











824 JOURNAL OF THE AERONAUTICAL SCIENCES NOVEMBER, 1957 
where G(x < 0, t) = 0 for all values of 8 > a > 0. 
Although this equation appears to be an equation 9, 
My = § 4 t = a if - . ee ° e oa . 
Ge > 8, 8) Sux > 0, 0%, f) Su(x > 0,0", #) Fredholm’s type,” it is actually quite different becaug 


Elimination of the unknown function G(x, f) from 
(10) and (11) will yield an integral relationship 
between S». and S;. With this relationship, S2 will be 


eliminated from Eq. (4) and an integral equation for the 


Eqs. 


pressure disturbance along the positive x-axis is ob- 
This approach is carried out as follows. 
method* * the characteristic 


tained. 


Following Evvard's 


coordinates are introduced, with 


Ct — x, Ct+x 


Eqs. (10) and (11) then become 
—|] in da’ pee 
. / , / , 
27 Jo Va—a'da VB—B 
l a. da’ *8 g(a’, B’)dB’ 
hoa, 8) = 9 | / ’ | / , 
-7TJ0 Va-— ada VB-B 


: ¢ da’ «’ 9(B’, a’)dp’ 
9 } / F = ( 13a) 
-7TJ0 Va-—a’'J0 VB-B 


= —hy(a, 8) + 


l [" da’ [ g(a’, B’)dp’ 
2a J 0 VB —_ a’ Ja’ Wa me B’ 

hi(a, B) = Sails. OF, t) 

hoa, B) = So(x, Ut) 


hi(a, 8B) = (12) 


ho( a, B) 
(13) 
where 


g(a, B) = G(x, ft) 


(13) holds for B > a > O (or 
(13) is actually 


Both Eq. (12) and Eq. 
x > 0) only. The second term in Eq. 
equivalent to —/,(8, a) [or —.S\(—-x, 0+, t)], therefore, 


hola, B) = —hy (a, B) — h(B, a) (14) 


and Eq. (4) becomes 


Phila, B) + hy(8, a) = —hij(a, B) (15) 


with fo(a, B) = So(x > 0, OF, tf) and B >a. h(8, a) 
can be expressed in terms of /1(a, 8) as follows:® 


l da’ ra g(a’, B’) dp’ 
hy(B, a) = 


7 ah VB — a’ Wa — ~’ 


--f{ da’ 
VB = a’ 


Oo @ hy(a", a)da” 


/ 
Oa’ Jo V a! — a’ 


V/B ta a hy(a’, a)da’ 
T Jo (8 — a’) va — a’ 


(15) becomes an integral equation for 


(16) 


and then Eq. 


hy(a, B), 
1 VB —* © hi(a’, ada’ 
T 0 (B — a’) Vae-e 
—ho(a, B) (17) 


2hy( a, B) 


y-axis has the effect of a mirror. The 


of the replacement of 8 by a in the function /), and th, 

appearance of the factor (8 — a’) in the integrand. 
Nevertheless the integral Eq. (17) is solved by th, 

successive substitution, as is the Fred 


for B > a, 


method of 
holm’s type—i.e., 


hy(a, B) = h(a, B) + h(a, B) + h@ (a, B) 


(18 
where 
h(a, B) [ho(a, B) 2] 
/ / 
ve ot" h (a’, ada’ (19 
h(a, B) = - 
2a J0 (B — a’ )Va-—a 
(n 1) , 
' ‘i A =o ra fj (a’, a) da 
h'’(a, B) = é (20 
2 (8B —a)Va-a 


<7 JV 


For most of the physical problem, /:"’(a, 8) is always 


bounded,* therefore, there exists a constant .1, so that 
h , (a, B) ~ A 
for all valuesof 8 >a. It follows that 
V3 = mae Ada’ 
h®(a, B)| < : ; = 
2a Jo (8B — a’)Va— a’ 


| ? 

: J 
tan~! \ < 

T  — 2 


a 


and siinilarly if 


it follows that 


h(a, B)| < (max/2)|h"~?(a, B)| < A/2’ 


Hence, the solution by successive substitution, Eq 
(18), converges faster than the power series of 1/2. 

The method of successive substitution for this prob- 
lem can be interpreted on physical ground. For the 
zero order approximation, the diffraction problem by 
the convex 


quadrant is approximated by the corresponding prob- 


right corner represented by the third 


lem with a flat plate along the negative x-axis. Due 
to the condition of antisymmetry, it is obvious that the 
pressure disturbances in the regions of y > 0 and y <0 
in pressure along the 
respectively, h(a, 8) 
(19) for 6 > a. Consequently 
(20) represents the corre- 


which balance the discontinuity 
positive x-axis are, and —h 
as defined by Eq. 
—h (8, a) defined by Eq. 
sponding pressure disturbance on the lower surface 
of the “‘imaginary”’ flat plate. 

The presence of the solid surface along the negativ 
image of the 
approximate solution in the third quadrant is added t 
the corresponding solution in the fourth quadrant. As 


* This is valid if Ao(a, B) 
if Ao possesses a logarithmic singularity or a singularity of a7! 


asa~>vU. 


or (So) is bounded and still true eve! | 
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DIFFRACTION 


i result h™ (a, B) 


is created along the positive x-axis and one cycle of the 


. new discontinuity in pressure Ap = 


successive substitution is completed. 

In reference 11, an integral equation for the normal 
derivative of pressure along the positive x-axis was ob- 
tained by substitution of Eqs. (12) and (13a) into Eq. 


t), 


"a g 3B". a : "6 
dp’ + 2 | 
0 \ » ee B’ Ja 


B’) 


g(a, 
is = 
Vv B _— B’ 


O [ h(a’, 8) 
Oa, 0 Va— a’ 


Omission of the first term will change the equation 
into the Abel's the 
method of substitution can be applied and has actually 
(1S). 


da’ (21 


one of type.” Consequently, 
the same physical meaning as the one for Eq. 
a proof of convergence of the solution for 
was lacking in reference 11 due to the diffi- 
the property of the 
field. In addition to the compli- 
operation in double 


However, 
Eq. (21) 
culty of 


assessing derivative of 


in the flow 
cacy of the 
integration is required to obtain the pressure distribu- 
tion from the solution of (21), 
tion of Eq. (17), which is the pressure distribution, can 
be obtained by simple quadrature. the new 
formulation, Eqs. (17)-(19), is well suited for obtaining 


pressure 
each substitution, a 


Eq. whereas the solu- 


Hence, 


numerical results. 


The convergence of the series of Eq. (18) will be 
faster if it is regrouped as follows: 
hyla, re = (2 3)h ” (Q@, B) + 
(1/3) \h (a, B) + 2h(a | + 
[h'? (a, B) + 2h a, B)| + = 
h(a, B) + (1/3) > [h\” (a, B) + 
1, 2,3 
2h'”’ (a, B)] (22) 


In the series of Eq. (18) the zero order approximation 
for Iya, B)—-.e., 


h(a, B) = —(1/2)ho(a, B) 


assuming 
22) 


(17), (15), by 


a; while in the series of Eq. ( 


is obtained from Eq. or Eq. 
0 for B > 


the zero order approximation for /;(a, 8)—4.e., 


hy(,, a)= 


(2/3)h (a, B) = —(1/3)ho(a, B) 
is obtained by assuming 
hy(a, B) hy(B, a) 


The latter assumption fulfills the condition that the 


value of the pressure at a = 6 is unique—+.e., 


lim [/,(a, 8) — M(B, a)| = 0 


B-a-—>0O0 


The former assumption leads to a discontinuity of 
(1/2)ho(8, 8) at a = 8. Hence, it is expected that the 
series of Eq. (22) converges faster than the series of 
Eq. (19). 
(22) is that the solution corresponding to the Nth sub- 
stitution—i.e 


Another reason in favor of the series of Eq. 


OF 


DISTURBANCES 825 
\ 
HT, RY oe (2/3)\A@ + (1/3 a \h 7 4 2h 


yields the correct behavior for h; at both end points, 


a = 0, and @ = g.* On the other hand, the series of 
Eq. (1S) converges to the correct value of —(1/3)/(3, 
8) at the end point where a = @ in the same manner 
as the series (1/2) — (1/4) + (1/8 does to 1/3. 


In case /4y(B, 8) = O for all values of 8 > 0, then there 
is no longer any reason to prefer the series of Eq. 
to the series of Eq. (18). 


After hy(a@ 3), or 1, ’ \a@, 


(22) 


8), is determined for 8 > a, 


the pressure distribution can be obtained from the 
following formulas: 
p(x, vy >0,t) = So + S; = So - 


1 oO £f5S.(k, 0, r)dédr 
7 OY {J R 


; *So(&, O, r)dédr . 
j= So = a (25) 
T OV R 


p(x, y < 0, 


where 
Si(x > 0, 0, t) = h(a, B (26) 
Si(x < 0, 0, t) = —So(—x, 0, t) — 2.S1(—x, 0, ¢ (27) 
Sex > 0 ,t) = $.(—x, 0-, 2) = 
So(x, OF, t) + Si(x, OF, ¢4) (28) 


FUNDAMENTAL SOLUTIONS 


(3) 


the 
So is a conical solution with apex located at a 
—L y = 0. 
So, in the first 


In many practical problems, incoming disturb- 
ance 
point other than the origin, say at x = 
The discontinuity in pressure, created by 


quadrant of the x-f plane, in general is expressed by a 


series, 
So(x, OF, t) = Ao + A1V 1 — [(x + 1)/(Ct + 1)] 4+ 
As\1 — [(x + 1)/(Ct + 1)]}*/?7 + (29) 
where Ao, Aj, 43, . are constants if So is bounded.t 
The sum of the first three terms usually differs less 
than 5 per cent from the value of Sy for 1 >t>x > 0 
[see Sections (4) and (5)]. It is, therefore, worthwhile 
to find the additional disturbances S;(x, y > 0, t) and 
So(x, y < 0, ¢) corresponding to 
: oe : 
So(x, OF, ¢) = 1, V (Ct — x)/(1 + Ct) 
* It can be shown that 
lim HH,“ (a, 8) = h® (0, 8) = —(1/2)h0(0, B 
a—) 
lim = Hy) (a, 8) = (2/3)h (8, 8) = —(1/3)ho(8, 8 
B-a—0 


+ Otherwise, the terms representing the singularities of Sp will 
appear in the series of Eq. (29) 
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Fic. 2. Pressure distribution corresponding to S)(x, 0, 4) = 1 
. +,\ 13/2 P , » 
and [(Ct — x)/(1 + Ct)]"’”, respectively. S$) and S» 
are called the fundamental solutions. 
(i) SJl(x > 0, 0+, 4) = 1 
he discontinuity in pressure So (x, O*, ¢) = 1 1s 


created when a plane weak shock (or sonic pulse), x = 
Ct, is passing over the 90° convex corner (Fig. 2). 
The pressure distribution can be obtained as a conical 


10 


solution® 


ee l (1 — 2)V3 ' 
i eee = 1 — -tan™! P (30) 
ta T (1 + Q*) — 42 cos w 


where tan~! is taken between 0 to z, 
Q= {1 - V1 — 7?) rs 
w = 20/3 

r= Vx? + y?/Ct 

oe = tan! (y/x) + (9/2) 


Therefore, the fundamental solution corresponding to 

Ry = ] is 

Si'(x, y > 0, t) = pe — So = p(x, y > 0, #) — 1 (31) 
S(x,y < 0, t) = px, y < 0, 2) (32) 


This problem is now re-solved by the method pre- 
sented in this paper in order to learn how fast the series 
of Eq. (23) converges. With So(x, 0*, t) = 1, Eqs. (19) 
and (20) yield 


h® = —(1/2) 
h(a/B) = (1/27) tan IW (a B)/{1 — (a@/B) | 
h®(a/B) = (—1/278) V a(8 —a)X (33) 
v1 h'? (da 
j (1 — A(a/8) V1 — A 


The last integral can be evaluated by numerical integra- 
tion. Fig. 2 shows the pressure distribution along 
y = 0, given by Eqs. (24), (26) and (27) where /y(a, 8) 
is approximated by the results of first and second sub- 
stitutions, J,“ and f°, respectively. It is clear 
that the result based upon the second substitution is 


NOVEMBER, 1957 


already very close to that of the exact linearized conical 


solution—i.e., Eq. (30). 


(ii) So'(x > 0, O+, 1) = V(Ct — x)/(Ct+ 1 


In order to be able to carry out the first and the 
second substitution without resorting to numerical 
; : II \ . 
integration, So (x > O, 0*, ¢) is expressed as a power 


series of (8 — a)/(2 + 28), 
Sol'(x > 0, OF, t) = Iola, 8) = 
V 2a/(2 + B+ a) = V2/[2(1 + B) — (8 — a)] = 
Va/(1 + 8) $1+ [(8 — a)/4(1 + B)] 4 
(3/32) [((@ — a)/A1+ 6)]? +... 5 (34) 


The error in using only the first three terms is 
less than 3 per cent for 2 > 6 > a > O. Cor- 
responding to these three terms, the integrals in Eg. 
(20) can be evaluated.*° The solutions for the first 
and second substitutions differ less than 2 per cent 
from each other. Because of Eq. (20), the series of 
Eq. (22) is an alternating series. The exact solution 
lies in between two successive substitutions. Hence, 
the result from the second approximation is satisfactory 
and there is no necessity to carry out the third sub- 
stitution. 

Since the denominator in the integrand in Eq. (20 
vanishes as a’ — a, the value of h(a, 8) depends pre- 
dominantly on the value of h(a’, a) asa’ >a. Itis 
clear that Sy'' approaches the first term on the right 
side of Eq. (34) asa—~ 6. When this first term is used 
for the evaluation of h“, the solution from the first sub- 
stitution is 
Hi (a, B) = h© + (2/3)h® = 

(—1/2)V 2a/(2 + a + 8) + [1/(6V1 4+ a)] Xx 


[V8 — VB— a] (35) 


The maximum difference between the simple solution 
of Eq. (35) and that from the second substitution cor- 
34), 


responding to the first three terms in series Eq. 
is 2.5 per cent. 

The additional pressure disturbance S,''(x, y > 0, f 
and S"'(x, y < 0, t) can be obtained from Eqs. (24) to 
(28) where /,(a, 8) is approximated by //,") in Eq. 
(35).°° The error due to this approximation is less 
than 2.5 per cent. 

(iii) SoM(x > 0, OF, t) = [(Ct — x)/(Ct + 1)]””" 

The additional disturbance pressure S,'""(«, y > 0, f) 
and S.""(x, y < 0, t) due to S,'"' can be obtained in the 
same manner as S,'' and 5S." (see reference 26). 


(4) DIFFRACTION OF WEAK SHOCK BY A RECTANGULAR 
BARRIER 

In Fig. 3 is shown a rectangular barrier of height L 

and width 2L. (The generality of the problem is not 

lost by specifying the ratio between the height and the 
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Fic. 3. Pressure distribution on the barrier 


The weak shock of strength ¢ hits the barrier 


During the interval 0 < t < L/C, the pres- 
» il 


width). 
at i - U0. 
sure distribution is given by the conical solution.'® 


(1 €)p l (x, y, t) =%9— pi, y, t) 


where p(x, v, ¢) is defined in Eq. (30). During the 
2L/C, an additional disturbance 
(1/e)p'?(x, y, t) = Pl—x, 
when p is reflected from the 90 
When ¢ > 2L/C, the pressure disturbance, p", is dif- 
fracted by the convex right corner at C and the dis- 
The continuations 


interval L/C < t < 
—y — 2L, t) is created 


concave corner at 4. 


turbance, Pp’, by the corner at B. 
of the solutions p" and p® for t > 2L/C create dis- 
continuities in pressure along the portion of the x-axis 
where x > 2L and along the positive y-axis, respec- 
tively. The discontinuities are expressed by the gen- 
eral form of Eq. (29)—e., 


(1/e)p(x > 2L, 0, t) = 1 + 0.17327 X 
V (Ch — x1)/(2L + Ch) + 0.0551 X 
[(Cty — Fil (OF 4. Ct) }* 2 4 
(1/e)p?(0, ¥ > 0,8) = 
—0.519798 V (Ch — y)/(2L + Ch) - 
0.06257 [(Ch — y)/(2L + Ch]? + 
where 


\ x—2L>0 and ¢, — 2L/C>0 


fhe additional pressure disturbances around the 90 
convex corners B and C can be obtained by super- 
position of the fundamental solutions in Section (3). 
Che pressure distribution on the barrier at the instants 
Ct/L 0:5, 1, 15, 2 235. ! is shown in 
After the instant ¢ = 
a positive integer, additional disturbances are 
created due to the reflection from the ground. After 
the instant ¢ = 2nL/C, the disturbances created pre- 
viously are diffracted by the right corners at A and C 


3, 3.0, and 
(2n + 1)L/C, where 


Fig. 3. 


nm 1S 


and the response for each case can be obtained by the 


method presented in Section (2). 


(5) INTERFERENCE OF A PRISMATIC BODY OF 
RECTANGULAR CROSS SECTION WITH A PLANAR WING 
IN SUPERSONIC FLIGHT 


The problem of wing-body interference at super 
sonic speeds has been subjected to intensive investi 
gation in recent years. A comprehensive study and 
survey was made by Lagerstrom and Van Dyke*® by 
Many 


results have been obtained for slender or conical wing 


Lawrence and Flax” and recently by Ferrari.”* 


body configurations. For the flow field about a wing 
body combination employing planar wings with super 


sonic edges and bodies deviating only slightly in shape 
from a right circular cylinder, theoretical methods are 
presented by Ferrari and Cramer,* by Nielsen,” and 
by Morikawa.*® The direct correlation between the 
first two methods is shown in reference 28 and the 
similarity between the third and the second method is 
pointed out in reference 27. These methods are feasible 
for practical computations due to the availability of 
the tabulation of the ‘“‘characteristic functions.’’** *! 

Besides the representation of fuselage by circular 
cylinders, another representation by prismatic bodies 
of rectangular cross sections was discussed by Lager- 
strom and Van Dyke.*? In this Section, the problem 
of the interference of a prismatic body of rectangular 
cross section with a planar wing (with supersonic or sub 
sonic edges) is reduced to two basic problems—the 
diffractions of disturbances by the subsonic edges of 
the wing and by the right corners of the body. They 
are solved, respectively, by Evvard’s method and the 
method presented in Section (2). 

This solution of the wing-body problem can be ap 
plied to study (i) the effect of “slight contouring” of 
the prismatic body, that is, the favorable shape of the 


body which differs slightly from a basic prismatic body 
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Fic. 4. Pressure distribution on body and wing [height/width = 1 


of rectangular cross section;** (11) the favorable inter- 
ference of wing-body configurations among the pris- 
matic bodies of rectangular cross section with the same 
volume (or cross-section area) but different height- 
width ratio; and (iii) the resultant effects of ‘‘slight con- 
touring,” (i), and variation of height-width ratio, (11). 

On the other hand, the solution of the interference 
problem for circular cylindrical bodies*~*! can be 
applied to study the effect of “‘slight contouring,”’ 
(i), but not the corresponding nonlinear effects of (ii) 


and (iti). 


(a) Wing-Body Interference—Wings With Supersonic 

Edges 

For wings with supersonic edges a basic configuration 
is shown in Fig. 4. An infinitely long prismatic body 
of rectangular cross section is mounted on a wing repre- 
sented by a semi-infinite wedge (or flat plate). The 
body is parallel to the free-stream direction and the 
upper surface of the wedge is inclined at a small angle 
5. The origin is located at the point where the plane 
of the weak shock at the leading edge (y = h, z = 
—Bh) hits the 90° convex corner B.* The discon- 
tinuity in pressure across the weak shock is denoted by 
e which is equal to J/*pC*6/B. When the shock is 
passing over the 90° corner at B—i.e., > 0—an addi- 
tional disturbance p(x, y, z) is created, 


p(x, y, 3) = ep(—y, x, 8) 


where p(x, y, 2) is the conical solution defined in Eq. 
(30) and Ct is replaced by 2/B. For z > hB, p is 
reflected from the surface of the wing (y = L). The 
additional disturbance due to reflection, which is 
p(x, y, 3) = p(x, —y + 2L,2) — « 
passes over the right corner B for s > 2Bh. The con- 
tinuation of p® into the region zs > 2Bh creates a dis- 
continuity in pressure across the y-z plane—.e., 


* It is sufficient to discuss only what happens on the right side 
of the plane of symmetry (x = —Z). 


Aply < 0,21 > 0) = p@(0, y < 0, s 
= —0.1734 e[(z, + By)/(a + 
2Bh) |! ae 0.55 e| (2 +- By 


> 


(s; + 2Bh)]’’" 4 
where 2} = z — 2Bh. The additional disturbances 
which balance the discontinuity are obtained by super- 
position of the second and the third fundamental solu- 
tions in Section (3). Similarly, the pressure disturb- 
ance p‘” is reflected from the plane of symmetry 
(x = —L) tor zs > BL. The reflected disturbance, 


p(x, y, 3) = p(—x — 2L, y, 2) 


is again diffracted by the right corner B for z > 2BL 
and the discontinuity in pressure across the x-z plane, 


which is 


Ap(x, 2) = p(x > 0, 0, 21) = 
0.51979Se V (sg; — xB)/(2LB + 2) + 
0.06257€[(s) — xB)/(2LB + 2) |° at 
for sz, = s — 2BL > 0, is removed by the superposition 
of the fundamental solutions. 

As in the problem of diffraction by a rectangular 
barrier, after the section of x = (mL + nh)B, where n 
and m are positive integers, additional disturbances 
will be produced due to the reflection of disturbances 
from the wing surface or the plane of symmetry or due 
to diffraction of disturbances by the right corner at 
B. The basic diffraction problem is treated in Section 
(3). 

Figs. 4, 4(a), 4(b), and 4(c) show the pressure dis- 
tribution at various sections on the surface of the body 
and the upper surface of the wing for width-height 
ratios equal to 1, 0.5, 2, and 1.5, respectively. The 
integrals of the pressure distribution fulfill the integral 


relationship,”*: **1.e., 


a0 
| p(x, 0, 2) dx { 
—T] 0 


Thus it serves as a check to the numerical results. 


[e — p(x, y = h,z+ BL)| dx 
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(b) Wing-Body Interference—Wings With Subsonic Edges 


For wings with subsonic edges, the pressure dis- 
tributions can be obtained by repeatedly using Ev- 
vard’s method* * and the method presented in Section 
(2) with or without modification. The required modi- 
fication can be illustrated by an example as shown in 
Fig. 5. The disturbance that originates at point B on 
the 90° convex corner of the body is reflected from the 
wing surface and passes over the subsonic leading 
edge at point D, the origin of the x, y, s coordinates. 
The incoming disturbance po(x, +y, 2) + pox, 
— y, Z) creates the discontinuity in pressure, 2/9(«, 0, z) in 
the x-z plane between the leading edge, 4 DG, and the 
hyperbola, CDE. The additional pressure disturbance 
can be split into two parts—namely, 


pi(x, y, 3) = pr“ (x, y, 3) + pi (x, y, 2) 


The discontinuity in the supersonic region—i.e., be- 
tween the sonic line DF and the hyperbola—is removed 
by pi‘, where 


ra 


pi (x, y > 0,2) = —~A (x, y < 0,2) = 


l if h(é, ¢), déd¢ 
t JJ V(s — £)? — [(« — &)? + y? |B 


h(x, s) = (0/Oy)po(x, OF, s) for x>2/B 
h(x,z) =0 for x«<2/B 


The discontinuity in the subsonic region—-that is, be- 
tween the leading edge and the sonic line, is now 2o(x, 
OF, 2) + 2p:(x, O*, 2). This discontinuity is re- 
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moved by fi, which can be obtained by Evvard’s 





method. 

Similar modification is necessary when a disturbance 
originates at the leading edge (point D) and _ passes 
over the 90° convex corner of the body (point // 


CONCLUSION 


The problem of the diffraction of a two-dimensional 
disturbance or an arbitrary pulse by a convex right- 
angled corner is solved by a method which permits 
final numerical results to be obtained conveniently 
The solutions obtained are applied to the analysis oj 
the two-dimensional unsteady flow field associated 
with the diffraction of a weak shock by a rectangular 
barrier and to the analysis of the equivalent three. 
dimensional steady supersonic flow field produced by 
the interference of wings and prismatic bodies of 


rectangular cross section. 
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Transonic Flow With a Detached Bow Wave 
Past a Wedge Between Two Parallel 
Plane Walls 


SHIGEKI MORIOKA* 
Yoshida College, University of Kyoto 


SUMMARY 


In the present paper the two-dimensional transonic flow with 
, detached bow wave past a doubly symmetrical double-wedge 
section placed midway between two parallel walls is investigated 
The subsonic flow field over the front half of the profile is de 
termined by applying the relaxation technique to the transenic 
perturbation equation in the hodograph plane. The method of 
characteristics is adopted for the supersonic flow over the rear 
half. Detailed computation has been carried out for three 
values of the ratio of the chord length to the tunnel width, as- 
suming the value of the transonic parameter to be 0.921 in every 
Cant 

Of the results, the following may be noted. When the tunnel 
width is gradually decreased, the bow wave moves forward and 
the pressure coefticients over the profile increase uniformly. 
However, a definite limit for the reduction of the tunnel width 
to maintain steady flow of a given (supersonic) Mach Number 
seems to exist. It is not clear what happens exceeding this limit. 
[he distance between the bow wave and the profile is com- 


paratively small even in this limiting stage. 


SYMBOLS 


= critical velocity 
= airfoil chord 


( = pressure coefficient 

® = reduced pressure coefficient (4 + 1! x 
(t/¢ 2/6 Cp 

,h,k = length of irregular lattice intervals relative to 

that of basic interval 

H = width of tunnel 

i = integral defined by Eq. (23 

bs = (7 = 1, 2, 3) component integral [see Eq. (20 

x. i = functions of e, f, g, h, and n, [see Eqs. (15) and 
(16 

k = numerical constant [see Eq. (9 

WV = Mach Number 

WV = Mach Number of uniform flow in front of bow 
wave 

p = pressure at exit of tunnel 

S(n = function defined by Eq. (17) 


= airfoil thickness 
= local speed of flow 
) = Cartesian coordinates in physical plane 

a = absolute value of 7 at return point of shock polar 
= absolute value of 7 at left-hand limit of lattice 
= ratio of specific heats (1.4 for air 

A = basic lattice interval 
= hodograph variable defined by Eq. (2 
= value of n corresponding to uniform flow ve 

locity in front of bow wave 

= reduced hodograph variable defined by Eq. (1 


= local inclination of flow relative to x-axis 
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6 = reduced hodograph variable defined by Eq l 
6 = variable of integration [see Eq. (8 
by, = half-angle of wedge 
By = value of 6, corresponding to wedge half-angle 
g = speed function (1/7? l vy +1) (t/e 
g = transonic similarity parameter (.1/,? l 
(¥ 1) (t/c)]? 
px = fluid density at critical speed 
v = stream function 
Yu = value of ¥ on tunnel wall 


Subscript 


0, 1,2,ete. = value at a prescribed lattice point 


INTRODUCTION 


A THEORETICAL INVESTIGATION of the aerodynamic 
characteristics has been performed by W. G. 
Vincenti and C. B. Wagoner! at zero angle of attack 
of a thin, doubly symmetrical double-wedge profile in 
the range of supersonic flight speed in which the bow 
wave is detached. They have determined the sub 
sonic flow field between the bow wave and the profile 
by applying the relaxation technique to a boundary 
value problem for the transonic small disturbance equa 
tion in the hodograph plane. On the other hand, they 
have used the method of characteristics in calculating 
the purely supersonic flow over the rear half of the wedge 
profile. 

The flow field as considered by Vincenti and Wagoner 
can be realized in a closed wind tunnel, the width of 
which is sufficiently large compared with the model 
size. The reflected shocks from the tunnel walls, in 
this case, do not affect the flow over the model profile. 
Now we consider what happens when the width of the 
wind tunnel is gradually decreased. As far as the 
shock angle at the wall of the bow wave lies within the 
range | of Fig. 1, the bow wave will be reflected regu- 
larly from the wall. When the tunnel is more nar- 
rowed and the shock angle enters the region II, the 
reflection will become of Mach’s type. Further, if the 
stage is reached in which the reflection occurs at a 
subsonic shock angle, the bow wave becomes normal 
at the wall and the reflected wave will disappear. It is 
this stage to be considered in the present paper. 

The flow field under consideration is in general 
governed by four parameters:? (1) Mach Number 
My of the uniform flow in front of the bow wave, (2) 
the pressure p at the exit of the tunnel, (3) wedge angle 
26.,, and (4) ratio cH] of the model length to the tunnel 


width. The parameter (2) has little influence upon 
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the flow upstream of the exit by taking p to be a super- 
sonic pressure. The parameter (4) represents the so- 
called tunnel effect and it is closely related, as is seen 
later, to the location of the return point of the shock 
polar in the hodograph plane (the return point is the 
image of the inflection point of the shock wave in the 
physical plane). As the return point moves to the 
point of normal shock condition, the width of the tunnel 
decreases to a minimum and then increases until the 
(subsonic) choked condition is reached when the two 
points coincide with each other. 

The return point is a singularity in the hodograph 
transformation. 
ally becomes zero at this point. 
shown® that this singularity is isolated provided it lies 


The transformation Jacobian actu- 
However, it has been 


in subsonic region some distance apart from the sonic 
line. (The necessary distance is smaller as the shock is 
weaker). 
physical and the hodograph plane is one to one. 
the flow field in question forms a regular boundary value 
problem for the Tricomi equation in the hodograph 
However, it seems difficult to treat the prob- 


Therefore the correspondence between the 
Thus, 


plane. 
lem analytically because of the complicated boundary 
conditions. In the present paper, a numerical method 
similar to that used in reference | is adopted—the flow 
over the front half of the wedge is determined by the 
relaxation method, while the flow over the rear half is 
calculated by the method of characteristics. Actual 
computation has been carried out for three cases of the 
flow, each for different value of the ratio c/HZJ. The 
transonic parameter for each case is assumed to have 
a value 0.921. (This value is adopted for the sake of 
comparison of the results with those for free flight con- 
dition due to Vincenti and Wagoner). 


BOUNDARY VALUE PROBLEM 


We now consider the two-dimensional transonic flow 
with a detached bow wave past a thin double-wedge 
profile placed midway between two parallel plane 
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walls (Fig. 2). (In the following description only the 


flow in the upper half-plane will be described, since the 


flow is symmetrical about the chord line.) The 





hodograph representation of the flow field is shown in 


Fig. 3. Fig. 4 shows the reduced hodograph plane in 
the transonic approximation. In these Figures jt 


may be noted that the shock polar returns back at some 
point B (return point) in the subsonic region. The 
characteristic DF corresponds to the so-called limiting 


Mach wave. 


Using the reduced hodograph variables 
n= 7/n, 8 = 2"?(6/n9%!?) (] 
with n= (y + 1)"8 [((V — ax)/ax] 9 
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boundary value problem in the n-plane can be 


the 
summarized as follows. 


Differential Equation 


¥ 2nvia = O 


Boundary Conditions 


On the basic streamline AOE, 


y 0 for 7» Ss —l, 6=0 
y 0) for n < 0, A 4, 
y>0 for n>—~, 0S 6<8 


On the shock polar ABC, 


for -l ins —-a, 6=(1—y7)vV1- n 


lor n = U, 
where 

k= 248. ¢/3'/*8T3(1/3) = 0.341 
On the solid wall CD, 


y vu for —-lsns0, 6=0 


REDUCTION TO THE FINITE-DIFFERENCE EQUATION 


Phe formulation of the finite-difference equation can 
be accomplished by grouping conveniently the lattice 


points into six categories similar to reference 1. 


arrangement of the relative lattice points and the dif- 
ference equation to each point in the present problem 


are as follows.7 


3) Points Adjacent to the Shock Polar 


/ a ed 
xk 
J 
( POP Sood. Sook 
Rl +k) (2+ k)"° 1+K " 24+k" 





t For the details of derivation of the difference equations see reference 1 












Regular Internal Points 





Points Far to the Left 
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(4) Points on the Shock Polar 
i yf #4 7 
‘ Re 
s t 4 ~ 
Fic. 8 
(1+ k) (2+ k) = 23 + 2k)L R(2+k) + 4(1+)L k(1 + k) + 2(1 + 2k)L 
= :T Y: = eghy, t 
2k l+k 2(2 + k) 
— 2+ 6k + 3k? 6L ete 
thKy; — + (eg —fh)K |} %o = 0 (lI4 
k(1 + k) (2 + k) k(2 + k) 4 
K= j[gh(e + f)/S(m)] + ef(g + hf, L = (eh — fe) [(gh/4i0) — (ef /2) JA 
K = [gh(e + f) + ef(g + h)S(m)]-, L = (eh — fg) [(gh/2) — efnl|K 6 
S(m0) = [11 + 7m0)/(B + Sm)] V1 + ao 7 
Where upper and lower signs in Eq. (14) correspond 
to the developments to the positive and negative direc- ¥>(0, 6) = (1/A) [(7/6)bo — (4/3); + 1/6): 
tions of the axes respectively. Also A and LZ in Eqs. 18 
(15) and (16) correspond respectively to horizontal 
and vertical developments. together with the relations 
* ea) ft 
(5) Point on the Sonic Line - 6 — 6) d0’ = 3.533 (Yo/A?’*) (19 
J bw (0’ — 6)?/8 
The difference equation for these points may be ob- ; . , , 
tained by putting the following expression into Eq. (7): for the first point below the upper boundary, an¢ 
=J,+J/.+ Js (20 
Jy = [1/(nd)28] SS [4/(4 + 3m)] [(2/3)m/m! n”] vn, 
n 0 
(2/3)5 = 1,. (2/3)... = (2/3) (5/3)... Kom — 1/3] 
Jo = (3/478) SS [(m + 1)¥8 — m3] (Wm — Wms) 
m=1 
Js = (1/A*/*) [(2/7) eo — (87/28)ri + (33/14) Yo + (13/28) ¥-1] 


for general point below the first lattice point. 


(6) Return-Point of the Shock Polar 


Accounting for the symmetry of the shape of shock 
wave in the vicinity of its inflection point in the physical 
plane, the value of y at the return point may be taken 
as the average of y at two neighboring points which are 
facing each other on both sides of the shock polar. 


SOLUTION IN HODOGRAPH PLANE 


As is described in the preceding section, a finite- 
difference equation can be set up for each lattice point 
in the hodograph plane. Thus, our boundary value 
problem is replaced by a large number of simultaneous 


algebraic equations involving an equal number of un- 


(See Fig. 9.) 


known values of y. Since the number of unknowns in 
each equation is small, they can be solved conveniently 
by the relaxation technique. The mechanics of the 
relaxation process have been well described by various 
authors‘ and need not be gone into here. 

Actual calculation 


cases, each having different position of the return 


has been carried out for three 
point—i.e., for three values of the ratio cH of the model 
length to the tunnel width in the physical plane. The 
value of the transonic parameter is assumed to be 4 
(This value is selected 


Figs. 5 to 7 show 


fixed value 0.921 for all cases. 
out of those treated in reference 1). 
calculated patterns of the streamlines in the reduced 
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hodograph plane, together with the distribution of the 


lattice points used in the relaxation calculation. 


See Figs. 10-12.) 


TRANSFORMATION TO PHYSICAL PLANE 


The relation between the physical x,y-plane and the 
reduced hodograph 76-plane is given by the following 


pair of equations (see reference 1): 


d(x/c) = (1/41) 2avzdn + 76) (21) 
d y/Cc) = (1 tT.) (26, (¥ + 1) (t/c)}} dy (22) 
with 
70 
Tip = pudyc/23/2(y + 1)'/3n9!/? = | nwa (n, 9.)d7 
(23) 
By integrating these equations numerically from — ¢ 


to the point considered, we can find the corresponding 
physical coordinates. 

Fig. 13 shows the calculated shapes of the bow wave 
and the sonic line relative to the wedge profile and the 
tunnel wall for three cases of flow for specific values of 
the Mach Number 1.164, y = 1. 
ratio t/c = 0.1. 


centi-Wagoner’s and Marschner’s results’ for compari- 


! and the thickness 
In this Figure are also shown Vin- 
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Fic. 1 

profile and tunnel wall 
son. The Vincenti-Wagoner case is that of the free 
flight condition and Marschner’s is the (subsonic) 


choked state which may be regarded as the limiting 
case with the bow wave at infinity upstream. 
The speed function & in the subsonic region can be 


found directly from the following relation: 


g= (M 1)/[(vy + 1) (t/c)]?8 = 


ton = (21/3/6,,2/3)n (24) 


The speed function in the supersonic region can be 
obtained, after determining the subsonic region, by 
the method of characteristics (specialized to the tran- 
sonic small disturbance theory). Then, the pressure 
coefficient C, can easily be calculated by the relation 


(E — &) (25) 


é. = —2(¥ + 1) (1/3 (t c)?/8 


Fig. 14 represents the distribution of the reduced 


pressure coefficient 
C, = (7 + 1)™* @/e)-@"C, 


along the front and the rear surfaces of the wedge pro- 
file. Here the corresponding results for free flight and 
choked conditions are also shown for the sake of com- 
parison. The total pressure drag can of course be 
found by integrating over the wedge surface the pres- 


sure coefficient just given. 
DISCUSSION OF RESULTS 


From the results obtained in the preceding section, 
we can understand the wall influence upon a transonic 
flight condition with a detached bow wave. Figs. 15 
and 14 show that the bow wave moves forward and the 
pressure coefficients over the model surface increase 
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uniformly as the wall approaches the model. Hoy 

ever, it is to be noted that there is a definite limit for the J 
approach of the wall (unexpected result to the writer 

In this limiting stage, the mean mass flow at the nar. | 
rowest passage attains its maximum which is clearly J 
larger than that for the (subsonic) choked conditio, 
(Marschner’s case). Thus, the flow comes to a kind oj 
choked condition (Supersonic choked condition, say), J 
It is not clear what happens beyond this limit. It js § 
probable that the steady flow of the type considere 
above no longer exists, and the bow wave may suddenly 
run out of the tunnel entrance. (A similar phenome 
non is familiar in a supersonic diffuser.) 





It may bh 
noticed here that, as is seen from Fig. 15, there are tw 
possible solutions for a range of wall positions betweer 
the lower limit of approach and Marschner’s choked f 
condition. One solution belongs to the branch _ just 
traced above and the other to the branch that bridges 
over the lower limit and Marschner’s case. It is clear 
however, that the latter branch may have a chance t 
be pursued by the continuous transition from the fre¢ 
flight condition only when the wall is moved back fron 
its lowest position. Moreover, this solution is likely t 
be unstable—as is seen from Fig. 13, a larger distance 
of the bow wave requires a wider tunnel. Now as 
sume that in the presumably unstable region one of the 
flow patterns belonging to the second branch has been 
established and that the distance of the walls is then 
reduced by asmall amount. This will reduce the mass 
flow through the narrowest cross section and will cause 
an accumulation of air in front of the body. Conse 
quently the bow wave will be pushed upstream. But 
then a flow field would arise which requires a wider 
cross section of the tunnel—i.e., the shock will be un- 
able to find a stable position. 
rence of two solutions for one given boundary condi- 
tion, one may remember the well known behavior of the 
attached shock to an infinite wedge when the wedge 
In this case, of two 
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angle is increased monotonically. 
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possible branches of the solution, the one that is re- 
lated to the weaker shock is pursued under the usual 
conditions up to the detachment angle. 
Thus, it seems natural to the writer that in the present 
case also the first branch above will be 
followed up to (or at least to the proximity of) its limit- 
ing point (lowest wall position) if the tunnel wall ap- 
Better understand- 


so-called 


mentioned 


proaches the model monotonically. 
ing of the occurrence of the two solutions and of the 
supersonic choked condition may be obtained by the 
consideration of the following example. Assume that 
the nose angle of the wedge is such that in a purely 
supersonic flow with an attached shock the velocity 
after the shock reflected from the wall is exactly sonic. 
Furthermore let the wedge side have such a length 
that the reflected shock hits exactly at the shoulder 
(Fig. 15). This flow pattern has certainly a maximum 
of the mass flow in the narrowest cross section between 
the model and the wall—i.e., for this flow pattern the 
distance of the walls has its absolute minimum. 

For a greater distance of the walls there may exist 
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two solutions—-a supersonic one and the solution with 
a detached shock. 
“lambda” shocks for wall distances slightly larger than 


(The latter may involve attached 
the absolute minimum.) The two branches correspond 
clearly to the first and second branches considered 
above. The flow pattern in Fig. 15 also corresponds to 
the supersonic choked condition mentioned before. 
Further, this example shows clearly that the supersonic 
choked condition can carry a greater mass flow than the 
subsonic choked condition. 

Fig. 14 shows that the pressure distribution over the 
profile has in general close resemblance to that for the 
free flight condition with the same value of the transonic 
parameter. It increases considerably, however, as the 
lowest wall position is approached and then tends to 
the distribution for the subsonic choked condition. It 
may be interesting to note that the distance between 
the bow wave and the profile also increases on approach- 
ing the limit, but the increment remains comparatively 
small. 

Finally, it may be remarked that, although in the 
present paper only the case where the width of the tun- 
nel is varied for a fixed Mach Number is considered still, 
it is possible to make a qualitatively similar argument 
for the case in which the Mach Number of the flow is 
varied in a given wind tunnel. 

This problem has been suggested to the writer by 
Prof. Z. Hasimoto, to whom the writer wishes to ex- 
press his hearty thanks. The writer's thanks are also 
due to Prof. K. Tamada for his advice and kind in- 
spection of the manuscript, and to Prof. M. Tada for 
his valuable discussions. 
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SUMMARY 


shc ck 
Various ar 


turbulent jet and a plane 


interaction of 
a shock tube was optically observed 


The 
traveling in 
rangements of both subsonic and supersonic jets were used 


typical shadowgraphs are reproduced and discussed 


Asa direct result of its passage through the jet, the shock wave 


developed a bulge. The growth and behavior of this bulge 


of instantaneous shadowgraphs 


recorded through a sequence 


The transverse spread and the changing profile of this bulge offer 


an interesting visual demonstration of the angle of spread and 


decay of the jet flow. It is shown that beyond an equivalent dis 


tance of about twelve jet diameters the bulge-profile reaches a 
state of similarity 

When the shock and the jet impinged against each other head 
acoustic radia 
shock front. 


re advanced, 


and high frequency 
field 


Some plausible explanations of this phenomenon a 


on, an intense, intermittent, 


tion manifested itself in the flow behind the 


INTRODUCTION 


HE RISING INTEREST in the development of faster 
j prbenente aircraft has accelerated the tempo of 
In this challenging research field, 
of Kuethe,” Corrsin,”: * 
and have 


jet-flow research. 
the contributions 
Pai,” ® Wilder,’ 
indirect bearing on the experiments described here. 


Tollmien,! 
Anderson* some direct or 
In recent years the problem of the intense noise gen- 
erated by modern jet-propelled aircraft, especially the 
uncomfortably high intensity at take-off and landing, 
Generally speak- 
to the 
Therefore a better 


has attracted considerable attention. 
ing, the 
vibrations and/or flow fluctuations. 
understanding of the basic mechanism of noise genera- 


noise originates either due structural 


tion processes in jet flow is essential for the successful 
planning of any extensive scientific program of noise 
During the past few years some valuable 
Powell® 1% 1! 


It has been 


abatement. 
light has been shed on the subject by 
Lighthill,'* !* Ribner,'? and Kovasznay.'” 
observed that the mean jet flow velocity, jet diameter, 
turbulence level, and interaction between any existing 
shock front and turbulent jet flow all contribute to the 
intensity and other characteristics of aerodynamic 
noise. '* 

It is well known that a free jet becomes turbulent 
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wave 


Some 


were 


short distance downstream of its exit.6 To jy 
the a shock 


in a turbulent jet flow, jets were mounted either clos 


only a 
vestigate interaction of front 
to the open end or inside a shock tube and the transien; 
passage of the shock front through the turbulent fh 
field was optically recorded.{ A brief description of th 
main units of the experimental equipment follows. 


THE SHOCK TUBE 


In these experiments the shock tube was used pr 
marily as a convenient source of plane normal shi 
The shock tube 


of a steel box of uniform cross section (approximatel 


waves of known strength. consiste 


tin. inside). To achieve the optimum sensitivit 
to disturbances in the 
the the shock tube 
The experimental observations were made in the 


6 by 


any flow under investigatior 


6-in. side of was laid horizontal 
12-11 


glass section, with top to bottom visibility, mount. 
at 15 ft 
Fig. 1). 


shock-tube theory, the bibliographies in references l6t 


downstream of the diaphragm location (se 


For any detailed information about the basi 
IS are suggested. 
SHOCK SPEED MEASUREMENTS 


The time taken by the shock wave to travel betwee 
the two timing hotwires, 
apart in the upper wall of the shock tube (see Fig. | 
19). The shock 
well-known rela 


mounted a known distance 


was clocked (for details, see reference 


speed Ul’, was calculated, and from 


tions, the shock strength 


pressure behind the shock wave 


pressure ahead of the shock wave 
and other pertinent flow parameters such as drift flow 
and the corresponding density ps, 


shock f 


velocity 2 tempera 


ture 72, and pressure p2 behind the front wert 


calculated. 
The hot-wire response pulse originating at the dow 
was simultaneously used to stop tht 


stream station 


chronograph as well as to actuate the delay system. 1 
decouple the two input loads, a cathode follower wa 
found helpful (see Fig. 2). 


ARRANGEMENT FOR TAKING SHADOWGRAPHS 


The arrangement used for taking sh: idowgr: uphs wit! 
the 


short duration high intensity ‘“‘point’’ source © 


t A few hot-wire oscillograms of this interaction were also tr 


corded but are not discussed here 





1 Turbulent Jet’ 


traveling § 





| 


light 1 
unit se 
shadow 
flow fie 
chroniz 
triggert 
The tri 
is dela’ 
R-C de 


trigger 


Arrange 

Diffe 
sional), 
ties, We 
shock t 
cording 
turbule: 
the trav 
ior the 
these e 


obsery ef 


Arrange: 


In Ar 
ing at t 
weakens 
this limi 
with ins 
bent at 
the shoe 
Shock fr 
the shoc 





| 
Jet| 


lo 1! 


raveling 


er cl § 


ransient 


ont fl 
ni of th 
WS. 


sed pri 
1 shoek 
isiste 
imatel 
sitiviti 
igation 
1zontal 


e 12-in 


1ounted 


on. (see 
1e basi 


es l6t 


yetweel 
listanet 
Fig. | 
> shock 


n rela 


ift flow 
mipera 


it were 


7 dow! 
op the 
m. 1 


rer Was 


hs with 
urce 0 


also re 





























Laboratory Pressure Line 
ee 
Regulator T 
— ae Pressure 
Silica - Jell /_ oo Tube 
6 J 16' Bed 7s a 
Stations Mounting Film f 
Timing Hot Wires or Jets ‘ines & a 
- op Cock 
“4 LU a ~ J 
Compression Chamber 4 Expansion Chamber | | 7 4 
! 
. — 
| Round Jet 
Diaphragm a —~ | KS 
’ ’ Glass 
Section 
cD » 
uzV — Triggering 7) nica 
—— or t) 
Cathode Unit — 
Pr | Pee aes Be stan “ Housing 
/ |Follower | tS 
a oe 
Stabilized A SS 4 ay 
Volt I \ 
oltage Supply \ ’ 7 
\ Delay t 
[Circuit | istert_ovep 5 ** 
| | 
=A 117 V—=—4 Chronograph 
Spark 
sty — 7 
Unit 


a 








Fic. | 
light is shown in Fig. 1. For details about the spark 
unit see reference 20. For recording the instantaneous 
shadowgraphs of the advancing shock front and/or the 
flow field behind it, the spark has to be properly syn- 
chronized. ‘The passage of the shock front over the 
triggering hot wire produces a sizeable response signal. 
The tripping of the spark gap, as a result of this signal, 
is delayed by a pre-estimated time through a simple 
R-C delay circuit introduced between the output of the 
triggering unit and the spark unit (see Fig. 1). 


THE JET ARRANGEMENTS 


Arrangement A 


Different types of jets (round as well as two-dimen- 
sional), operating at subsonic or supersonic exit veloci- 
ties, were mounted just outside the open end of the 
shock tube and a systematic set of shadowgraphs re- 
cording the passage of the shock fronts through the 
turbulent jet flow was taken. A general distortion of 
the traveling shock front always occurred. However, 
lor the shock strengths and the jet velocities used in 
these experiments, nc auxiliary sound emission was 


»)) 


observed (see Fig. 


Arrangement B 


In Arrangement A, the rarefaction wave originat- 
ing at the open end of the shock tube progressively 
Weakens the shock wave (see Fig. 2). To get around 
this limitation, an axisymmetric jet (circular brass tube 
with inside diameter = 3/16 in., horizontal stem 8 in., 
bent at gentle right angle) was centrally mounted inside 
the shock tube (Fig. 3). The jet flow and the advancing 
shock front were both directed toward the open end of 
the shock tube. 


The air was supplied from the labora- 


839 


High Voitage Cable 


Block diagram of the experimental arrangement 


tory pressure line (100 Ibs. per sq.in.), and to reduce its 
water-vapor content, it was made to pass through a 
heavy steel chamber packed with silica jell. The exit 
velocity was controlled by a regulator near the pressure 
line and/or by a stopcock close to the shock tube (Fig. 
1: 

The maximum mean velocity U, 
measured with the total head tube and/or with the hot- 


of the jet flow was 


wire technique at various stations downstream from the 
exit of the jet. It is known that the maximum velocity 


U’,, on the axis of the jet decreases approximately fol- 





Fic. 2. Jet mounted outside the open end of shock tube (ar- 
rangement A Two-dimensional jet (exit 1/8 in. by 6 in 
exit velocity U, = 200 ft./sec., shock strength S = 2.5, delay 


S4 microsec 
uncollimated 


open end corresponds to zero delay), spark light 
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Fic. 3. Jet directed towards the open end of shock tube (ar 


rangement B). A = round jet, B = open end of shock tube, 


C = triggering hot-wire, D = hot-wire probe for recording jet 
shock interaction 





Fic. 4. Interruption of jet flow due to passage of shock front 
Round jet with plug, Ul supersonic, S = 1.58, delay = 114 
microsec., spark light uncollimated 





Fic. 5. Shock jet interaction (arrangement B). Round jet, 
d = 3/16in., U. = 340 ft./sec., S = 1.58, delay = 248 microsec., 
x,/d = 23.6, spark light uncollimated with shock front in line 


with optical axis of divergent light beam. 


‘ 


lowing a hyperbolic power law’ 
U,./U, = 6:5(d/x 


where x is the distance along the jet axis downstream 
the exit and d is the diameter of the circular jet. Th; 
relation can roughly be extended to the present cay 
and the approximate exit velocity (’, can be estimate 
The following observations were made froin the y 


corded shadowgraphs. After the traveling shock fron 


advances just past the jet exit, the jet flow is inter 
rupted for a short interval of time (the same happene 
in Arrangement A). This seems to be caused by th 
sudden rise in the ambient pressure surrounding th 
jet. An effort was made to overcome this difficulty. 

1/4-in. long solid plug with five symmetrically drilk 
holes (diameter = 0.041 in.) was silver soldered in th 
exit of the jet. The higher pressure ratio across th 
plug resulted in a choked jet. Mild acoustic radiatio; 


eminated from the jet flow near its exit. However, the | 


passage of the shock front still interrupted the jet flo 
(Fig. 4). 

Since the jet flow ahead of the shock front has no way 
of knowing the state of the ‘‘shocked” jet behind th 
shock wave, the advancing shock front overtakes ar 
passes essentially through the domain of the normal pr 


shock jet flow. Some of the already well-known aspects 


of the jet flow can be optically demonstrated from the 


sequence of the shadowgraphs of a shock front advane- 


ing through jet flow (Fig. 5). For instance, the progres 


sive increase of the transverse distance of the ‘‘bulge 


in the shock front was exploited to find the angle oi 
spread of the subsonic jet exhausting into quiescent air 


as follows. 

The successive locations x, of the main shock front 
relative to the jet exit, were determined from x, = U 
by using the experimentally measured shock speed | 
and the delay time ¢. The delay time corresponding t 


the arrival of the shock front at the jet exit was taken as 


zero. From the corresponding shadowgraphs, the max! 


mum magnitude of the transverse distance (measured t 


to the point where the bulge profile just becomes tan 

gential to the main shock front) was scaled and plotted 

at the corresponding locations x, of the shock front. 
For x;/d > 10 and with apex at the center of the jet 


exit, a semiangle of spread @ = 10.5° + 1° was meas 
ured. If one laid the transverse distance rp (wher 
/ = 1/2 Imax.), and with the same apex as before, 7% lies | 
on a cone of semiangle (6),, = 6° + 0.5° (Fig. 6 


These values agree rather well with those quoted 
literature for the semiangles of spread of turbulent 
subsonic jets exhausting in quiescent surroundings.® I! 
was observed that in most shadowgraphs the shape 0! 
the bulge was symmetrical with respect to the x axis 


The occasional exceptions seemed to be due to the loca 


distortion caused by some strong local flow fluctuations 

The depth Jax. of the bulge reduced in magnitude the 
farther downstream from the jet exit the shock front 
advanced. In the usual manner of the jet-flow ™ 
vestigations,*® ® 8 r/ry vs. 1/lmax. was plotted for suc 
cessive shadowgraphs of the same (strengthwise, 5 = 
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Fic. 8. Shock-jet-interaction with shock front and jet flow 

directed against each other (arrangement C Round jet, d = 
16 in ,O. = 1.58 





, delay if shock front was 
delay corresponding to shock 
d = 3.46, spark light colli 


a) Supersonic jet U. & 1,200 ft./sec 
to be at jet exist = 414 microsec., 
front location = 378 microsec., x 
mated 


SHOCK WAVE 





(b) Supersonic jet. U, & 1,200 ft./sec., delay if shock front was 


to be at jet exist = 385 microsec., delay corresponding to shock 
location = 377 microsec., x,/d = 1, spark light collimated 





775 ft./sec., delay if shock front was to 
be at jet exist = 418 microsec., delay corresponding to shock 
location = 402 microsec., x,/d = 1.46, spark light divergent with 
shock front in line with the optical axis of divergent beam 


c) Subsonic jet. U, 


fe eS 
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1.58) shock front. It was found that for x, d > 12, th 
similarity plot agreed with the one for velocity dis 
tribution derived by Tollmien! using Prandtl’s mome; 
tum transfer theory see Fig. 7 

Any physical interpretation of the observations 
complicated by the fact that the shock bulge woul 


have spread in transverse direction and decreased 
depth even if the jet was removed part way out 


Arrangement C 


To circumvent the complications caused by the in 
terruption of a shocked jet 
the 
centrally mounted through the open end and directe 
Fig. | The 
jet assembly and the air supply was the same as the one 
The jet flow and the 


interaction has been recorded), 


toward the diaphragm of the shock tube 


described under Arrangement B. 
and therefore its drift flow 
The 
of the jet flow was measured with a pitot tube at ; 


traveling shock front wert 


directed against each other maximum velocity 
number of stations downstream of the jet exit and the 
approximate exit velocity was estimated by relation | 
Both supersonic (Figs. Sa and Sb) and subsonie (Fig. ¢ 
jets were tried. Since relation 1 applied only to th 
subsonic domain of the jet flow, the estimation of th 
supersonic exit velocity is somewhat arbitrary. 

Since the expansion chamber of the shock tube is 
closed by a diaphragm and the jet is directed upstrean 
(Fig. 1), the jet discharge flows, with a relatively low 
back-flow velocity, in the opposite direction to the flow 
in the central jet core 

The following observations were made. The per 
tinent data are listed under each shadowgraph. 

The instantaneous shadowgraphs of the shock front 
reveal that the central portion of the shock front de 
velops a bubble-top bulge at its rear and seemingly there 
isa conspicuous splitting or doubling of the shock front 
(Fig. 8). 


These distortions of the shock front are ol 
viously related to the relative direction of advance 
the shock front and the local jet flow. 

If we choose the coordinates such that the shock front 
is stationary, and it is assumed that the shock front 1s 
plane and the back flow velocity is negligible, then th 
problem reduces to the one with two supersonic parallel 
the 
flow region 1) with different Mach Numbers such that 
M; = (U, + u;)/a;> M, = U 
through a plane stationary shock front. 


streams (jet flow domain in region 3 and outer 
a; > 1, being convected 
From. the 


well-known one- (or two-) dimensional shock analysis 


, 


one may write p2/p; as a function of \/; and p4/ p; as 


function of ./; and a. Appiving the necessary condi 
tion that po/ pi = ps/ps;, one gets sin a = J/,/ 17/3. 

If M3; — Al;—1e., jet velocity u; > 0, then a > 90 
Since /; > 1/1, a < 90 
should exist in the jet flow region. 
shock front U’, 
the jet flow w; near the exit (potential core region) was 


therefore an oblique shock 
The speed of the 
was known. The approximate speed 0 
roughly estimated as described before: 1/,/.\/; was 
calculated and the expected value of the inclination @ 


came out to be approximately 48°. As measured from 


at least till the shock-jet J 
round jet was 


eee 





ema yomn 8 





ELING SHOCKE 





the sl weraph a & 40° (Fig. Sa This simplified 
interpretation ignores the complexity (1.e., three-di- 
mens! lity and additional wave pattern) of the bulge 
shape 

When the turbulent jet, be it subsonic or super 


solic 


the traveling shock front, an auxiliary wave system 


origil 
shock wave [The wave system which resembles an in 
tense acoustic emission seems to have two distinct 
nodes. One is comprised of very high frequency emis- 
sion (almost continuous) and the other, apparently the 
stronger of the two, is distinctly intermittent and sym 
metrical in character (Fig. 8). The apex of these wave 
lets seem to be coincident with apex of the bulge. The 
wave system is not clear in the turbulent domain of the 
jet flow, but is abundant in the outer region of the drift 
flow Due to the subsonic drift flow behind the shock 
the emitted wave system piles up on the main shock 
front at its rear!’ while ahead of the shock front, the 
quiescent region surrounding the jet flow (back flow 
velocity is overlooked) is devoid, within the sensitivity 
limitation of the shadowgraphic technique, of any sound 
field emitted due to the jet flow alone. 

By proper choice of the coordinate system, the shock 
front can be considered to be stationary. In the in- 
stantaneous shadowgraphs, the acoustic emission can 
be analyzed as if eminating from a stationary point 
source located at the apex of the shock bulge and being 
radiated into the moving stream behind the shock front 
Fig. 9). It is assumed?! that the propagation velocity 
normal to the wave fronts is equal to the local speed of 
sound a» in the drift flow region behind the shock front 
which is related to the frequency f and wavelength d of 
the acoustic emission by do = fX. If the shock speed U 
is known, then a2 and the velocity l’, —  (stationary- 
shock coordinates) in region 2 can be calculated, and 
with the help of the vectorial diagram as shown in Fig. 9, 
the wavelength A can be estimated.*! For the case 
when the jet was choked and the shock wave (S = 
1.58) has advanced to a location x/d = 3.5 downstream 


from the jet exit (Fig. Sa 
U, — us)/a2 = M2 = 081, Neve & 0.24 in., 


a. = 1,230 ft./sec., fire B® 60 XK 10° cycles/sec. 


where \,,,, and f,,, are, respectively, the average order of 
magnitude of the wavelength and frequency of the in- 
termittent acoustic wavelets. It was observed that the 
average frequency is somewhat higher (fr, @ SO X 10 
cycles sec.) for the same operating conditions, but with 
shock front at x,/d = 1 and is somewhat lower at dis- 
tances farther downstream from the exit (f,,. A 50 X 


10° eyeles/sec. at x,/d = 13). For the same shock 
strength, the frequency and the character of the emitted 
sound field are strikingly similar in both the supersonic 
and subsonic jets (Fig. 8). However for the super- 
sonic jet (UL & 1,200 ft./sec.), the emitted wave 


system first appeared at x,/d = 24 + | and for the 


subsonic jet gs RA 4640 ft. sec.) at x, ad = 15 = 1. 
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is flowing against the direction of advance of 


tes and propagates in the flow region behind the 
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Since the maximum axial jet velocity l’,, was measured 


at various stations downstream of the exit it is, there 


fore, possible to estimate [ 


for both supersonic and 
subsonic jets at the respective locations at which the 
sound emission first clearly manifested itself. Such an 
estimate revealed that for both supersonic and subsonic 
jets, the corresponding local / A 340 + 30 ft./sec. 
Based on normal shock considerations, the drift flow 
velocity u# = 38S ft./sec. 

The drift and the jet flows impinge against each other 
and develop a common stagnation point close to the 
shock front. The shocked jet flow 1s forced to reverse in 
direction and is continually being pushed outward in 
the radial direction resulting in a clear swelling of the jet 
flow (see Fig. 8). The drift flow is called upon to adjust 
itself around this swelling and any resulting imbalance 
between the pressure in the jet-flow region and the 
outer-flow region seems to contribute toward the ob 
served auxiliary wave system whose intermittency is, 
most likely, due to the oscillating or buffeting of the 
shock segments constituting the bulge. These oscilla- 
tions could be caused by the intense velocity or pressure 
fluctuations convected through the shock front.’ 

The characteristics of turbulence—1.e., its intensity and 
spectrum—will undoubtedly have a direct influence on 
the nature of the emitted sound field. The present 
data were not sufficient to estimate the level of the 


emitted acoustic power. However, it seems apparent 
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from the shadowgraphs that the level of emission in- 
creases as the shock front progresses toward the jet exit. 

All these observations were recorded with the shock 
wave of the same strength (S = 1.58). The shock 
strength may have a direct effect on the characteristics 
(location of origin and intensity) of the acoustic radia- 
tion. 
the operating conditions of these experiments, Ay», A d, 
Therefore, the 


In passing, it ought to be pointed out that for 


where d is the diameter of the round jet. 
dependence of \ on the diameter of the jet nozzle (or 
local diameter of the jet flow) seems likely. 

The other mode of the observed acoustic radiation 
consisting of very high frequency components possibly 
is due to the convection of the corresponding small 
turbulent disturbances 


scale (higher frequencies) 


through the shock front as expected by Lighthill,'* '% 


Ribner,'* and Kovasznay." 
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Reply to ‘‘Further Comments 
Diffuser Instability’ ”’ 


on ‘Supersonic 


Robert L. Trimpi 
Aeronautical Research Engineer, Langley Aeronautica 
Laboratory, NACA, Langley Field, Va 
April 29, 1957 
- statements! made by C. L. Dailey indicate incorrect 
interpretations of the comments? and original publications 
of the writer on the subject of ram-jet buzz his brief di 
cussion will attempt to clarify a few of the points so misinter 
preted 
Phe validity of the application of a quasi-one-dimensional un 


steady flow approach to the buzzing problem is questioned in 
Consequently an examination of the 
ippro 


permits cre 


reference | imphications 


which are associated with the use of such an ich is in 


The 


tional area variation but considers only 


theory 


ISS-SECC 


order quasi-one-dimensional 


one velocity component 


(axial) and averages fluid properties at a cross section. First, 
in regards to the use of quasi-one-dimensional theory inside the 
ram-jet, many experiments have proved that the cross-sectional 
area versus distance distribution of practical ram-jets is such 
that a (both 
gives good correlation 
of the 


inside the ram-jet 


one-dimensional analysis steady and nonsteady 


Thus, regardless of what happens ahead 


ram-jet cowling, the one-dimensional approach is valid 


Now the difference between a three-dimensional and a quasi 


one-dimensional approach lies principally in the normal vari 


stream tube In particular, 


the 


ation of flow parameters inside a 


when the quasi-one-dimensional approach is applied to 


stream tube extending from the ‘‘normal shock”’ to the cowling 
the 


gradients of pressure and entropy (the latter is mainly boundary 


lip at instant of incidence of inlet instability, the normal 


layer effect) are ignored by averaging across the stream tubs 
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flexible walls (to account for spillage around the cowling) can The phenomenon of inlet buzz is one of the more gross, a) f 
yield qualitative information regarding the initiation of buzz hence more obvious, aspects of the general subject of propulsj 
It might be noted that, contrary to the contention of reference 1, system transient behavior. While there is no point labori al 
no negative values of this stream tube area are required to handle the description of this phenomenon further, it is worth no te if 
outflow but merely an increased timewise growth of area that its cause is not yet thoroughly understood and falls in ¢ j 
In reference 1 the statement that “. . . a buzz cycle can be general class of things to be cleared up before transient syste 
generated by the spontaneous ‘outward’ motion of a shock in a analysis will enjoy the same rigor that is common in stea 
one-dimensional channel . . ” was erroneously attributed to flow analysis 
reference 2. No such statement was made in reference 2; in p 
fact, the following excerpt from reference 3 expresses the views | 
of the writer in this respect: ‘The entrance of the vortex sheet, as 
increasing separation, or boundary laver growth (behind the + 
shock), either on the central body or on the inner surface of the On the Variational Theorem for Creep 
cowling, generates a compression wave which forces the normal 
shock outward.” T td Bien 
The opinion? that the “fill up” portion of the pressure-time Assistant Professor, Department of Aeronautical Engineerin: 
curve of a buzz cycle was a series of steps rather than a relaxation Massachusetts Institute of Technology, Cambridge, Ma ¢ : 
cycle was also challenged.!. The writer interprets a relaxation July 8, 1957 of E« 
curve to be a smooth exponential curve and in addition the terms It 1s 
“simple exponential solution” are used in reference 5 to describe N A RECENT NACA PUBLICATION,' Sanders, McComb ng 
the solution to the theoretical pressure-time variation. The I Schlechte presented a variational theorem for a body under a 
experimental pressure-time curves of references 4 and 5 do in- going creep. The authors stated that the integral to be varied w 
deed show such steps and not a smooth exponential trend obtained by trial and error. It appears to the present writer th 
The analysis of reference 1 is inconsistent in regard to another the same integral can be derived in a more rigorous manner using | 
aspect of the “fill up” cycle. In one paragraph it is acknowl the generalized Reissner's Principle for bodies with initial stress 1003 
edged that excellent pressure versus time results may be ob- and initial strains.2 W 
tained by the use of the quasi-one-dimensional theory which As shown by Wang and Prager® the present boundary proble 1.1.1 
includes the reflections both at the exit nozzle and at the swal- may be stated in the following manner: In the interval between | ‘ a 
lowed shock to satisfy prescribed boundary conditions. Then the instant ¢ and dt, the surface stress is assumed to be change ing I 7 
in the next paragraph the question is raised whether the simple by d7; on S; and the displacement by da; on S». Creep deforma — 5.1 
mechanisin of reflected compressions are able to build the plenum tions will occur during this interval. What are the corresponding 
chamber pressures back up. Such contradictions should not changes of stresses do,;; and the corresponding displacements d 
appear if comments result from careful consideration of the throughout the body? 
problem under discussion For this problem the state of stress at the beginning of the tin 
interval may be considered as initial stress and the creep strains 
REFERENCES : : : ae ; 
during the interval dt may be considered as initial strains 
Dailey Charles L., Furtier Comments on “Supersonic Diffuser Insta For a body with initial stresses and initial strains the integr Ont 
bility Journal of the Aeronautical Sciences, Readers’ Forum, Vol. 24, No s s 
+. 19 OU Renee, 1008 7p to be varied may be written as for 
? Trimpi, R. L., Comments on “Sup mic Diffuser Instability,’ Journal of P . : 
the Aeronautical Sciences, Readers’ Forum, Vol. 23 6, pp. 611, 612, Z=™ f,. (a); + dojj)de;, iW dojjei;| dl Y/alt 
June, 1956 e e A 
Trimpi, Robert I 1 Theory for Stability and Buz ulsation Amplitude | (7. 4 dl du.dS — | (du adaVT7T; + d7 dS 
in Ram Jets and An Experimental Investigation Includ Scale Effects, NACA Js JS 0 
Report 1265, 1956 (formerly NACA R mt L53G28 ites 
' Trimpi, Robert I {n Analysis of Buzzing in Supersonic Ram Jets by a 
Molified One-Dimensional Nonstationary Wave Theory, NACA TN 3695 oi; = stress at the beginning of the interval dt 
july, 1066 (formerly NACA RM LSZAIG, 1962) do = increment in stress during dt 
Dailey, C. L., Supe tic Diffuser Instability, Journal of the Aeronautical ; HI 
Sciences, Vol. 22, No. 11, pp. 733-749, November, 1955 de;; = increment in strain during d/ T.. 
ej’ = initial strain eae 
HW” = a function of incremental stresses defined by de;; = lites 
+ oW/ddo,, ~ 


Reply For problems of finite elastic deformations the strains can | - t 


expressed as 











\ fr 
Cc. L. Dailey € = (1/2)(u tT U T 2 il - 
grange 
Research Associate, University of Southern California, Los Angeles nin ; 
Calif ence, plates 
July 23, 1957 de = (1/2 du ° oe du T @ du T U du T i F re 
the thi 
du du ve 
rt IS HOPED that this discussion will not be continued pub- Teue ; saa 
‘ HE. al : rhe integral wp, may be separated into two parts, indwi 
licly The many statements and _ clarifications thereof 
made to date must surely be sufficient to permit the interested rn sy TT 81 : si 
. . . ° . . ini 
reader to arrive at a point of view o his own acai me 
It is difficult, however, to be entirely satisfied with the present e 
' ; : , 1 
state of knowledge of propulsion system transient analysis rp! = { (1/2)0;;(du 1 du L du Loy, du dV - 
vans *y cs - oo : J : : ‘ prool 
Mr. Trimpi’s treatment of inlet buzz affords a convincing ex- . P - 
ample of the extent to which muddled thinking can be seriously | | as — f (du du;)7T; dS (3 vgs 
° . eve ] 
regarded as an apparently meaningful explanation . oa R 
° . . ° ° fe AC 
It is sobering to realize that no method exists at present for =" = { ((1/2)doi(du;, ; + du, La du iK~ du 4 a 
predicting either damping and frequency of coupled engine-inlet Ji 
modes of oscillation or internal pressure noise spectrum of (1/2)o;jdu du W — dajje d| 
, . ‘ste ve *h inf ation i bvious practical Fm ¥ 
coupled systems. Yet such inform ition is of obvious practic dTdudS (du: — da:)dTdS (6 ab 
importance toa propulsion engineer. J Si ny 7 
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gross, and & It wn that the variation of zg’ vanishes because the for instance, that problems of plane stress and membrane deforma 
t 
ropuls; ‘nti o;; satisfy the equations of equilibrium and tion will lead to designs whose weight is an absolute minimum 
laboring hour mditions under the criterion derived in reference 1; in the latter paper, 
th noting If t te of creep strain is expressed as a function of stress in the present one, the minima derived are the relative ones in the 
IS in the § onl in the case of secondary creep, where é = f( J? sense of the calculus of variations, and which of these is absolute 
it syster S.. the creep strain may be written as €;;°° = €j;“ dt. Further has to be decided in each particular case. For beams and plates, 
In steay mort rements in stresses and displacements may be ex however, the method of references 5 and 6 can do no better than 
terms of the corresponding time derivatives, do;; = reference | No general method for determining an absolute 
| 6:4 ] 7, dt, du; = u; dt minimum is vet available 
j I iting integral rp” Is Further progress in the field is likely to accrue from the two 
. o parallel lines of attack—-the method of limit analysis (see ref 
) ae a ' : 62 — 
i r | Oi \¢ 1/2 )oiju ’ i Gije d\ erences 5, 6, and 7) on the one hand, and the calculus of variations 
e . ° 1 1 ‘ 
i ° bg ) reference 1 and the present paper) on the other 
i Ti; dS { (it ‘)7; dS dt? (7 
» * \ 
: ; , SOLUTION OF THE PROBLEM 
wher is of the same form as I] except that da,; is replaced by 
neering o I integral to be varied is now reduced to that in the brace Phe sandwich shell of length 21, radius R, is to have a core of 
; of Ea. (7 It is equivalent to the integral I given in reference 1 negligible weight and prescribed constant thickness //, and the 
| It is « r that the variational operations should be applied to thickness / of its faces 1s to be designed for minimum weight. The 
' } incremental stresses and displacements only—i.e., to time shell faces are to flow in the plastic-rigid mode; the vield condi 
mp . . ° 
derivatives of quantities but not to the quantities themselves tion can be written 
Vv under 
‘ * is 
iried w REFERENCES P( Ne, Neo 52) | 
iter that J ; 
3 Sanders. J. L.. Ir.. McComb, H. G., Jr., and Schlechte, F. R., A Varia where V,, Vy, ./, are the stress and moment resultants and o is 
1er using § , F tn Plot j ’ ’ : e 
\ tional é vith A sic Plates and Colun NACA TN the vield stress in pure tension. The von Mises criterion, general 
| Stresse 1003. Mav 957 . itotl It t . 
; iz » three resultants, is 
Washizu, K., On the Variational Princif Elasticity and Plasticity — i 
probler M.I.T. Aeroelastic and Structures Research Lab. Tech. Report 25-18, March F= N,-4 N.? a W 1 3/4H2)M o 2h? » 
betweer ge > = » Effect a : . . 
: Wang, A. J., and Prager, W., Thermal and Creep Effects in Vi Harden rhe generalization of Tresca’s criterion for shells was given inde 
change ing E Plastic Solids, Journal of the Aeronautical Sciences, Vol. 21, No — 
12 , pendently 
leforma 5, pp. 343, 344, 360, May, 1954 : ’ . 
For homogeneous functions F, Euler's theorem applies; when 
ponding . 
3 k is of second degree, 
ents du 
j V,(OF/ON + NZ(OF/ON,) + M,(0F/OM 1/2)F 3 
er + 
the time | a . 
rhe rates of extension and curvature art 
) Strains 
q en = w/R, é du/dx, h d*w/dx } 


integra On the Minimum Weight Design Problem where uw and w are rates of displacement and x ts the coordinate 
for Cylindrical Sandwich Shells* in the axial shell direction. From the theory of the plastic po 


tential, the following flow rule may be associated with the vield 
Walter F. Fre berger condition (1 
Assistant Professor, Brown University, Providence, R.| o 4 
iS (1 Aeonpags >i igaais eo = adF/ONy, « = adF/QN,, « = adF/d 5) 
luly 8, 1957 


where a is a function of position, expressible in terms of the yield 
condition and the energy dissipation rate 


INTRODUCTION hell surface is 


f 


The rate of energy dissipation per unit are 


4 pa PROBLEM of the plastic design of a structure for minimum oe ae oe ae a ) 
weight is a problem of constrained optimalization. The un 
known functions—stress-resultants and thickness-parameters for 


] 
de = <a : 
ll 


plates and shells—-are constrained to lie in a domain defined by D 1/2) aF vo ” 


the plastic vield condition of the material and the equilibrium itl 
and thus 








can | equations of the structure; the objective function to be optimized 
s the weight of the structure, usually in integral form OF/ON; h?/2D)(w/R), OF/ON h=/2D\(du/dy 
\ frontal attack on this problem is most naturally by the La dF/aM h?/2D) (d2w /dx? 2 
. grange multiplier method of the caleulus of variations. For 
plates, both solid and sandwich, this has been done in reference 1 Phe equations of equilibrium of the shell are 
| For cylindrical shells the problem is complicated by the fact that Vv." 4 N./R re ind \ . Q 
the thickness parameter no longer enters homogeneously into the : 
vield-condition; this homogeneity can, however, be preserved for where primes denote differertiation with respect to x. and p and q 
sandwich shells, and this latter problem is solved in the present are the radial and axial loads on the shell per unit area 
} paper The weight of the sheil will be a constant magnitude (the core 
\n intuitive criterion applicable to structures obeying a piece plus a quantity proportional to 
Wise smooth yield condition led in earlier papers to a °] 
aV- minimum weight design procedure for shells. Although formal | Pa dx 
proof of the optimality of this procedure is as yet lacking, the 
S (5 heuristic criterion has been successful in all cases studied How The variational problem thus reduces to the following: to 
ever, it is by its nature not applicable to smooth yield surfaces minimize 
Recently, the theorems of limit analysis have been applied to 7, 
. i general discussion of minimum weight design.’ ® It is shown, { P hax 10) 
. 
ts presented in this paper were obtained in the course of re under the conditions 
>» (t earch sponsored by the Ballistic Research Laboratories of Aberdeen Proving 
roun r Contract DA-19-020-ORD-798 F(N:z, Ve a. V 11 
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V + (1/R)NZ —p=0 12 
\ —g=0 13 


With Lagrange multipliers A, uw, v for Eqs. (11), (12), and (13), 


respectively, the Euler equations become 


1 — 2yh = 0 14 
NOF/OM - ps’ = 0 15 
\MOF/ON,) — v’ = 0 16 

MOF/ON,) + (u/R) = O 17 


Eliminating A and w from Eqs. (14), (15), and (17), we obtain 
R/h)(OF/ON,)|"" = (1/h)(OF/OM 1S 
and substituting Eq. (8) into Eq. (18 


hw’/D)"’ = hw'’/D 19 


Putting D/h = A, the rate of energy dissipation per unit volume, 


Eq. (19) can be written 


or, integrating, 
w?/A?)A’ = B 2») 


where 8 isa constant of integration. It can be shown by consider- 


ing the boundary conditions at the ends of the shell, that 8 = 0 
at x = +L, and thus for all x 


We conclude, then, that we must, for minimum weight, have 
A’=0 or A = D/k = constant 21) 


A relative minimum for the shell weight in the sense of the 
calculus of variations is thus fouiid for sandwich shells designed 
to give constant rate of energy dissipation per unit volume of face 
material. In view of the local plane stress conditions prevailing 
in the faces, this condition will also give an absolute minimum 
and thus an overall optimum design according to reference 6 

For the von Mises yield condition (2), Eq. (18) gives 


h]’’= (8R/2H?)(M,/h) 22 


and it can readily be shown that in this case the functions 1/,(x) 
and h(x) are determined by the following two differential equa 


tions: 


—2M,""" + 4t°M/" — 2(t"? — t'")M,"' — (38/H?R2)M, = 


(M,'')? + d(x) M."’ + (8/4H?)M,? + e(x) = oe? 


where we have put ¢ = log h, and the quantities a, 5, c, d, e are 
known functions of x depending on the applied loads p(x) and 


qQ\x 
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Past a Flat Plate 
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r i SHE PROBLEM of the laminar boundary layer on a flat plat 
in a stream of vorticity w, has been the subject of a rec 
series of notes by Li.' 2 The appropriate bound 


equations for steady incompressible flow with no pressure gr 


dient are 
uU(OuU/OX) + V(OuU/OY) = v(A°U/O4 


Ou / Ox) 





with boundary conditions 





t+=v=0 at y=O0 | 
! 
“u=t ayy at =U 
if 
where “# and we are constants ' 
In place of condition (4), Li adopts in his first note i 
' 
u-—> ut tT Wov j ' 
as ' 
as the condition at the outer edge of the boundary layer. 0; 
writing 
n= VV Mo/wx FE = wo wx/u f 
Li obtains a solution im series for the stream function y in th 
form 
¥ = V morx ji foln) + Efi(n) 4 | 7) | 
i 
1 
where 
u = (d¥/dy) v = —(dy/da 8 
in accordance with Eq. (2). On equating like powers of € it 
Eq. (1), 
of,’ + fif = () ) 
ele SE Nal ff ort a ] 
the boundary conditions being 
f, =f =0 f, =f,’ =0 at n=0 
fy’ > 1 fi)" ~>1 as nm 12 
The required solution of Eq. (9) is the well-known Blasius fun 
tion, and Li gives a solution of Eq. (10) obtained by an approxi 
mate method. 

In a later note® Li retracts his previous solution. He points 
out that, as » > tomy c, where c is a constant (with th 
known value 1.7208). It follows that, at the edge of the bound 
ary layer, 

— 2)¢ V uov/x 13 
and hence, in Eq. (1), P 

v(Ou/Ay) — (1/2)woe VY uov/a 

Li deduces that this term must be balanced by a pressure gra 
dient since the other terms of Eq (1) become zero outside th 
boundary layer, and on introducing this he obtains a modiiit 
Eq. (10) with —c on the right-hand side. He does not carry 
a solution of the modified equation 

However, this concept of an induced pressure gradient 1s 
r . . ‘ . ; 
fallacious one. It arises from the use of the boundary condit! 

5) in place of (4). The correct condition at the outer edge 
the boundary layer is 
u(Ou/Ox) + v(du/dv) = O 
since neither pressure nor viscous forces are significant, 
so u is constant along each streamline. Since the streamlines 
* Lecturer in Mathematics i 
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contributions to the velocity and the shearing stress 


due to stream vorticity 


e displaced by an amount equal to the displacement thickness 





§, this implies that, at the edge of the boundary layer, 

u > Uy + wooly — 4) 16 
In (5), Li requires that w shall be constant at fixed y. This in 
volves an increase of velocity of fluid elements (for wy > 0), and 
hence an accelerating pressure gradient is called for. There is 
no mechanism within the boundary layer to induce such a pres 
For all boundary-layer flows (16) is the appro 
(5 " but the dif 


sure gradient 
priate type of boundary condition rather than 
ference between the two is apparent only when the external 
stream has vorticity In fact, the solution of the boundary 


laver Eq. (1) under the condition (4) automatically satisfies 


the revised boundary condition (16) For large n, fo ~ 7 ( 


If f, ~ (1/2)n? + dn + e {so that Eq. (12) is satisfied], substitu 


tion in Eq. (10) shows that d = —c. Hence, from Eq. (7 
u = (O~/OVv) > to + @w ly CV wx/to) 17) 
and thus 
u(Ou/OX) — (1/2)woe VY uov/x 18) 


cancelling with the value of Eq. (14) and so satisfying Eq. (15 





This is only as was to be expected. The boundary-layer equa- 
tions allow for continuity and for the fall to zero of viscous 
stress with increasing distance from the surface, and so the con 
ditions at the outer edge of the boundary layer must be correctly 
Eqs. (9)-(12) are indeed the 


described. Thus Li’s original 


ippropriate ones 


e method used by Li! to solve Eq. (10) is to replace it by 





two approximate equations, valid for small and large values of 


respectively, and to fix the undetermined constants by match 


ing the solutions at » = 3.0. For 7 small, fo is replaced by 
1/2)an? (where a = fy"(0) = 0.33206), and for » large, fy is re- 
placed by n c. These approximations are far too crude for 


to expect the solutions to have any range of overlap. In 
iddition, the solution Li gives for small y does not in fact satisfy 

the appropriate equation 
\ far more accurate procedure is to integrate Eq. (10) nu- 


merically in a straightforward manner For 7 small, 


= a(n?/2! 1/2)a2(m°/5!) + (11/4)a%(n8/8!) +... 19) 


ind hence from Eqs. (10) and (11) 


= 92 


= } 2/91 1/2)a(m*/5!) 4+ 23 /4)a?(n* 8!) +4 20) 


vhere b is some constant With 6 = 1, f, and its derivatives were 
caiculated at a set of small values of n, and the solution was then 


continued by use of Milne’s integration formulae,* 


1/3)h(2) vy’ +2y), 2 —w = 


Using intervals of 0.2 in » and working to four decimal places 


ALDERS’ 





FORUM 


with a desk calculating machine, the integration was swift and 


presented no difficulties. It was found that ipproached the 


limiting value 1.258, showing that the correct value of # is 


1/1.258 = 0.795 The value of obtained bv Li was 0.278, 


very serious underestimate 


The variation with 7 of f;’ and f," is shown in Fig Since 
Uo = Uni f n r 27 Y t : 21 
and 
rT = p(OUu/ ON = ply” VY v/ tox if rT 2 ”n T } 22 
these functions give the contributions to the velocity and shear 


ing stress in the boundary layer, due to the vorticity in the 
stream. In particular, the skin-friction 7) on the plate is given 


by 


=z ().332 ply VV v/tyx F 0.795 pw + 93 


The first term represents the skin-friction due to a uniform stream 
u \lso outside the boundary layer r = pw rhus the pres 
ence of the boundary layer reduces the effect on the plate of the 


shear in the external stream 
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ytecnnic 


¢ kw PROBLEM Of simple shear flow past a flat plate in an 
incompressible fluid of small viscosity is the subject matte 


of Glauert’s note above He has made some interesting observa 


tions on the results by the writer published in this JoURNAI! 
For the sake of reaching a good understanding of the problem, it 
seems proper to submit the following comments 

The writer is not convinced that the demonstration given in 
Glauert’s note proves the nonexistence of an induced pressure 


gradient in the present problem In fact, for large n, if we take 


nn“ + Gn -+ ¢ 


ada ae A a t's . Ofte «x , 1 
vields ~ 2)1 +e(ie.,d =0 Hence, for large 


higher order tern 


Consequently, 


u(oQu/ox) + Ou /Oyv) —> ats —™ (c/2)wo V yin 3 


Now, in the inviscid free stream the Euler’s equation of motion 


in the x direction (allowing for induced pressure effects) is 


u(Ou/Ox) + v(Ou/dOY) = 1/p) (Op/On } 
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U*Ue 
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4° x 
U=Uoe Oo U=0O 
Fic. la. 6* = displacement thickness 


Comparing Eqs. (3) and (4) we obtain the induced pressure as 
(1/p) (Op/Ox) = v5w (5) 


Therefore, Glauert’s type of demonstration is here used to demon 
strate the consistency of the induced pressure gradient concept 
with Eq. (1) above. This seems to show that Glauert’s Eq 
(15) represents a more restrictive inviscid flow equation than 
Eq. (4), (it disallows the induced pressure effects). He can 
obtain the necessary consistency only by adopting his revised 


boundary condition (16) which reads 


u—> Uo + wooly — 5*) 6 
for » > The physical interpretation of Eq. (6) can be de 
scribed as in Fig. la. For x < 0, the streamline corresponding 
to u = wm lies on the line y = 0. For x > O, this streamline 
becomes y = 6* on which wu = u This approximate treatment 
implies that there is no mass flowing into the region 0 < y < 6* 


In incompressible flow this approximate treatment yields reason 
able accuracy, as the writer believes, only when the displace 
ment thickness 6* is very small—i.e., when the displacement 
effects are negligible. This physical picture leads to the formula 
tion reported in the writer’s first note (Glauert’s reference 1 
As the displacement effects are not always negligible, an alter 
native approximation is, as Glauert also recognized, to take u 
as constant at fixed y (Fig. 1b). This physical picture leads to 
the formulation reported in the writer’s third note (Glauert’s 
reference 3 In support of this latter formulation, one may say 
that since the viscous effects are confined within the boundary 
layer, 0) y < 6, a streamline from the free stream at y = 6 will 
maintain its constant velocity until it enters the boundary 
layer; therefore, the boundary condition (2) should be applied in 
integrating the boundary-layer equations 

The approximate method used to solve the boundary-layer 
equations in the first note by the writer was essentially the method 
used by Kuo For » ~ 0, the appropriate equation was further 
simplified by dropping higher order terms in y as was done in 
Kuo’s paper. In this manner, the solution (18) in the writer’s 
first note was obtained. The writer himself was not quite 
satisfied with this crude method of computations \bout a 
year ago, Mr. R. E. Duffy helped the writer recomputing the 
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solution on a REAC machine. The results were mack 





able to certain groups but were not published so far because 
the writer’s opinion, the formulation given in his first 
provides only partial answers to the problem. However, t 
writer is happy to report almost complete agreement betw 
the REAC results and Glauert’s data. Mr. Duffy has also 
puted the velocity profile inside the boundary layer for vari 
values of free stream vorticity, using these REAC results 

The nonhomogeneous Eq. (1) has recently been numeri 
integrated with the help of Mr. W. Lick. The results are sh 
in Fig. 2. It is noted that here f,"(0) = 3.11 which is consi 
ably higher than Glauert’s value. This seems to indicate 
importance of the induced pressure term 

Recently, Sakurai? investigated the problem of shear flow al 
a flat plate with uniform suction. Sakurai’s Eq. (8) reads 


his notations, 
p[v(d?u/dy? vp(du/dy)| = Op/Ox = a@ = —w wtp 

In the case of uniform flow past a flat plate with uniform s 

tion, the governing equation is 
p|v(d?u/dy*?) — u(du/dy)| = 0 = Op/ox 
This example illustrates clearly that the effect of a free-strear 
vorticity together with v) is equivalent to a pressure gradient 
the x direction 
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A Series-Expansion Method for Computing 
Random Gust Responses* 


Kenneth A. Foss 

Research Engineer, Aeroelastic and Structures Research Laborator 
Massachusetts Institute of Technology, Cambridge, Mas 

July 8, 1957 


— RECENTLY PRESENTED a general unification of t 
statistical approach to the dynamic response of aircraft 
turbulent air. In a follow-up note in the Readers’ Forum 
Richardson? pointed out that, although Ribner’s formulation 
simple and straightforward, the actual computation of the st 
tistics of a response becomes exceedingly complex and unwiel 
when the two-dimensional characteristics of turbulence are ¢ 
sidered in conjunction with a lifting-surface aerodynamic theor 
A series-expansion method which may alleviate some of the cot 
putational difficulties will now be described 

When an airplane is flying through homogeneous turbulen 


* This method was developed under research contract AF 33(61t 
sponsored by the Aircraft Laboratory, Wright Air Development C enter of th 


USAF 
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the s density of the time variation of a stress o@ In Its 
fe ) ture can be written as 
. 
f= Velde, be) POe(ks, heldks, (wo = Uh 
m 
ert spectral function for the variations of vertical gust 
elo horizontal plane; and IT, is a transfer function 
it Y computed as the steady-state stress response to the 
two-dil ional sinusoidal upwash pattern 
Ut—x _ pikey piw[t—(x/l » 
i, ~) 
| yt for very simple cases, Eq. (1) cannot be integrated in 
closes 1. Furthermore, the asymptotic behavior of the 
integrand usually resembles that of a damped sinusoid making 


yerv diflicult any numerical integration technique However, 
since practical cases of gust response involve large-scale atmos 


nheric turbulence, the effect of spanwise variations on Eq. (1 


nsidered as a small correction to its one-dimensional 


p ye 

equi t 
Qa\w Vo(ki, 0) * bela } 
ere @ is the power spectrum of w(x, 0, ¢) and is related to 

Mi) k,, k by 

. 
$,(w) = (2/l / x(k, ke) dh 
+ 

This suggests the expansion of Ty (k), Rs) into an infinite series 
ibout I’, (Rk, 0) which can be accomplished by expanding e’**” in 


Eq. (2) intoa series 
lhe Maclaurin series for this exponential when applied to Eq 


1) would involve integrations of the type, 


f ko" bw(Ri, ko) dk 


ind the asymptotic behavior of ¢,. is such that these integrals may 


not exist for the higher powers of k Thus, in this method e” 


will be expanded into a series of functions each of which ap 
proad 
One such set of functions can be obtained from powers of 


hes a finite value as k» becomes infinite 


x k = K 


where « bk./2, and \ is an arbitrary constant 3v substituting 


the Maclaurin series for 


into the Maclaurin series for e” the required series expansion 


of Eq. (2) is obtained in terms of powers of x(k 


If we now confine our discussion to symmetric responses of the 


urplane, only the real part of e’”*’ is effective; and Eq. (2) can be 


expanded into 





C.F, g~ x/l + 3in(v) x2” (R 6 
wher 
} a 4 % = $1) — (n?/2!d\? 
B; = —[(n°/6!) — (2n4/4!A2) + (9?/2!A4)], 7 
nd where » = 2y/b. Since power spectral methods are restricted 





to linear systems, the transfer function obtained from Eq. (6) can 


ve Written as 


Where each [’,, is obtained as the response to 


y,t = Bly eur 9 


e generalized forces due to each w,, as well as those due to dis 


ements of the generalized coordinates can be obtained from 
the linearized unsteady-flow theories used in flutter analyses 

By squaring the absolute magnitude of Eq. (8) one obtains 
trom E¢ ] 














A “+ 
$,(w 
= 2 
L/U)w Term 
2 Terms 
3 Terms 
4 and 5Terms 
Fic. 1 Series expansion of span-averaged gust spectrum 
Qq\W le k ) Py \@ T S Ww) y ¢ vy 
a 
, l 
where represents the collecting of terms with like powers of 


x, and where 


For a particular type of homogeneous turbulence each y,, can 
be evaluated once and for all 
The first term of Eq. (10) is the output spectrum for one 


dimensional turbulence, and the remaining terms of the series are 


a correction which accounts for the spanwise variations of up 


wash. In large-scale atmospheric turbulence this correction is 


small and the series should converge rapidly; and the amount of 


computational work required is roughly the same as that for solv 


ing an equivalent one-dimensional problem as many times as 
the number of terms needed in the series 

To illustrate the convergence of this method, let us consider 
the simple problem of computing the lift on a rigid uniform wing 


} 


dy been shown by Liepmann, 


by strip theory It has alre: 
Diederich,‘ and Hakkinen and Richardson’ that the power spec 


trum of this lift can be written in the form 


k ,0 * oy (aw 12 


OL(w I 


where ¢, is a span-averaged gust spectrum which, for the purpose 


of the present method, can be defined by 


o,(w) = 2 uf [ 1/b J : e* ay]"o ki, ky) dk 13 


By making use of the approximate expression for isotropic 
turbulence given in reference 3, this series-expansion method 
with A = 4 


For this example, the series has converged for 


was applied to Eq. (13); and the results are shown in 


Fig | ill prac 
tical purposes with four terms (the first term is the one-dimen 
when the of-turbulence 
This value of 6/L 


than would be experienced in atmospheric turbulence, 


sional gust spectrum span to scale 


ratio, b/L, equals 0.6 is probably larger 
ind other 
calculations have shown that this example converges with only 


» 


two terms of the series when //Z is less than 0.2 
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Effect of Air Temperature on Vortex Shedding 
Frequency of Cylinders 


William D. Clifford* 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, Ohio 


July 15, 1957 


age S INVESTIGATIONS have correlated the Strouhal Num 
ber of evlinders as a function of Reynolds Number. The 
Strouhal Number is S = fd/U, where f is the vortex shedding fre 


quency, d is the cylinder diameter, and L is the free-stream ve 
locity. Rayleigh! first suggested that the Strouhal Number is a 
function of the Reynolds Number based on cylinder diameter 
Among others, Kovasznay? and Roshko*® have found consistent 
correlation between Strouhal Number and Reynolds Number in 


three distinct Reynolds Number ranges: (1) a 40 to 150 Reynolds 


Number lower range of the classic viscous vortex street region, (2 
a 150 to 300 Reynolds Number transition period, and (3) a 800 to 
100,000 Reynolds Number upper range where shedding vortices 
are quickly dissipated 

At the time of this investigation, previous research had not in 
vestigated the effect of temperature on Strouhal Number. Stream 
velocity and ecvlinder diameter have been common variables, 
while stream temperature has been held constant at room condi 
tions. It was the purpose of this investigation to determine 
whether the Reynolds Number completely accounted for the ef 
fect of air temperature on Strouhal Number 

Apparatus used in this work consisted of an air supply metered 
through a pressure regulator and rotameter. The air, heated by 
in electric heater, then passed through a calming chamber and a 
convergent nozzle into a 0.875-in. by 1.5-in. long cylindrical 
duct exhausting into the room. Temperature was varied from 
room temperature to 300° F. and was measured by a thermo 
couple located at the exit of the duct. A 0.0625-in. diam 
eter rod was placed 0.125 in. downstream from the nozzle 
throat. Shedding frequencies were sensed by a hot-wire anemom 
eter at a position 4 diameters downstream of the cylinder, and 
the signal was relayed through an amplifier into the vertical axis 
of an oscilloscope. Frequency was determined by connecting an 
oscillator to the horizontal axis so that the vortex shedding fre 
quency could be determined by the Lissajous figure method. The 
effect of cylinder blockage on the approach stream velocity was 
compensated for by the method described in reference 4. 

As shown in Fig. 1, complete agreement was not obtained 
with the data of Roshko and Kovasznay. Consistent variations 
of the Strouhal-Reynolds Number correlation were observed as 
the air temperature was varied from room temperature to 300° F 
From these observations, Reynolds Number apparently does not 
completely account for the effect of air temperature on Strouhal 
Number. The results of this investigation suggest further re 
search into effects of temperature and pressure on Strouhal num 


ber 
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Summary of Low-Prandtl-Number Heat- 
Transfer Results for Forced Convection ona 
Flat Plate 


E. M. Sparrow and J. L. Gregg 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, O} 
July 19, 1957 


INTRODUCTION 


| pegs ii the heat-transfer characteristics of liquid met 
has been stimulated by nuclear-reactor applications. 1 
present note is concerned with laminar foreced-convection 
transfer over a flat plate for the Prandtl Number range a 
propriate to liquid metals (0.002-0.03 Exact solutions of t 
boundary-layer energy equation* have been obtained for tw 
situations: uniform wall temperature and uniform wall 
flux. Space limitations permit us to deal only with heat-transfer 
results here. Analytical details, plots of velocity and temper 
ature profiles, and tabulations of the new solutions will be giver 
elsewhere 

The previously reported exact solutions of the boundary-layer 
energy equation have been carried out for Prandtl Numbers 
greater than 0.6. For the Prandt! Number range of liquid metals 
only approximate formulations have been made, the results 
which are also presented here 

A sketch of the physical model showing pertinent coordinat 


and physical quantities is shown in Fig. 1 


RESULTS FOR UNIFORM WALL TEMPERATURE 


The heat-transfer results are reported in terms of the local heat 
transfer coefficient, 4, and local Nusselt Number, .Vu \ls 
involved are the Reynolds Number, Re,, and the Prandt! Nun 
ber, P» 
which are as follows: 


These several quantities have their usual definitions 


h g/(T.,, — T..), Nu hx/k, Re U.x/v,P Cpua/R 


where g is the local heat flux 

The heat-transfer results based on exact solutions of the 
boundary-layer equations are shown as the solid curve on Fig. - 
and are tabulated in Table 1 

Fig. 2 contains two other curves depicting the results of af 
proximate formulations, One of these is based on a Karma! 
Pohlhausen type approach 
polynomial approximation for temperature and velocity profiles 


i.e., integral formulation — plus 


is ust 


* The phrase ‘‘exact solutions of the boundary-layer equations 
to describe numerically computed similar solutions of the boundary-lay® 
differential equations his terminology is common in boundary-layé 


work —e.g.. Schlichting 
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i rence 2.) From that analysis, the local heat-transfer 


givell 


BS result { low Prandtl Number range is 


Vu,/(Re,Pr)'/? = 0.529/(1 + 0.82 Pr!/? (2 
The third curve of Fig. 2 is due to Morgan, Pipkin, and 
Warner Chey obtained a solution of the boundary-layer energy 
quation utilizing an approximation for the velocity profile. The 
jasis of their analysis is the fact that the viscous boundary layer 
‘< much thinner than the thermal boundary layer for very low 
Prandt] Numbers. So, as a first step, the velocity profile was 
simply taken as a uniform potential flow, and then an approxi 
nate correction was made for the viscous effects. The local 
eat-transfer results are expressed by 
Vu,/(Re,Pr)'/? = 0.564 — 0.547 Pr'/? (3 
\s the Prandtl Number decreases, it would be expected that 


the exact boundary-layer solution should approach more and 
more closely to Eq. (3), and this is verified by the curves of Fig 
29 The heat-transfer results of both approximate formulations 
agree very The 
generally no more than 5 per cent, with Eq. (3 


well with the exact solution deviations are 
showing even 
better agreement \s a final observation, it may be noted that 
the group Nu,/(Re,Pr)' 

Number in the liquid metal range 


a unique function of the Peclet Number (= 


2 is relatively independent of Prandtl 
In other words, the Nusselt 
Number is almost 
Re,P 


at 


RESULTS FOR UNIFORM HEAT FLUX 

The parameters defined in Eq. (1) are also used in presenting 
the results for the uniform-heat-flux case. The solid line of Fig. 3 
represents the findings of the exact boundary-layer solutions 
see Table 1 


m the approximate 


The other curves of Fig. 3 depict results based 


formulations which have already been 


described above. The equations from which these latter curves 
have been plotted are 


) 


Vu,/(Re,Pr)'/2 = 0.816/(1 + 1.064 Pr!/?) (4 


for the Karman-Pohlhausen formulation, and 


Nu,/(Re,Pr)'/? = 0.886 — 0.491 Pr!)/? (5 


for Morgan’s formulatior. 

Inspection of Fig. 3 shows that the approximate analyses give 
good results for the uniform-heat-flux case, although the spread 
umong the curves of this figure is slightly greater than on Fig. 2 
In particular, the deviation between Morgan’s results and the 
exact solution is both greater and of different sign for the uni 


form-heat-flux problem than for the uniform-wall-temperature 
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. 2. Local heat-transfer results for uniform-wall-temperature 
situation (hav = 2c). 
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Heat-Transfer Results From Exact Solutions 


Re, Pr) \/2 





P) Uniform g 
0.08 0.7170 
0.01 0.7756 
0.006 0.7969 


EXACT SOLUTION 
—-—— KARMAN- POHLHAUSEN METHOD 
—-—-—-— MORGAN VELOCITY APPROXIMATION (REF. 3) 
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lic. 3. Local heat-transfer results for uniform-heat-flux 
situation 
problem. The explanation may partly lie in the fact that the 
thermal boundary layer is thinner for the uniform-heat-flux case 
The Karman-Pohlhausen method appears to yield results of 


similar accuracy for the uniform-heat-flux case as for the uniform- 
wall-temperature cas 

3, it is again seen that the Nusselt Number is almost 
Number for the 


From Fig 
a unique function of the Peclet liquid metal 


range. 
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An Exact Solution of the Navier-Stokes 
Equations 


W. Chester 
Visiting Research Fellow, 
Pasadena, Calif. 


July 18, 1957 


California Institute of Technology, 


te NAVIER-STOKES EQUATIONS for a perfect gas with constant 
thermal conductivity and constant coefficients of viscosity 


may be written, using the conventional notation, as follows: 


V-(pV) = 0 
p(DV/Dt) + Vp = wV°V + (A w)V(V-V) 
pcp(DT/Dt) (Dp/Dt) = ® kVv?T 
p/bo = pT /poT 
If we assume that the Cartesian components of the velocity V = 
u(y), 0, O], that the pressure p = p(x, ¢), and that the density 


p = p, is constant, then the continuity equation is satisfied, and 
the momentum and energy equations reduce to 


[((Op/dt) + u(Op/dx)| = 
(4 ] 


Op/ox = p(d*u/dy?), 
u(du/dy)? 
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When 7 = 1-5, both these equations are satisfied by 
u = U[(1 — y?)/a?} 
T/T. = p/po = 1 — (2Up/a*p,.)(x — Ut) 
A similar result holds for the circular pipe provided » = 2. The 


last equation is then modified to 
T/T, = p/po = 1 


The solutions correspond respectively to steady flow with a 


- (4U'p/a*p,) (x lt 


parabolic velocity profile in an insulated two-dimensional channel 
of width 2a and in an insulated circular pipe of radius a. The 
velocity along the axis of the channel or pipe is U’, and the pres- 
sure remains constant for a particle of fluid moving with this 
The motion is produced by a linear pressure gradient; 
a fixed 


velocity. 
the pressure, itself, also increases linearly with time at 
point in the channel or pipe, as does the temperature. This could 
presumably be attained by connecting one end of the channel or 
pipe to a large reservoir into which heat is being added at a con 


stant rate. 


On One-Dimensional Heat Conduction With an 
Arbitrary Heating Rate 


George W. Sutton 
Research Engineer, Missile and Ordnance Systems Department, 
General Electric Company, Philadelphia, Pa. 


July 24, 1957 


ECENTLY THERE HAS BEEN INTEREST in one-dimensional heat 


conduction with heat transfer which varies arbitrarily 


with time at one boundary when the other boundary is insulated 
(see Fig. 1). This situation may exist, for example, on the skin 
of a rapidly accelerating aircraft or on a missile re-entering the 


earth’s atmosphere. For this case, the heat conduction equation 


1s 
zx ~— 1, = 0 (1) 
where 
X = 2/L 
r = at/L? 
t = time 
a = thermal diffusivity, k/pC 
7 = temperature 
The boundary conditions to Eq. (1) are 
Tx (0, 7) = 0 (2a) 
Tx(1, 7) = Lko'g(t) whent > 0 (2b) 
Tx(1, 7) = 0 when?t < 0 (2c) 


where q(t) is an arbitrary heating rate per unit area to the surface. 
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T (x,t) aT 


q(t) —__—_» 
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If g(t) may be represented by a polynomial—that is, 


q(t) = 


then by use of an obvious extension of the method of Laply 
transform given in reference 1 the temperature distributijoy 


given by 


> 


’ 


where 7(x, 0) = 7), a constant. In connection with a stud 


where s = 1, 2, design charts of the functions 
Zm(X, vr) = 2”"[(m — 1)/2]!r'/? } erfe [(27 + 1 — 
r=0 


Xr"/2/2] + i” erfe [(2r + 1 + Xr-"2)/2 


2s + 1 were prepared for values of s up to 5; that js 
m = 11. Thus, if the heating rate may be represented by 
polynomial up to the fifth degree, then the temperature dis 


where m 


tribution becomes very simply: 


Lk p A st*Zoss1 [(x/L), (at/L2)] + 7 


s=0 


T= 


With the use of the design charts of reference 2, the temperatur 
distribution is easily found. The accuracy of this solution i: 
generally limited to small times where the convergence of th 
series of Eq. (5) is rapid. 

We have, therefore, sought a solution to Eq. (1) which is muc 


simpler to calculate at large times. This solution is given by 


© 2(s+1—n) s+] 
T(X,7) = Do Al YO YS Bua X cy 70)" + 
s=0 m=0 n=0 
el 
| T(d, 0) dX + +» 8 cos (prX ee 7 
oO p=1 
where A,’ = L**'k-la~%A,. Eq. (7) is also an exact solutiont 


the problem. 

The coefficients B,,,,° 
tion equation and the boundary conditions for the corresponding 
The cosine series represents 


are chosen to satisfy the heat conduc 


s term in the heating rate of Eq. (3). 
the difference between the initial temperature distribution an 
Thus th 


the first two terms on the right side of Eq. (7) at ¢ = 0. 
coefficients C,, are the Fourier cosine coefficients of the function 


o/ x 2(s+1 n) 
T(X. 0) = | T(X,6)dX — DO As’ YS Bao X™ 6 
J 0 s=0 m=0 


The Fourier series decreases rapidly with time so that whe 
r > 0.5, its contribution is less than 2 per cent to 7(X, 7) at 
therefore may be neglected. Thus the solution is given simp! 
by a polynomial in XY and + 

The coefficients B,,,,“°) are simple to obtain 


requires that for any value of 


9 


First, Eq. (2 


B ) = () 


Next, substitution of Eq. (7) into Eq. (1) yields the following rt 


cursion formula for the coefficients: 


(m + 2) (m + LIBwrys, ni = (2 + 1L)Bm, nr” 


Eqs. (9) and (10) show that for all odd values of m 


Binn® =O (m odd 


Substitution of Eq. (7) into Eq. (2b) gives 
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TABLE I 
Values of the Coefficients B, 




























































n 
m 0 1 2 3 { 5 
) 0 0. 166667 1. QOOO0O 
» 0. 500000 
0 0.019444 0. 166667 0. 500000 
» 0.083333 0. 500000 
} 0.041667 
z 0 0.004101 0. OB8889 0. 166667 0.333333 
i » 0.019444 0. 166667 0. 500000 
: } 0.013889 0.083333 
6 0.002778 
2 0) 0.001594 0. 142857 0.058333 0. 166667 0. 250000 
» 0.007143 0.058333 0. 250000 0. 500000 
4 } 0.004861 —(). 041667 0. 125000 
6 0.001389 0.008333 
4 8 0.000175 
4 0 -~0.001559 0.010595 0. 126984 0. IS8S8S89 0. 166667 0. 200000 
» 0.005297 -0. 126984 0. 283333 0.333333 0.50000 
} 0.001058 0.047222 0). 166667 0). 166667 
6 0.001574 ~Q.011111 0.016667 
7 8 0.001984 0.000595 
10 0. 000007 
Bz, = 1/2 (12 Natural Frequencies of Rectangular Plates 
With Edges Elastically Restrained Against 
while for m # s, Rotation 
2(s+1 N 
m Br tha (13) C. V. Joga Rao and C. Lakshmi Kantham 
m=2 


Assistant Professor and Research Student, Respectively, Department 
of Aeronautical Engineering, Indian Institute of Science, 


a "sn ge Bo” Thi ni ; Bangalore, India 
cients may be calculated except a lis coefficient is ob- August 2, 1957 


By use of the value of By, ,“’, and Eqs. (10) and (13), all coeffi 


tained from the initial condition at ¢ = 0; that is 


2(s+1) 7 _—- with attachments to heavier members along the 

Bry X™dX = 0 (14 edges can be described as having edges elastically restrained 
sii against rotation, in many cases uniformly along each edge. At 
for each value of s the edges, setting slope a = 6.\J, when JM is the edge bending 
Values of the coefficients B,,,,"’ which were obtained by the use moment with 6 always positive, the elastic restraint can be an- 
of this method are tabulated in Table I for values of s up to 4 alytically defined en eS describing respectively, 
The coefficients for s = 2 given in Table I do not agree with those clamped and simply supported edges. In this note natural fre- 
of Eq. (I-31) in reference 3. However, it is noted that Eq. (1-31) quencies of such plates are calculated mainly following the 
does not satisfy the differential equation nomenclature of Dana ¥ a ; ; ; 
If the heating rate g(t) cannot conveniently be represented by a The differential equation for a uniform vibrating beam of stiff- 
polynomial, one may use in exact solution to Eq (1) with bound- ness E] and length / wah claaticalty panes rage Aner ~ 
characteristic functions. Taking the origin at the middle and 

. ‘4 . . 1 

setting 9, = VW p;2Ap/EI 1, where p,; = r" angular frequency A = 


: . , es area of cross section, and p = mass density, we have for sym- 
T(x, t) = g(t )G(x, t — th) dt; + 7 (15 
0 


ary conditions given by Eq. (2), which is 


where the Green’s function G(x, ¢ — 4), which has been obtained 








by solving Eq. (1) for g(t) = 6(t,), is given by: 
G(x,t — th) = (1/pCL) x 
| 3rd. MODE 
1+ 2 = (—1)" cos (nxx/L)e~ (Prta/L*) (th) 16 Sep erar  rerarar ay 
n=1 
2no. MODE 
Thus, the temperature distribution for an arbitrary heating rate T T 
may be calculated by simple quadrature of Eq. (15) 
| ao alas ae am hp oo wal 
40|—_+————_—_+—————_+ 4 =— ™ 4 + 4 
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metric vibrations with 


X, = [cosh (0,x//)/cosh (6,/2)| leos (0,%/1)/cos (6,/2)| — ( coth (9,/2 cot (6,/2) + 28( EI/1)6, 0 (4) 
with tanh (0,/2) + tan (9,/2) + 28 (EJ/l)o, = 0 It can be easily verified that, for the two limiting cases @ +g 
8 — oo, Eqs. (2) and (4) together yield the well-knowg 


equations cos 6, cosh #6, = Land sin #, = 0. Putting \ = 8(EI/pj 


and for antisymmetric vibrations 
these characteristic functions satisfy the following relationshipgs 


XY, = [sinh (0,x/l)/sinh (0,/2)] lsin (0,x/1)/sin (6,/2)| : get : 
for symmetric vibrations 


(1 nf - X,X,dx = Oforr # s 


= }d0,(tanh (6,/2 tan (6,/: + 20 +1} for 


©) /2 
if (d\,/dx) (dX 5/dx )dx = |40,70,2/(0,4 — 0,4)] }6,[tanh (6,/2 tan (6,/2)| — @,[tanh (@,/2) — tan (0,/2)|;} forr # s 
—i/2 


4 2) tan (6,/2) } forr = s 


= }(6,2/2) |tanh? (0,/2) + tan? (0,/2)| — 8,[tanh (6, 
PT/2 
rf (d2.X,/dx?) (d2.X,./dx?)dx = 8\0,76,2 for r + s 
1/2 
» 


= 6,4) d0,[tanh (6,/2) tan (6,/2 6\ + 1} forr =s 


for antisymmetric vibrations 


(1 nf ; a X,Y, dx = Oforr # s 


= }\0,[coth (0-/2) + cot (6,-/2)|] + 2A +4 1} forr = s 


; 2)] 9.\coth (8,/2 + cot (8,/2) { for 7 


el/2 
if - (dX,/dx) (dX ,/dx )dx [40,70,2/(0,' — 0,4)] ,@-leoth (6/2) + cot (6,/2) 2 


2/2) [eoth? (0,/2) + cot? (0,/2)] 6,/coth (6,/2) + cot (6, 2)]} forr =s 


vf ne (d?2X,/dx?) (d?X,/dx*)dx = 8\0,70,? for r # s 
» 


= 64 1)80,[coth (0,/2) + cot (0,/2 


Taking the expression for the deflection of the plate as wo = 
AmnXn(x) Y,(y), where X,,(x), Yn(y) are the character- ' ; 
u X . OV./OArs = (8A?/B) (b/a? 2. > A mn[0m70-27T ns 4+ 
m uN 


istic functions of vibrating beams with elastically restrained ends. 
The strain energy of bending V >, of the plate is given by? 
e strain energy of bending Vp, of the plate 1s given b) (a3/b3)6,,20,27 (2/8) (b “ik tp ee ee oe (18) 
m nN 


ea/2  (b/2 
Vp = (D 2) f f [ (O72) /Ox?) + 
a/2 —b;2 i : : ‘ 
and from consideration of kinetic energy 


(O7w/Oy?)|? dxdy (11) 
A OT /0A,, = abphp? >. > A, 
for the deflection is zero along the boundary; the energy stored rr = — 


in the constraints is given by 
abphp* =, 


Ve = (1/2) g Ma ds y - 


where / is the thickness of the plate. 


the contour being the boundary of the plate 
From principle of minimization, we have for a square plate 


Substituting for wo in terms of beam functions, 


OVp/0Ay = (D/a2)d>, YS Amn {(b/a)RmrT ns + > S> Ans (Canre + M1 — ¥*)Baore 


n 


(20) 
m nN 
(4°/6®)RusT mr + (4/B) (SmeSan + SrmSns) 3 where ¢? = a‘php?/D, D being the flexural rigidity of the plate 
The system of simultaneous equations (20) has been solved for 
Pa /2 @ for various values of \ for a square plate to obtain the lowest 
Rr = a f a/2 (d?.X',,/dx?) (d?X,/dx*) dx three frequencies. The convergence has been found to be rapi 
The results are presented in Fig. 1 


where 


b/2 
ms f , , (d?) n/dy*) (d?Y,/dy?) dy The results agree satisfactorily with the experimental value 
Js 


Hoppmann and 


for the fundamental frequency obtained by 
Greenspon.? This procedure can be modified to evaluate fre- 
quencies of rectangular plates with one value of elastic restraift 
factor along two parallel edges and with another value along the 
other two parallel edges. This is of interest in evaluation @f 
pase — - natural frequencies of rectangular wing panels 
a NX» (d?2X,/dx?) dx 
a/2 
‘i - REFERENCES 
= 5 a n(d?¥",/dy?) dy 
e b/2 1 Young, Dana, Vibration of Rectangular Plates by the Rits Method, J. Appl. 
Hence we write Mech., Vol. 17, pp. 448-453, December, 1950 
2 Wang, C. T., Applied Elasticitv, pp. 286, 287 


~ ‘ ments a 

(OV p/0A,.) = (D a)>, > ee oe Co., Inc., N.Y., 1953 
m ! * Hoppmann, W. H., II, and Greenspon, |] in Experimental Denice for 
Obtaining Elastic Rotational Constraint on Boundary of a Plate, Proc of the 


Second U.S. Nat. Congress of Applied Mechanics, pp. 187 
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Similarly we have 191. 1054 








